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Abstract 



A spin 1/2 system on a honeycomb lattice is studied. The interactions between nearest 
neighbors are of XX, YY or ZZ type, depending on the direction of the link; different types 
of interactions may differ in strength. The model is solved exactly by a reduction to free 
fermions in a static Z2 gauge field. A phase diagram in the parameter space is obtained. 
One of the phases has an energy gap and carries excitations that are Abelian anyons. 
The other phase is gapless, but acquires a gap in the presence of magnetic field. In the 
latter case excitations are non- Abelian anyons whose braiding rules coincide with those of 
conformal blocks for the Ising model. We also consider a general theory of free fermions 
with a gapped spectrum, which is characterized by a spectral Chern number v. The 
Abelian and non-Abelian phases of the original model correspond to v = and v = ±1, 
respectively. The anyonic properties of excitation depend on v mod 16, whereas v itself 
governs edge thermal transport. The paper also provides mathematical background on 
anyons as well as an elementary theory of Chern number for quasidiagonal matrices. 
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Comments to the contents: What is this paper about? 

Certainly, the main result of the paper is an exact solution of a particular two-dimensional 
quantum model. However, I was sitting on that result for too long, trying to perfect it, derive 
some properties of the model, and put them into a more general framework. Thus many 
ramifications have come along. Some of them stem from the desire to avoid the use of conformal 
field theory, which is more relevant to edge excitations rather than the bulk physics. This 
program has been partially successful, but some rudiments of conformal field theory (namely, 
the topological spin 9 a = e 2m ( h a- h a) anc [ ^he chiral central charge c_ = c — c) are still used. 

The paper is self-contained and provides an introduction into the subject. For most read- 
ers, a good strategy is to follow the exposition through the beginning of Sec. El Nonabelian 
anyons and take a glance at the rest of that section, where things become more technical. But 
the mathematically inclined reader may be interested in those details, as well as some of the 
appendices. I have tried to make the paper modular so that some parts of it can be understood 
without detailed reading of the other. This has caused some redundancy though. 

Appendix [El Algebraic theory of anyons is an elementary introduction into unitary modular 
categories, which generalizes the discussion in Sec. [HJ 

Appendix [Cl Quasidiagonal matrices is also mostly expository but some of the arguments 
may be new. It begins with a simplified treatment of "operator flow" and "noncommutative 
Chern number" (the latter has been used to prove the quantization of Hall conductivity in 
disordered systems PQ), but the main goal is to explain "unpaired Majorana modes", a certain 
parity phenomenon related to the Chern number. 

Appendix[D]on the chiral central charge and Appendix [F] on weak symmetry breaking contain 
some raw ideas that might eventually develop into interesting theories. 

Introduction 

Overview of the subject. Anyons are particles with unusual statistics (neither Bose nor 
Fermi), which can only occur in two dimensions. Quantum statistics may be understood as a 
special kind of interaction: when two particles interchange along some specified trajectories, the 
overall quantum state is multiplied by e %ip . In three dimensions, there is only one topologically 
distinct way to swap two particles. Two swaps are equivalent to the identity transformation, 
hence e lLp = ±1. On the contrary, in two dimensions the double swap corresponds to one 
particle making a full turn around the other; this process is topologically nontrivial. Therefore 
the exchange phase <p can, in principle, have any value — hence the name "anyon" . (However, 
a stability consideration requires that <p be a rational multiple of 2n.) Of course, the real 
question is whether such particles exist in nature or can be built somehow, but we will follow 
the historic path, approaching the problem from the mathematical end. 

The study of anyons was initiated by Wilczek [21 [3] in early 1980's. He proposed a simple 
but rather abstract model, which was based on (2 + l)-dimensional electrodynamics. This 
theory has integer electric charges and vortices carrying magnetic flux (which is a real number 
defined up to an integer). Considered separately, both kinds of particles are bosons. But when 
a charge q goes around a vortex v, it picks up the phase 2irqv, known as the Aharonov-Bohm 
phase. Thus, charges and vortices have nontrivial mutual statistics and therefore must be called 
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anyons when considered together. Moreover, composite objects (q, v) consisting of a charge and 
a vortex are anyons by themselves because they have nontrivial exchange phase (p( q v ) = lisqv. 

A general way to describe quantum statistics is to consider particle worldlines in the (2 + 1)- 
dimensional space-time. Such worldlines form a braid, therefore the statistics is characterized 
by a representation of the braid group. In the preceding discussion we assumed that braiding 
is characterized just by phase factors, i.e., that the representation is one- dimensional. The 
corresponding anyons are called Abelian. But one can also consider multidimensional represen- 
tations of the braid group; in this case the anyons are called non- Abelian. Actually, it may not 
be so important how the braid group acts, but the very existence of a multidimensional space 
associated with several particles is a key feature. Vectors in this space are quantum states that 
have almost the same energy (see discussion of topological quantum computation below). 

Historically, the theory of non- Abelian anyons emerged from conformal theory (CFT). How- 
ever, only topological and algebraic structure in CFT is relevant to anyons. Different pieces of 
this structure were discovered in a colossal work of many people, culminating in the paper by 
Moore and Seiberg |4j. Witten's work on quantum Chern-Simons theory [5] was also very influ- 
ential. A more abstract approach (based on local field theory) was developed by Fredenhagen, 
Rehren, and Schroer [6] and by Frohlich and Gabbiani [TJ. 

The most amazing thing about anyons is that they actually exist as excitations in some 
condensed matter systems. Such systems also have highly nontrivial ground states that are 
described as possessing topological order. The best studied example (both theoretically and 
experimentally) is the Laughlin state [8] in the fractional quantum Hall system at the filling 
factor v — 1/3. It carries Abelian anyons with exchange phase <p = ir/3 and electric charge 
±1/3. It is the fractional value of the charge that was predicted in original Laughlin's paper and 
confirmed by several methods, in particular by a shot noise measurement [9j [10] • The statistical 
phase is a subtler property which is deduced theoretically [Til [T2] ; a nontrivial experimental 
test has been performed recently using quasiparticle tunneling [T3] . 

A different kind of state is observed at the filling factor v — 5/2, though it is more fragile and 
less studied experimentally. There is much evidence suggesting that this system is described by 
a beautiful theory proposed by Moore and Read [HI [15]. The Moore- Read state admits non- 
Abelian anyons with charge ±1/4. If 2n such particles are present, the associated Hilbert space 
has dimension 2 n_1 . (The non- Abelian anyons studied in this paper have similar properties, 
though there is no electric charge.) 

The notion of anyons assumes that the underlying state has an energy gap (at least for 
topologically nontrivial quasiparticles). Otherwise excitations are not localizable and braiding 
may not be defined. Note that if all excitations are gapped, then all equal-time correlators 
decay exponentially with distance [IB]. 

An example of anyons in a spin- 1/2 system originates from the theory of resonating valence 
bond (RVB). The idea of RVB was put forward by Anderson [TTJ and used later as a model 
of the undoped insulating phase in high-T c cuprate superconductors [18]. Without electrically 
charged holes, the problem seems to be described adequately by a Heisenberg-like Hamiltonian, 
but its solution has proved very difficult. Several variants of an RVB state have been proposed, 
both gapless and gapped. Here we discuss a particular gapped RVB phase, namely the one 
which is realized on the triangular lattice [19] , but which apparently exists on the square lattice 
as well. This phase admits quasiparticles of four types: trivial excitations (such as spin waves), 
spinous (spin-1/2 fermions, which are conserved modulo 2), Z2-vortices (spinless bosons, also 
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a) b) 
Figure 1: A classical vortex (a) distorted by fluctuations (b). 

called visons), and spinon-vison composites [20]. The mutual statistics of spinons and visons is 
characterized by the Aharonov-Bohm factor —1, therefore the composite particles are bosons. 
Note that the relevance of this theory to cuprate superconductors is under debate. Senthil and 
Fisher proposed an interesting way to detect visons in these materials [21] , but the experiment 
gave a negative result [22J. However, some kind of RVB state is likely to realize in a different 
material, CS2CUCI4. This conclusion is drawn from neutron-scattering experiments that have 
shown the presence of spin-1/2 excitations [23] . 

Anyonic particles are best viewed as a kind of topological defects that reveal non-trivial 
properties of the ground state. Thus anyons carry some topological quantum numbers which 
make them stable: a single particle cannot be annihilated locally but only through the fusion 
with an antiparticle. An intuitive way to picture an anyon is to imagine a vortex in a medium 
with a local order parameter (see Fig. [1^,). Now suppose that quantum fluctuations are so 
strong that the order parameter is completely washed out and only the topology remains (see 
Fig.[Tb). Of course, that is only a rough illustration. It resembles the Kosterlitz-Thouless phase 
with power-law correlation decay, while in anyonic systems correlations decay exponentially due 
to the energy gap. 

A real example can be constructed with spins on the edges of a square lattice. Basis states 
of the spins are described by the variables Sj = ±1, which may be regarded as a Z 2 gauge field 
(i.e., "vector potential"), whereas the "magnetic field intensity" on plaquette p is given by 

w p = Y[ Sj . (1) 

j ^boundary (p) 

We say that there is a vortex on plaquette p if w p = — 1. Now we may define the vortex-free 
state: 

\^f) = c 2J |s), where s = (si, s^) (2) 

s: w p (s)=l 
for all p 

(c is a normalization factor). The state with a single vortex on a given plaquette is defined 
similarly. It is clear that the vortex can be detected by measuring the observable Y\j e i o~j for 
any enclosing path /, though no local order parameter exists. 

The state ([2]) can be represented as the ground state of the following Hamiltonian with 
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four-body interaction [23]: 

h = - j e a * - j ™ s b *> where As = n a i ' ^ = n 4 ( 3 ) 

vertices plaquettes star(s) boundary(p) 

Its elementary excitations are ^-charges with energy 2 J e and vortices with energy 2 J m . Certain 
essential features of this model are stable to small local perturbations (such as external magnetic 
field or Heisenberg interaction between neighboring spins). Note that the robust characteristic 
of excitations is not the energy or the property of being elementary, but rather supers election 
sector. It is defined as a class of states that can be transformed one to another by local operators. 
This particular model has the vacuum sector 1, the charge sector e, the vortex sector m, and 
a charge-vortex composite e. Particles of type e and m are bosons with nontrivial mutual 
statistics, whereas e is a fermion. Thus, the model represents a universality class of topological 
order — actually, the same class as RVB. 

Anyonic superselection sectors may or may not be linked to conventional quantum num- 
bers, like spin or electric charge. Most studies have been focused on the case where anyons 
carry fractional electric charge or half-integer spin. Such anyons are potentially easier to find 
experimentally because they contribute to collective effects (in particular, electric current) or 
have characteristic selection rules for spin-dependent scattering. Chargeless and spinless quasi- 
particles are generally harder to identify. But anyons, by virtue of their topological stability, 
must have some observable signatures. For example, anyons can be trapped by impurities and 
stay there for sufficiently long time, modifying the spectrum of local modes (magnons, excitons, 
etc.). However, effective methods to observe anyons are yet to be found. 

Thus, the hunt for anyons and topological order is a difficult endeavor. Why do we care? 
First, because these are conceptually important phenomena, breaking some paradigms. In 
particular, consider these principles (which work well and provide important guidance in many 
cases) : 

1. Conservation laws come from symmetries (by Noether's theorem or its quantum ana- 
logue); 

2. Symmetries are initially present in the Hamiltonian (or Lagrangian), but may be sponta- 
neously broken. 

Let us limit our discussion to the case of gauge symmetries and local conservation laws, which 
are described by fusion rules between superselection sectors. A profound understanding of the 
first principle and its underlying assumptions is due to Doplicher and Roberts [251 [26]. They 
proved that any consistent system of fusion rules for bosons is equivalent to the multiplication 
rules for irreducible representations of some compact group. Fermions also fit into this frame- 
work. However, anyonic fusion rules are not generally described by a group! As far as the 
second principle is concerned, topological order does not require any preexisting symmetry but 
leads to new conservation laws. Thus, the formation of topological order is exactly the opposite 
of symmetry breaking! 

Topological quantum computation. A more practical reason to look for anyons is their 
potential use in quantum computing. In Ref. [24J I suggested that topologically ordered states 
can serve as a physical analogue of error-correcting quantum codes. Thus, anyonic systems 



6 



provide a realization of quantum memory that is protected from decoherence. Some quantum 
gates can be implemented by braiding; this implementation is exact and does not require 
explicit error correction. Freedman, Larsen, and Wang [27] proved that for certain types of 
non-Abelian anyons braiding enables one to perform universal quantum computation. This 
scheme is usually referred to as topological quantum computation (TQC). 

Let us outline some basic principles of TQC. First, topologically ordered systems have 
degenerate ground states under certain circumstances. In particular, the existence of Abelian 
anyons implies the ground state degeneracy on the torus [28]. Indeed, consider a process in 
which a particle-antiparticle pair is created, one of the particles winds around the torus, and 
the pair is annihilated. Such a process corresponds to an operator acting on the ground state. 
If A and B are such operators corresponding to two basic loops on the torus, then ABA~ 1 B~ 1 
describes a process in which none of the particles effectively crosses the torus, but one of them 
winds around the other. If the Aharonov-Bohm phase is nontrivial, then A and B do not 
commute. Therefore they act on a multidimensional space. 

Actually, the degeneracy is not absolute but very precise. It is lifted due to virtual particle 
tunneling across the torus, but this process is exponentially suppressed. Therefore the distance 
between ground energy levels is proportional to exp(— L/£), where L the linear size of the torus 
and £ is some characteristic length, which is related to the gap in the excitation spectrum. 

In non-Abelian systems, degeneracy occurs even in the planar geometry when several anyons 
are localized in some places far apart from each other (it is this space of quantum states the 
braid group acts on). The underlying elementary property may be described as follows: two 
given non-Abelian particles can fuse in several ways (like multi-dimensional representations 
of a non-Abelian group). For example, the non-Abelian phase studied in this paper has the 
following fusion rules: 

e X s = 1, e X a = a, axa = l+e, 

where 1 is the vacuum sector, and e and a are some other superselection sectors. The last 
rule is especially interesting: it means that two cr-particles may either annihilate or fuse into 
an e-particle. But when the cr-particles stay apart, 1 and e correspond to two quantum states, 
\il>1 a ) and \ip° a )- These states are persistent. For example, if we create {ip™) by splitting an e 
into two cr's, wait some time, and fuse the cr-particles back, we will still get an e-particle. 

Here is a subtler property: the fusion states and \ip^ a ) are practically indistinguish- 

able and have almost the same energy. In fact, a natural process that "distinguishes" them by 
multiplying by different factors is tunneling of a virtual e-particle between the fixed cr-particles 
(which is possible since cr x e = a). However, e-particles are gapped, therefore this process 
is exponentially suppressed. Of course, this explanation depends on many details, but it is a 
general principle that different fusion states can only be distinguished by transporting a quasi- 
particle. Such processes are unlikely even in the presence of thermal bath and external noise, 
as long as the temperature and the noise frequency are much smaller than the gap. 

In the above example, the two-particle fusion states |V'i <T ) and \ip° u ) cannot form coherent 
superpositions because they belong to different superselection sectors (1 and e, resp.). In order 
to implement a qubit, one needs four cr-particles. A logical |0) is represented by the quantum 
state that is obtained by creating the pairs (1, 2) and (3, 4) from the vacuum (see Fig. [2]). A 
logical |1) is encoded by the complementary state |£ £ ): we first create a pair of e-particles, and 
then split each of them into a crcr-pair. Note that both states belong to the vacuum sector and 
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16) \Q \m) \ve) 

Figure 2: Four ways to initialize an anyonic qubit. 

therefore can form superpositions. Also shown in Fig. [2] are two alternative ways to initialize 
the qubit, \rji) and \r] £ ). The detailed analysis presented in Sees. 18.41 and I8"31 implies that 

b7i> = ^(l6> + l&>), ta> = ^(l&> -&>)• 

Therefore we can perform the following gedanken experiment. We create the state and then 
measure the qubit in the { \rji), |?7 e )}-basis by fusing the pairs (1, 3) and (2, 4). With probability 
1/2 both pairs annihilate, and with probability 1/2 we get two e-particles. One can also think 
of a simple robustness test for quantum states: if there is no decoherence, then both and 
\r]i) are persistent. 

As already mentioned, braiding is described by operators that are exact (up to virtual 
quasiparticle tunneling). Indeed, the operators of counterclockwise exchange between two par- 
ticles (.R-matrices) are related to the fusion rules by so-called hexagon equations and pentagon 
equation. We will see on concrete exapmles that these equations have only a finite number of 
solutions and therefore do not admit small deformations. In general, it is a nontrivial theorem 
known as Ocneanu rigidity [2H1 ED], see Sec. IE. 61 

Thus, we have all essential elements of a quantum computer implemented in a robust fashion: 
an initial state is made by creating pairs and/or by splitting particles, unitary gates are realized 
by braiding, and measurements are performed by fusion. This "purely topological" scheme is 
universal for sufficiently complicated phases such as the k = 3 parafermion state [31], lattice 
models based on some finite groups (e.g. 5*5 [24J, A 5 [321 [33] and S3 [M]), and double Chern- 
Simons models [33, [3B1 EZ]. Unfortunately, the model studied in this paper is not universal 
in this sense. One can, however, combine a topologically protected quantum memory with a 
nontopological realization of gates (using explicit error correction). Note that some weak form 
of topological protection is possible even in one- dimensional Josephson junction arrays [38J, 
which is due to the build-in U(l)-symmetry. Several other schemes of Josephson junction- 
based topological quantum memory have been proposed recently [311 SHI HI]- 

Unlike many other quantum computation proposals, TQC should not have serious scalability 
issues. What is usually considered an initial step, i.e., implementing a single gate, may actually 
be close to the solution of the whole problem. It is an extremely challenging task, though. It 
demands the ability to control individual quasiparticles, which is beyond the reach of present 
technology. One should however keep in mind that the ultimate goal is to build a practical 
quantum computer, which will contain at least a few hundred logical qubits and involve error- 
correcting coding: either in software (with considerable overhead) or by topological protection 
or maybe by some other means. At any rate, that is a task for the technology of the future. 
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But for the meantime, finding and studying topological phases seems to be a very reasonable 
goal, also attractive from the fundamental science perspective. 

Comparison with earlier work and a summary of the results. In this paper we study 
a particular exactly solvable spin model on a two-dimensional lattice. It only involves two-body 
interactions and therefore is simpler than Hamiltonian ([3]) considered in [24] . but the solution 
is less trivial. It is not clear how to realize this model in solid state, but an optical lattice 
implementation has been proposed [12]. 

The model has two phases (denoted by A and B) which occur at different values of param- 
eters. The exact solution is obtained by a reduction to free real fermions. Thus quasiparticles 
in the system may be characterized as fermions and Z 2 -vortices. Vortices and fermions interact 
by the Aharonov-Bohm factor equal to —1. In phase A the fermions have an energy gap, and 
the vortices are bosons that fall into two distinct superselection sectors. (Interestingly enough, 
the two types of vortices have identical physical properties and are related to each other by a 
lattice translation.) The overall particle classification, fusion rules, and statistics are the same 
as in the model (T5]) or RVB. In phase B the fermions are gapless and there is only one type of 
vortices with undefined statistics. Adding a magnetic field to the Hamiltonian opens a gap in 
the fermionic spectrum, and the vortices become non-Abelian anyons. The difference between 
the vortex statistics in phase A and phase B with the magnetic field may be attributed to 
different topology of fermionic pairing. 

Topological properties of Fermi-systems were first studied in the theory of integer quantum 
Hall effect [43, 44J. Let us outline the main result. To begin with, the Hall conductivity of 
noninteracting electrons in a periodic potential (e.g., in the Hofstadter model with m/n flux 
quanta per plaquette) is expressed in terms of a single-electron Hamiltonian in the Fourier basis. 
Such a Hamiltonian is an n x n matrix that depends on the momentum q. For each value of q 
one can define a subspace £(q) C C n that is associated with negative-energy states, i.e., ones 
that are occupied by electrons. Thus, a vector bundle over the momentum space is defined. 
The quantized Hall conductivity is proportional to the Chern number of this bundle. Bellissard 
at al [1] have generalized this theory to disordered systems by using a powerful mathematical 
theory called noncommutative geometry [45J. 

Even more interesting topological phenomena occur when the number of particles is not 
conserved (due to the presence of terms like ajaj., as in the mean-field description of supercon- 
ductors). In this case the single-electron Hamiltonian is replaced by a more general object, the 
Bogolyubov-Nambu matrix. It also has an associated Chern number u, which is twice the num- 
ber defined above when the previous definition is applicable. But in general v is an arbitrary 
integer. The first physical example of this kind, the 3 He-A film, was studied by Volovik [46J. 
Volovik and Yakovenko |47j showed that the Chern number in this system determines the statis- 
tics of solitons. More recently, Read and Green [45] considered BCS pairing of spinless particles 
with angular momentum I = —1. They identified a "strong pairing phase" with zero Chern 
number and a "weak pairing phase" with v — 1. The latter is closely related to the Moore- Read 
state and has non-Abelian vortices and chiral edge modes. 

In the present paper, these results are generalized to an arbitrary Fermi-system described 
by a quadratic Hamiltonian on a two-dimensional lattice. We show that Z 2 -vortices are Abelian 
particles when the Chern number v is even and non-Abelian anyons when v is odd. The non- 
Abelian statistics is due to unpaired Majorana modes associated with vortices. Our method 
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relies on a quasidiagonal matrix formalism (see Appendix [C]), which is similar to, but more 
elementary than, noncommutative geometry. It can also be applied to disordered systems. 

Furthermore, we find that there are actually 16 (8 Abelian and 8 non-Abelian) types of 
vortex-fermion statistics, which correspond to different values of v mod 16. Only three of them 
(for v — 0, ±1) are realized in the original spin model. We give a complete algebraic description 
of all 16 cases, see tables on pages 1301 HU andl42l 



1 The model 



We study a spin-1/2 system in which spins are located at the vertices of a honeycomb lattice, 
see Fig. [3h, This lattice consists of two equivalent simple sublattices, referred to as "even" and 
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a) b) 
Figure 3: Three types of links in the honeycomb lattice. 

"odd" (they are shown by empty and full circles in the figure). A unit cell of the lattice contains 
one vertex of each kind. Links are divided into three types, depending on their direction (see 
Fig. [3b); we call them "x-links", u y- links", and "z-links". The Hamiltonian is as follows: 



y-Ymks 



(4) 



z-links 



-links 



where J x , J y , J z are model parameters. 

Let us introduce a special notation for the individual terms in the Hamiltonian: 

{crjcr^, if (j, k) is an x-link; 
a?a y k , if (j, k) is an y-link; 
crjcr^, if (j, k) is an z-link. 

Remarkably, all operators Kjf. commute with the following operators W p , which are associated 
to lattice plaquettes (i.e., hexagons): 



(5) 




2 4' 
P 



W n = a x x a\alalalal = K l2 K 23 K u K A5 K 56 K 6l . 




(6) 
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Here p is a label of the plaquette. Note that different operators W p commute with each other. 

Thus Hamiltonian (j3j) has the set of "integrals of motion" W p , which greatly simplifies the 
problem. In order to find eigenstates of the Hamiltonian, we first divide the total Hilbert space 
C into sectors — eigenspaces of W p , which are also invariant subspaces of H. This can be 
written as follows: 

^ = ^M)l,...,W ro ) (7) 

Wi,...,W m 

where m is the number of plaquettes. Each operator W p has eigenvalues +1 and —1, therefore 
each sector corresponds to a choice of w p — ±1 for each plaquette p. Then we need to solve for 
the eigenvalues of the Hamiltonian restricted to a particular sector C Wl) .. Mm . 

The honeycomb lattice has 1/2 plaquette per vertex, therefore m ~ n/2, where n is the 
number of vertices. It follows that the dimension of each sector is ~ 2 n /2 m ~ 2 n / 2 (we will 
in fact see that all these dimensions are equal). Thus splitting into sectors does not solve the 
problem yet. Fortunately, it turns out that the degrees of freedom within each sector can be 
described as real (Majorana) fermions, and the restricted Hamiltonian is simply a quadratic 
form in Majorana operators. This makes an exact solution possible. 



2 Representing spins by Majorana operators 
2.1 A general spin-fermion transformation 

Let us remind the reader some general formalism pertaining to Fermi systems. A system 
with n fermionic modes is usually described by the annihilation and creation operators a/-, at 
(k — 1, . . . , n). Instead, one can use their linear combinations, 

C2fc-i — Ofc + a, , C2k — ; , 

i 

which are called Majorana operators. The operators Cj (j = 1, . . . , 2n) are Hermitian and obey 
the following relations: 

c) = 1, CjCi = -cicj if j ^ I. (8) 

Note that all operators Cj can be treated on equal basis. 

We now describe a representation of a spin by two fermionic modes, i.e., by four Majorana 
operators. Let us denote these operators by b x , b y , b z , and c (instead of ci, C2, C3, C4). The 
Majorana operators act on the 4-dimensional Fock space_^M, whereas the Hilbert space of a 
spin is identified with a two-dimensional subspace M. C M. defined by this condition: 

\i) EM if and only if D|f) = \0, where D = b x b y b z c. (9) 

We call M. and Ai the physical subspace and the extended space, respectively; the operator D 
may be thought of as a gauge transformation for the group Z 2 . 

The Pauli operator o~ x , a y , a z can be represented by some operators a x , a y , a z acting on 
the extended space. Such a representation must satisfy two conditions: (1) a x , a y , a z preserve 
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the subspace M.; (2) when restricted to M., the operators a x , a y , a z obey the same algebraic 
relations as a x , a y , a z . We will use the following particular representation: 

b z 

x * y a x = ib x c, a y = ib y c, a z = ib z c. (10) 



(We have associated the Majorana operators with four dots for a reason that will be clear later.) 
This representation is correct since a a (a = x, y, z) commutes with D (so that A4 is preserved), 
(a Q )t = a a , (a a ) 2 = 1, and 

a x a y a z = ib x b y b z c = iD. 

The last equation is consistent with the formula a x a y a z = i because D acts as the identity 
operator on the subspace M.. 

A multi-spin system is described by four Majorana operators per spin. The corresponding 
operators a aj , D? and the physical subspace C C £ are defined as follows: 

a a i = ibfcj, Dj = t>jb ! 'bjCj: 

' ' ' " (11) 
\£)e£ if and only if Dj\£) = |£) for all j. 

Any spin Hamiltonian H{a'j} can be replaced by the fermionic Hamiltonian iJ{6",Cj} = 
H{a aj } the action of which is restricted to the physical subspace. (The resulting Hamilto- 
nian H is rather special; in particular, it commutes with the operators Dj.) 

Remark 2.1. The substitution cr" a? = ib^Cj is gauge-equivalent to a more familiar one 
(see [H] and references therein): 

rr'/ - i.e., a x ^ -ibfi*, o] ■ * ,\>M]. a z i-> -ibffi. (12) 

Thus one can represent a spin by only 3 Majorana operators without imposing gauge constraints. 
However, this is not sufficient for our purposes. 



2.2 Application to the concrete model 



Let us apply the general procedure to the the spin Hamiltonian (jlj). Each term Kj k = <J®a% 

becomes Kjk = (ib'jCj)(ib'^Ck) = —i (ib'jb^) CjC^. The operator in parentheses, Ujk = ib^bf, is 
Hermitian; we associated it with the link (j, k). (The index a takes values x, y or z depending 
on the direction of the link, i.e., a = (Xjk-) Thus we get: 



H — - AjkCjC k} A jk — < 



2J a . k Ujk if j and k are connected, 
otherwise, 



hk I uuuuw**, (13) 

u jk = ib? k b a *. 

Note that each pair of connected sites is counted twice, and Ukj — — %fc- The structure of this 
Hamiltonian is shown in Fig. |H 
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Figure 4: Graphic representation of Hamiltonian ([131) . 



Remarkably, the operators Ujk commute with the Hamiltonian and with each other. There- 
fore the Hilbert space £ splits into common eigenspaces of Ujk, which are indexed by the 
corresponding eigenvalues Ujk = ±1. Similarly to ([7j) we may write C = where u 

stands for the collection of all Ujk- The restriction of Hamiltonian ({TBI to the subspace C u is 
obtained by "removing hats", i.e., replacing operators by numbers. This procedure results in 
the Hamiltonian H u = | . k Aj k CjCk, which corresponds to free fermions. The ground state of 

H u can be found exactly; let us denote it by \^? u )- 

Note, however, that the subspace C u is not gauge-invariant: applying the gauge operator 
Dj changes the values of Ujk on the links connecting the vertex j with three adjacent vertices k. 
Therefore the state \^f u ) does not belong to the physical subspace. To obtain a physical space 
we must symmetrize over all gauge transformations. Specifically, we construct the following 
state: 

i*«> = n i*«> e c - ( 14 ) 

Here w denotes the equivalence class of u under the gauge transformations. For the planar 
lattice (but not on the torus) w is characterized by numbers w p = ±1 defined as products 
of Ujk around hexagons. To avoid ambiguity (due to the relation Ukj = — we choose a 
particular direction for each link: 

w p = I^l u i k (j G even sublattice, k G odd sublattice) . (15) 

(j',fc)£boundary(p) 

The corresponding operator W p = Yl Ujk commutes with the gauge transformations as well as 
the Hamiltonian. The restriction of this operator to the physical subspace coincides with the 
integral of motion W p defined earlier (see (jSJ)). 

Notation change: From now on, we will not make distinction between operators acting in 
the extended space and their restrictions to the physical subspace e.g., W p versus W p . The tilde 
mark will be used for other purposes. 
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2.3 Path and loop operators 



One may think of the variables 11 kj clS db Z 2 gauge field. The number w v is interpreted as the 
magnetic flux through the plaquette p. If w p = —1, we say that the plaquette carries a vortex. 

The product of Uj k along an arbitrary path corresponds to the transfer of a fermion between 
the initial and the final point. However, this product is not gauge- invariant. One can define an 
invariant fermionic path operator in terms of spins or in terms of fermions: 

W(j ,...,j n ) = K jnj , 

where Kj k is given by ((Sj). If the path is closed, i.e., j n = j , the factors c n and c cancel each 
other. In this case the path operator is called the Wilson loop; it generalizes the notion of 
magnetic flux. 

On the honeycomb lattice all loops have even length, and formula ( fT6l) agrees with the sign 
convention based on the partition into the even and odd sublattice. However, the spin model 
can be generalized to any trivalent graph, in which case the loop length is arbitrary. For an odd 
loop I the operator W(l) has eigenvalues W\ = ±i. This is not just an artifact of the definition: 
odd loops are special in that they cause spontaneous breaking of the time-reversal symmetry. 

The time-reversal operator is a conjugate-linear unitary operator T such that 

TaJT- 1 = -of, TbjT- 1 = bj, TcjT~ l = Cj . (17) 

(The first equation is a physical requirement; the other two represent the action of T in the 
extended space.) Therefore T commutes with the Hamiltonian (Jlj) and the Wilson loop. Mul- 
tiplying the equation = wi\^) by T, we get WiT|\P) = wfT\^f). Thus the time-reversal 
operator changes Wi to w^. For a bipartite graph u> ; * = wi for all loops, therefore fixing the 
variables Wi does not break the time-reversal symmetry. On the contrary, for a similar model 
on a non-bipartite graph (e.g., a lattice containing triangles) the operator T does not preserve 
the field configuration, which is defined by the values of wi on all loops. But T is a symmetry 
operator, therefore all Hamiltonian eigenstates are (at least) two-fold degenerate. 



n 



-iu 



Jsjs-l 



C n CQ, 



(16) 



3 Quadratic Hamiltonians 

In the previous section we transformed the spin model (jlj) to a quadratic fermionic Hamiltonian 
of the general form 

H(A)= l -Y,A jk c 3 c k , (18) 

where A is a real skew-symmetric matrix of size n = 2m. Let us briefly state some general 
properties of such Hamiltonians and fix the terminology. 

First, we comment on the normalization factor 1/4 in Eq. ( |T8l) . It is chosen so that 

[-iH(A),-iH(B)\ = -iH([A,B\). (19) 

Thus the Lie algebra of quadratic operators —iH(A) (acting on the 2 m -dimensional Fock space) 
is identified with so (2m). Operators of the form e ~ lH ( A ) constitute the Lie group Spin(2m). 
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The center of this group consists of phase factors ±1 (e.g., e 7rClC2 = —1). The quotient group 
Spin (2m) /{+1, —1} = SO(2m) describes the action of e~ lH<yA ^> on Majorana operators by con- 
jugation: 

'^^QjkCj, where Q = e A (20) 



(Q is a real orthogonal matrix with determinant 1). 

Note that the sum in Eq. ([2"0l corresponds to the multiplication of the row vector (ci, . . . , C2 m ) 
by the matrix Q. On the other hand, when we consider a linear combination of Majorana op- 
erators, 

F{x)=Y< x i c v ( 21 ) 

3 

we prefer to view x as a column vector. If the coefficients Xj are real, we call F(x) (or x itself) 
a Majorana mode. 

To find the ground state of the Hamiltonian (ITS]) , one needs to reduce it to a canonical form 



-^canonical — p ^^^^'kK ~ e k( a l a k ~ |) , ^ 0, (22) 

k=l k=l 

where b' k , b k are normal modes, and a\ = \(b' k — ib k ), a k = \{b' k + ib k ) are the corresponding cre- 
ation and annihilation operators. The ground state of -^canonical is characterized by the condition 
Ofcl^) =0 for all k. The reduction to the canonical form is achieved by the transformation 

(P'n K, ■■■,b' m , 6* ) = (ci, c 2 , . . . , c 2m _i, c 2m ) Q, Q E 0(2m), (23) 



such that 



A = Q 



/ £1 

-ei 



V 







(24) 



The numbers constitute the spectrum of the Hermitian matrix iA, whereas odd (even) 
columns of Q are equal to the real (resp., imaginary) part of the eigenvectors. The ground 
state of the Majorana system has energy 



E 



9T, 



£k 



k=l 



4 



Tr \iA\ 



(25) 



where the function | ■ | acts on the eigenvalues, the eigenvectors being fixed. (In fact, any 
function of a real variable can be applied to Hermitian matrices.) 

Note that different quadratic Hamiltonians may give rise to the same ground state. The 
latter actually depends on 

/0 1 \ 

-1 



B = — isgn(iA) = Q 



1 
-1 0/ 



Q 1 



(26) 
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(We assume that A is not degenerate.) The matrix B is real skew-symmetric and satisfies 
B 2 = — 1. It determines the ground state through the condition 

^ PjkCj 1^) = for all k, where P jk = -(8 jk - iB jk ). (27) 

j 

Loosely speaking, B corresponds to a pairing^ between Majorana modes. The operators b' = 
F(x') and b" = F(x") are paired if x" = ±Bx' . 

The matrix P in Eq. fl27|) is called the spectral projector. It projects the 2m-dimensional 
complex space C 2m onto the m-dimensional subspace L spanned by the eigenvectors of %A 
corresponding to negative eigenvalues. For any vector z G L the corresponding operator F(z) 
annihilates the ground state, so we may call L the space of annihilation operators. Note that if 
z,z' G L then ^ . ZjZj = 0. The choice of an m-dimensional subspace L C C 2m satisfying this 
condition is equivalent to the choice of matrix B. 

The ground state of a quadratic Hamiltonian can also be characterized by correlation func- 
tions. The second-order correlator is (^/\cjCk\^) = 2Pjy; higher-order correlators can be found 
using Wick's formula. 

4 The spectrum of fermions and the phase diagram 

We now study the system of Majorana fermions on the honeycomb lattice. It is described 
by the quadratic Hamiltonian H u = H(A), where Aj k = 2J a . k Uj k , Uj k = ±1. Although the 
Hamiltonian is parametrized by Uj k , the corresponding gauge-invariant state (or the state of 
the spin system) actually depends on the variables w p , see (fT51) . 

First, we remark that the global ground state energy does not depend on the signs of the 
exchange constants J x , J y , J z since changing the signs can be compensated by changing the 
corresponding variables Uj k . We further notice that the ground state energy for H u does not 
depend on these signs even if u is fixed. Suppose, for instance, that we replace J z by — J z . 
Such a change is equivalent to altering Uj k for all z-links. But the gauge-invariant quantities 
w p remain constant, so we may apply a gauge transformation that returns Uj k to their original 
values. The net effect is that the Majorana operators at some sites are transformed clS Cj I ^ Cj . 
Specifically, the transformation acts on the set of sites Q z that lie in the shaded area in the 
picture below. In terms of spins, this action is induced by the unitary operator 



Rz = n ° z i ( 2g ) 

So, for the purpose of finding the ground state energy and the excitation spectrum, the signs 
of exchange constants do not matter (but other physical quantities may depend on them). 

^^Note the nonstandard use of terminology: in our sense, electron half-modes in an insulator are paired as 
well as in a superconductor. The only difference is that the insulating pairing commutes with the number 
of electrons while the superconducting one doesn't. This distinction is irrelevant to our model because the 
Hamiltonian doesn't preserve any integral charge, and the number of fermions is only conserved modulo 2. 
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The most interesting choice of Ujk is the one that minimizes the ground state energy. It turns 
out that the energy minimum is achieved by the vortex- free field configuration, i.e., w p — 1 



for all plaquettes p. This statement follows from a beautiful theorem proved by Lieb [50J. 
(Not knowing about Lieb's result, I did some numerical study suggesting the same answer, see 
Appendix |Al) Thus we may assume that Ujk = 1 for all links (j,k), where j belongs to the 
even sublattice, and k belongs to the odd sublattice. This field configuration (denoted by uf k d ) 
possesses a translational symmetry, therefore the fermionic spectrum can be found analytically 
using the Fourier transform. 

The general procedure is as follows. Let us represent the site index j as (s, A), where s refers 
to a unit cell, and A to a position type inside the cell (we choose the unit cell as shown in the 
figure accompanying Eq. (1321)). The Hamiltonian becomes H = (i/4) ^ s A t A^t^AC^, where 
A s \ t tv actually depends on A, /i, and t — s. Then we pass to the momentum representation: 



e- 4{q ' rs) c sA , 



(29) 
(30) 



where N is the total number of the unit cells. (Here and on, operators in the momentum 
representation are marked with tilde.) Note that a} 



q,A 



-q,A aIld a p,M a q,A + a q,A a p,M 



^pq^A- 



The spectrum e(q) is given by the eigenvalues of the matrix iA(q). One may call it a "double 
spectrum" because of its redundancy: e(— q) = — e(— q). The "single spectrum" can be obtained 
by taking only positive eigenvalues (if none of the eigenvalues is zero). 
We now apply this procedure to the concrete Hamiltonian 



vortex- 



free ^ ^ ^ AjfcCjCk] 



Ajk — 2J a uf k d . 



(31) 



We choose a basis (ni, n 2 ) of the translation group and obtain the following result: 
unit cell 

/(q) = 2(j x e i ^ + Jye^) + J z ), 




(32) 



where nx = (|, ^), n 2 = (— |, ^) in the standard xy-coordinates. 

An important property of the spectrum is whether it is gapless, i.e., whether e(q) is zero 
for some q. The equation J x e l ^' ni ^ + J 2/ e^ q ' n2 - ) + J z = has solutions if and only if \J, 
\J Z \ satisfy the triangle inequalities: 



r | ■ | Jy \ i 



\Jx\ ^ \ Jv\ \Jz 



\Jv\ ^ \ Jx\ \ Jz 



\Jz\ ^ \ Jx\ ~t~ \ Jy\ 



(33) 



If the inequalities are strict ("<" instead of '%"), there are exactly 2 solutions: q = ±q*. 
The region defined by inequalities ( 1331) is marked by B in Fig. [5j this phase is gapless. The 
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Figure 5: Phase diagram of the model. The triangle is the section of the positive octant 
(J x , J y , J z ^ 0) by the plane J x + J y + J z = 1. The diagrams for the other octants are similar. 

gapped phases, A x , A y , and A z , are algebraically distinct, though related to each other by 
rotational symmetry. They differ in the way lattice translations act on anyonic states (see 
Section l5\2l) . Therefore a continuous transition from one gapped phase to another is impossible, 
even if we introduce new terms in the Hamiltonian. On the other hand, the 8 copies of each 
phase (corresponding to different sign combinations of J x , J y , J z ) have the same translational 
properties. It is unknown whether the 8 copies of the gapless phase are algebraically different. 

We now consider the zeros of the spectrum that exist in the gapless phase. The momentum q 
is defined modulo the reciprocal lattice, i.e., it belongs to a torus. We represent the momentum 
space by the parallelogram spanned by (qi, q 2 ) — the basis dual to (n 1; n 2 ). In the symmetric 
case (J x = J y = J z ) the zeros of the spectrum are given by 



If \J X \ and \ J y \ decrease while \ J Z \ remains constant, q* and — q* move toward each other (within 
the parallelogram) until they fuse and disappear. This happens when \J X \ + \J y \ = \J Z \. The 
points q* and — q* can also effectively fuse at opposite sides of the parallelogram. (Note that 
the equation q* = — q* has three nonzero solutions on the torus). 

At the points ±q* the spectrum has conic singularities (assuming that q* ^ — q*): 




q* = §qi + §q2 (mod qi, q 2 ) 
q* = fqi + |q2 (mod qi, q 2 ) 



(34) 




where 5q = q — q* or <5q = q + q* . 




(35) 
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5 Properties of the gapped phases 



In a gapped phase, spin correlations decay exponentially with distance, therefore spatially 
separated quasiparticles cannot interact directly That is, a small displacement or another 
local action on one particle does not influence the other. However, the particles can interact 
topologically if they move around each other. This phenomenon is described by braiding rules. 
(We refer to braids that are formed by particle worldlines in the 3-dimensional space-time.) In 
our case the particles are vortices and fermions. When a fermion moves around a vortex, the 
overall quantum state is multiplied by —1. As mentioned in the introduction, such particles 
(with braiding characterized simply by phase factors) are called Abelian anyons. 

The description of anyons begins with identifying supers election sectors, i.e., excitation types 
defined up to local operations. (An "excitation" is assumed to be localized in space, but it may 
have uncertain energy or be composed of several unbound particles.) The trivial superselection 
sector is that of the vacuum; it also contains all excitations that can be obtained from the 
vacuum by the action of local operators. 

At first sight, each gapped phase in our model has three superselection sectors: a fermion, 
a vortex and the vacuum. However, we will see that there are actually two types of vortices 
that live on different subsets of plaquettes. They have the same energy and other physical 
characteristics, yet they belong to different superselection sectors: to transform one type of 
vortex into the other one has to create or annihilate a fermion. 

To understand the particle types and other algebraic properties of the gapped phases, we 
will map our model to an already known one [24|. Let us focus on the phase A z , which occurs 
when | J x \ + | J y \ < \J Z \. Since we are only interested in discrete characteristics of the phase, we 
may set | J x \, \ J y \ -C \J Z \ and apply the perturbation theory. 



5.1 Perturbation theory study 

The Hamiltonian is H = H + V, where H is the main part and V is the perturbation: 

h = —j z j2 v = -j, £ °n - j vH ° y A 

z-links x-links j/-links 

Let us assume that J z > (the opposite case is studied analogously). 

We first set J x = J y = and find the ground state. It is highly degenerate: each two spins 
connected by a z-link are aligned (ft or ID) but their common direction is not fixed. We 
regard each such pair as an effective spin. The transition from physical spins to effective spins 
is shown in Fig. Ek,b. The ground state energy is E = —NJ Z , where N is the number of unit 
cells, i.e., half the number of spins. 

Our goal is to find an effective Hamiltonian that would act in the space of effective spins 
£ e ff- One way to solve the problem is to choose a basis in £ eff and compute the matrix elements 

w£\ b) = { a\v\ b) , < |i,s>|6) = £'ra^m etc . 



J 



However, we will use the more general Green function formalism. 

Let T : C eS — > £ be the embedding that maps the effective Hilbert space onto the ground 
subspace of H . The map T simply doubles each spin: T|m) = \mm), where m = j or 
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b) 




c) 



Figure 6: Reduction of the model. Strong links in the original model (a) become effective spins 
(b), which are associated with the links of a new lattice (c). 



m = {. The eigenvalues of the "effective Hamiltonian" (if one exists) are supposed to be the 
energy levels of H that originate from ground states of Hq. These levels can be unambiguously 
defined as poles of the Green function G(E) = T'(E — H)~ 1 T, which is an operator acting 
on £ e ff and depending on the parameter E. The Green function is conventionally expressed as 
(E — Eq — £(-£?) ) , where £(.£>) is called self-energy, so the energy levels in question are the 
values of E for which the operator E — Eq — H{E) is degenerate. Neglecting the the dependence 
of S(-E') on E (for E m Eq), we define the effective Hamiltonian as H e s = Eq + T,(Eq). 

The self-energy is computed by the standard method. Let G' (E) = ((E — Ho)" 1 ) be the 
unperturbed Green function for exited states of H . The ' sign indicates that the operator 
[[E — Hq)^ 1 ) acts on excited states in the natural way but vanishes on ground states. Then 



£(£) = T f (y + VG' {E)V + VG' {E)VG' {E)V + ■ ■ ■ ) T, 



(36) 



We set E = E and compute H e g = E + £ in the zeroth order (H^ = Eq), first order, second 
order, and on, until we find a nonconstant termj§ The calculation follows. 



1. H$ = XVT = 0. 



2. H. 



(2) 



eff 



72 72 



4J, 



J Q£ l J U 

—N -. Indeed, consider the action 

4J Z 



z-links y-links 

of the second V in the expression T^VG' VT . Each term a^a^ or o v -o\ flips two spins, 
increasing the energy by 4J 2 . The other V must flip them back. 



3 H {3) 

°- ""eff 



r*vG' vG' vr = o. 



4. H. 



(4) 
eff 



J 2 J 2 

T^VG'qVG'qVG' g VT = const - 



^^<5 P , where Q P = {W p ) e g is the effective 



spin representation of the operator ([6]). The factor ^ j s obtained by summing 24 terms 



2 Higher order terms may be less significant than the dependence of S(-E) on E in the first nonconstant term. 
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each of which corresponds to flipping 4 spin pairs in a particular order: 



11-11 

— = 8- — + 8 + 8 . 

16 64 64 128 

The above arguments can easily be adapted to the case J z < 0. Now we have T : \]) i— > 
|ti)> II) l— *■ lit)- The result turns out to be the same, with J z replaced by \J Z \. 
Thus the effective Hamiltonian has the form 

J 2 J 2 

HeS = ~ IQ\J V \1 y^^P' Q'P = CT fcft(p) <T fight(p)° r up(p) Cr down(p) (37) 

I z\ p 

(the geometric arrangement of the spins corresponds to Fig. |6b). 



5.2 Abelian anyons 

Hamiltonian (157j) already lends itself to direct analysis. However, let us first reduce it to the 
more familiar form (j5J. We construct a new lattice A' so that the effective spins lie on its links 
(see Fig. [6b). This is a sublattice of index 2 in the original lattice A (here "lattice" means 
"translational group"). The basis vectors of A' are mi = ni — 112 and rri2 = ni + 112. The 
plaquettes of the effective spin lattice become plaquettes and vertices of the new lattice, so the 
Hamiltonian can be written as follows: 



H. 



cff 



\vertices plaquettes / 



where J cS = J 2 J 2 /(1Q\J Z \ 3 ). 

Now we apply the unitary transformation 



u = n x j n y * ^ 



horizontal links vertical links 

for suitably chosen spin rotations X and Y so that the Hamiltonian becomes 




H' eB = UH eS U^ = -Je ff YA S + V B p ), (39) 



where A s and B p are defined in Eq. (Caution: transformation (1551) breaks the translational 
symmetry of the original model.) 

The last Hamiltonian has been studied in detail [21]. Its key properties are that all the 
terms A s , B p commute, and that the ground state minimizes each term separately. Thus the 
ground state satisfies these conditions: 

= +|tt>, B p |tf) = +|tf>. (40) 

Excited states can be obtained by replacing the + sign to a — sign for a few vertices and plaque- 
ttes. Those vertices and plaquettes are the locations of anyons. We call them "electric charges" 
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and "magnetic vortices" , or e-particles and m-particles, respectively. When an e-particle moves 
around a m-particle, the overall state of the system is multiplied by —1. This property is stable 
with respect to small local perturbations of the Hamiltonian. (A local operator is a sum of 
terms each of which acts on a small number of neighboring spins.) It is also a robust property 
that the number of particles of each type is conserved modulo 2. 

The model has 4 superselection sectors: 1 (the vacuum), e, m, and e = e x m. The 
latter expression denotes a composite object consisting of an "electric charge" and a "magnetic 
vortex" . These are the fusion rules: 

exe — mxm — exe — 1, 

(41) 

e x m = e, e x e = m, m x e = e. 

(In general, fusion rules must be supplemented by associativity relations, or 6j-symbols, but 
they are trivial in our case.) 

Let us discuss the braiding rules. One special case has been mentioned: moving an e-particle 
around an m-particle yields — 1. This fact can be represented pictorially: 




The fist diagram shows the "top view" of the process. The diagrams in the second equation 
correspond to the "front view": the "up" direction is time. 

It is easy to show that e-particles are bosons with respect to themselves (though they clearly 
do not behave like bosons with respect to m-particles); m-particles are also bosons. However 
e-particles are fermions. To see this, consider two processes. In the first process two ee-pairs 
are created, two of the four e-objects are exchanged by a 180° counterclockwise rotation, then 
the pairs are annihilated. (Each e-object is represented by e and m, so there are 8 elementary 
particles involved.) In the second process the two pairs are annihilated immediately. It does 
not matter how exactly we create and annihilate the pairs, but we should do it the same way 
in both cases. For example, we may use this definition: 

\m e\ \m e\ A A 

creation = \ V / , annihilation = / / ■ \ . (43) 

V V |m e| \m e\ 

Now we compare the two processes. Each one effects the multiplication of the ground state by 
a number, but the two numbers differ by —1: 




Indeed, in the left diagram the dashed line is linked with the solid line. This corresponds to an 
e-particle going around an m particle, which yields the minus sign. 
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Figure 7: Weak breaking of the translational symmetry: e-vortices and m-vortices live on 
alternating rows of hexagons. 

Braiding an e-particle with an e- or m-particle also gives —1. This completes the description 
of braiding rules. 

It remains to interpret these properties in terms of vortices and fermions in the original 
model. Tracing sites and plaquettes of the reduced model (j3J) back to the original model, we 
conclude that e-particles and m-particles are vortices that live on alternating rows of hexagons, 
see Fig. Note that the two types of vortices have the same energy and other physical 
properties, yet they cannot be transformed one to another without creating or absorbing a 
fermion. A general view on this kind of phenomenon is given in Appendix [Fl 

The fermions in the original model belong to the superselection sector s, although they 
are not composed of e and m. In the perturbation-theoretic limit, the energy of a fermion is 
about 2\J Z \ whereas an em-pair has energy 4J e g <C \J Z \. The fermions are stable due to the 
conservation of W p (and also due to the conservation of energy). However, they will decay into 
e and m if we let the spins interact with a zero-temperature bath, i.e., another system that can 
absorb the energy released in the decay. 

6 Phase B acquires a gap in the presence of magnetic 
field 

6.1 The conic singularity and the time-reversal symmetry 

Phase B (cf. Fig. [5]) carries gapped vortices and gapless fermions. Note that vortices in this 
phase do not have well-defined statistics, i.e., the effect of transporting one vortex around the 
other depends on details of the process. Indeed, a pair of vortices separated by distance L is 
strongly coupled to fermionic modes near the singularity of the spectrum, |q — q*| ~ L^ 1 . This 
coupling results in effective interaction between the vortices that is proportional to e(q) ~ L^ 1 
and oscillates with characteristic wavevector 2q*. When one vortex moves around the other, 
the quantum state picks up a nonuniversal phase <p ~ L~H, where t is the duration of the 
process. Since the vortex velocity v = L/t must be small to ensure adiabaticity (or, at least, 
to prevent the emission of fermionic pairs), the nonuniversal phase (p is always large. 

The conic singularity in the spectrum is, in fact, a robust feature that is related to time- 
reversal symmetry. As discussed in the end of Section (2], this symmetry is not broken by fixing 
the gauge sector (i.e., the variables w p ) since the honeycomb lattice is bipartite. Let us show 
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that a small perturbation commuting with the time-reversal operator T can not open a spectral 
gap. 

We may perform the perturbation theory expansion relative to the gauge sector of the 
ground state (i.e., any term that changes the field configuration is taken into account in the 
second and higher orders). This procedure yields an effective Hamiltonian which acts in the 
fixed gauge sector and therefore can be represented in terms of Cj and Ujk- The operator T 
defined by (ITT]) changes the sign of Ujk = ib^ h b k ]k , but this change can be compensated by a 
gauge transformation. Thus we get a physically equivalent operator T' such that 

rpi * ( T i\-i _ * T , (T i\-x _ S Cj if j G even sublattice, 

Ujk{ > ~ Ujkl ' ~\ - Cj if j e odd sublattice. 

A T'-invariant perturbation to the fermionic Hamiltonian (I3TT) (corresponding to a fixed gauge) 
can not contain terms like iCjCk, where j and k belong to the same sublattice. Thus the 
perturbed matrix A(q) in Eq. ( 1321) still has zeros on the diagonal, though the exact form of 
the function /(q) may be different. However, a zero of a complex- valued function in two real 
variables is a topological feature, therefore it survives the perturbation. 



6.2 Derivation of an effective Hamiltonian 

What if the perturbation does not respect the time-reversal symmetry? We will now show that 
the simplest perturbation of this kind, 

V = ~ J2 ( h *°j + h v a 3 + ' ( 45 ) 

j 

does open a spectral gap. (Physically, the vector h = (h x , h y , h z ) is an external magnetic field 
acting on all spins.) For simplicity we will assume 

Let us use the perturbation theory to construct an effective Hamiltonian H e g acting on the 
vortex-free sector. One can easily see that = 0. Although the second-order term does 
not vanish, it preserves the time-reversal symmetry. Therefore we must consider the third-order 
term, which can be written as follows: 

H® = U VG' (E )VG' (E )VU , 

where Ho is the projector onto the vortex-free sector, and G' Q is the unperturbed Green function 
with the vortex-free sector excluded. In principle, the Green function can be computed for 
each gauge sector using the formula Gq{E) = —i e^ E ~ Ho+l5 ^ t dt (where 5 is an infinitely 
small number). For fixed values of the field variables Ujk the unperturbed Hamiltonian may be 
represented in the form (1181) and exponentiated implicitly by exponentiating the corresponding 
matrix A; the final result may be written as a normal-ordered expansion up to the second order. 
However, it is a rather difficult calculation, so we will use a qualitative argument instead. 

Let us assume that all intermediate states involved in the calculation have energy AE ~ \J\ 
above the ground state. (Actually, AE w 0.27| J\ for the lowest energy state with two adjacent 
vortices, see Appendix lAl) Then G' (E ) can be replaced by —(1 — n )/|J|. The effective 
Hamiltonian becomes 

H$ ~ -^£<^K. (46) 
i,k,i 
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where the summation takes place over spin triples arranged as follows: 




Configuration (a) corresponds to the term a^a^af = —iDiu^u^iCjC^ (where Di may be omitted 
as we work in the physical subspace), or simply —iCjCk in the standard gauge. Configuration 
(b) corresponds to a four-fermion term and therefore does not directly influence the spectrum. 
Thus we arrive at this effective Hamiltonian: 

j,k 

A = 2J(^)+2k( ), ( 48 ) 

h x hy h z 



Here ( — — ) is just another notation for u std , i.e., the matrix whose entry ( — — )jf t is equal to 
1 if there is a solid arrow from k to j in the figure, —1 if an arrow goes from j to k, and 
otherwise. ( — — ) is defined similarly. 

6.3 The spectrum and the Chern number 

The fermionic spectrum e(q) of the Hamiltonian (14*81) is given by the eigenvalues of a modified 
matrix iA(q) (cf. Eq. (152]) ): 

^(q) = (_f / ( ( q q ) r .1^) , e(q) = ±V / l/(q)l 2 + A(q) 2 , (49) 

where /(q) = 2j(e^ q ' ni ^ + e^ q,n2 ^ + l) and A(q) = An (sin(q, ni) +sin(q, — n 2 ) + sin(q, n 2 — ni)J . 
Actually, the exact form of the function A(q) does not matter; the important parameter is 

A = A(q,) = -A(-q,) = 6V3K ~ (50) 
which determines the energy gap. The conic singularities are resolved as follows: 

e(q) w ±v/3J 2 |5q| 2 + A 2 , 

(51) 

where Sq = q — q* or 5q = q + q*. 

Remark 6.1. The magnetic field also gives nontrivial dispersion to vortices. Indeed, the 
operators W p are no longer conserved, therefore a vortex can hop to an adjacent hexagon. 
Thus the vortex energy depends on the momentum. This effect is linear in h, but it is not so 
important as the change in the fermionic spectrum. 
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Let us also find the fermionic spectral projector, which determines the ground state. The 
global spectral projector P is defined by Eq. f l27|) ; we now consider its Fourier component: 



^(q) = \ t 1 " sgn(zA(q)) = - (l + m x {q)a x + m y {q)a y + m z {q)a z ) , (52) 



(-Sq y , -Sq x , -^j) if q ~ q*, 



v/(^) 2 +(<5g a ) 2 +AV(3J 2 ) 

m(q) « <( (53) 

TCTri^vc^ H*» Sq *> A) if <* « 

The function m maps the torus to the unit sphere. If A > 0, then this map has degree 1. 
Indeed, the neighborhood of q* is mapped onto the lower hemisphere, the neighborhood of — q* 
is mapped onto the upper hemisphere; in both cases the orientation is preserved. (The rest of 
the torus is mapped onto the equator.) For negative A the map has degree —1. 

An important topological quantity characterizing a two-dimensional system of noninteract- 
ing (or weakly interacting) fermions with an energy gap is the spectral Chern number. It plays 
a central role in the theory of the integer quantum Hall effect [4~3| HU [1] . In our model there 
is no analogue of Hall conductivity (because the number of fermions is not conserved), but the 
Chern number determines the edge mode chirality and anyonic properties of vortices (cf. |48j). 

The spectral Chern number is defined as follows. For each value of the momentum q we 
consider the space L(q) of annihilation operators, i.e., fermionic modes with negative energy; 
this is the subspace the matrix P(q) projects onto. Thus we obtain a complex vector bundle 
over the momentum space. (In our case L(q) is a one-dimensional subspace of C 2 , so the 
bundle is one-dimensional.) The first Chern number of this bundle is denoted by v and can be 
expressed as follows (cf. |44J): 

" - lih^ p ^ - i/^ES-SS)^- (54) 

The Chern number is always an integer. If the spectral projector -P(q) is given by Eqs. (152|) . 
(J53D, then 

1 f ( dm dm \ , , . , 1 , KKN 

V = Tn J [dg X Wy "J ^ % = Sg = ^ 

We will use the notation B v (where v = ±1) to designate phase B in the magnetic field. In the 
Abelian phases A x , A y , A z the Chern number is zero. 



7 Edge modes and thermal transport 

Remarkably, any system with nonzero Chern number possesses gapless edge modes. Such modes 
were first discovered in the integer quantum Hall effect [51]; they are chiral, i.e., propagate only 
in one direction (see Fig. [8]). In fact, left-moving and right-moving modes may coexist, but the 
following relation holds [52J: 

Vedge = (# of left-movers — # of right-movers) = v. (56) 
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a) 



b) 



Figure 8: Chiral edge modes: left-moving (a) and right-moving (b). 



In the absence of special symmetry, counterpropagating modes usually cancel each other, so 
the surviving modes have the same chirality. A calculation of the edge spectrum in phases B v 
(for some specific boundary conditions) and a simple proof of Eq. (j56p are given in Appendix IB1 
More rigorous and general results, which even apply to disordered systems, can be found in 



It is important to note that the analogy to the quantum Hall effect is not exact. In our 
model (like in two-dimensional superfluid and superconducting systems [551 HH]) edge modes 
are described by real fermions, as contrasted to complex fermions in the quantum Hall effect. 
Therefore each quantum Hall mode is equivalent to two modes in our system. 

Chiral edge modes can carry energy, leading to potentially measurable thermal transport. 
(The temperature T is assumed to be much smaller than the energy gap in the bulk, so that 
the effect of bulk excitations is negligible.) For quantum Hall systems, this phenomenon was 
discussed in Refs. [511 EZ]- The energy current along the edge in the left (counterclockwise) 
direction is given by the following formula: 



where c_ is some real number. (The factor tt/12 is introduced to make a connection to conformal 
field theory, see below.) It is remarkable that c_ does not depend on particular conditions at the 
edge, but rather on the bulk state. Indeed, the energy current is conserved, therefore it remains 
constant even if some conditions change along the edge. The effect is invariant with respect 
to time rescaling. Since the energy current has dimension (time) -2 , it must be proportional to 
T 2 . But the value of the dimensionless proportionality coefficient cannot be found using such 
simple arguments. 

There are two standard ways to calculate the coefficient c_. They both rely on certain 
assumptions but can be applied to our model, yielding this result: 



The first argument [HE] (adapted to real fermions) assumes translational invariance and the 

absence of interaction. Each edge mode is described by a free fermion with an energy spectrum 
s(q) such that s(—q) = —e(q) and e(q) — > ±oo as q — > ±oo. The signs in the last two expressions 
agree if the mode propagates in the direction of positive q (for simplicity we may assume that 
e(q) > when q is positive and e(q) < when q is negative). Thus the Hamiltonian has the 
form 



Refs. [531 EH- 



v 







g 



q:e(q)>0 
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If e(q) > 0, then a q is an annihilation operator and a_ 9 = a J is the corresponding creation 
operator. The mode propagates with group velocity v (q) = de/dq, and the occupation number 
n(q) is given by the Fermi distribution. The energy flow due to each mode propagating in the 
positive direction can be calculated as follows: 

dq f e(q) de dq 1 f°° ede n 2 



Each mode propagating in the opposite direction contributes —1%, therefore I = jg^T 2 . 

The second derivation [57] is based on the assumption that the edge modes can be described 
by a conformal field theory (CFT). In this case, 

c_ = c — c, (59) 

where c and c are the Virasoro central charges. Thus, c_ is called the chiral central charge. Left- 
moving fermions have (c, c) = (|,0) whereas right-moving fermions have (c, c) = (0, |), which 
implies Eq. ( |58|) . More generally, c and c are some rational numbers, and so is c_. The chiral 
central charge parametrizes a two-dimensional gravitational anomaly [59] of the corresponding 
CFT; it can also be identified with the coefficient in front of the gravitational Chern-Simons 
action in a three-dimensional theory [5j. Volovik [50] suggested that for 3 He-A films, the role 
of gravitational field is played by the order parameter interacting with fermions. However, it 
is not obvious how to define a "gravitational field" for lattice models. 

It remains a bit mysterious how the chiral central charge is related to the ground state 
and spin correlators in the bulk. This question is partially answered in Appendix [Dl but the 
obtained expression for c_ is not easy to use, nor can we demonstrate that c_ is rational. Note 
that there is a beautiful relation between the chiral central charge and algebraic properties of 
anyons [TJEI], which does imply the rationality of c_. We discuss that relation in Appendix [El 
(see Eq. (11721) on page [66]), though it is unclear how to deduce it only considering the bulk. The 
only known argument is to assume that edge modes are described by a CFT, then one can use 
modular invariance [22]. In fact, the modular invariance alone would suffice. In Appendix [Dl 
we try to derive it from general principles, but again, encounter a problem. 



8 Non-Abelian anyons 



We continue the study of phase B in the magnetic field. Now that all bulk excitations are 
gapped, their braiding rules must be well-defined. Of course, this is only true if the particles 
are separated by distances that are much larger than the correlation length associated with the 
spectrum (1511) . The correlation length may be defined as follows: £ = | Iraqi" 1 , where q is a 
complex solution to the equation e(q) = 0. Thus 







J 3 


A 







(60) 



The braiding rules for vortices depend on the spectral Chern number v. Although v is 
actually equal to +1 or —1 (depending on the direction of the magnetic field), one may formally 
consider a model with an arbitrary gapped fermionic spectrum, in which case v may take any 
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integer value. We will see that vortices behave as non-Abelian anyons for any odd value of v, 
but their exact statistics depends on v mod 16. 

The properties of the anyons are summarized in Table [j] on page [301 The notation and 
underlying concepts are explained below; see also Appendix [El Let us first show a quick 
way of deriving those properties from conformal filed theory (CFT) in the most important 
case, v = ±1. (For a general reference on CFT, see [S3]-) Then we will give an alternative 
derivation, which uses only rudimentary CFT and refers to the operational meaning of braiding 
and fusion. 

8.1 Bulk-edge correspondence 

The properties listed in Tabled] form the same type of algebraic structure that was described by 
Moore and Seiberg [I] in the CFT context. However, the actual connection to CFT is indirect: 
anyons are related to edge modes, which in turn can be described by a field theory in 1 + 1 
dimensions. More concretely, the space-time may be represented as a cylinder (see Fig. [9]). It 
is convenient to use the imaginary time formalism (t = —ir, where r G 1R), so that we have a 
two-dimensional Euclidean field theory on the side surface of the cylinder. The surface may be 
parametrized by a complex variable z = r + ix, where x is the spatial coordinate. 

The two-dimensional field theory describes the physics of the edge, which is generally richer 
than that of the bulk. The theory possesses both local and nonlocal fields. The insertion of 
a nonlocal field <f>(r + ix) corresponds to an anyonic particle emerging on the edge or sinking 
into the bulk at point (r,x). The correlation function of several nonlocal fields has nontrivial 
monodromy which coincides with the anyonic braiding. Specifically, the counterclockwise ex- 
change of anyons in the bulk is equivalent to moving the fields counterclockwise on surface (if 
we look at the cylinder from outside). Moreover, the value of the correlator is not a number, but 
rather an operator transforming the initial anyonic state (on the bottom of the cylinder) into 
the final one (on the top of the cylinder). For non-Abelian anyons, the space of such operators 
is multidimensional. 

The anyon-CFT correspondence has been successfully used in the study of quantum Hall 
systems [T^l [T5l |3T] . The correspondence is well-understood if all boundary fields are either 
holomorphic or antiholomorphic, which is the case for our model. We have seen that the edge 
carries a left-moving (holomorphic) fermion for v — 1, or a right-moving (antiholomorphic) 
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Superselection sectors: 1 (vacuum), e (fermion), a (vortex). 

Quantum dimension: d± — 1, d £ — 1 d a — 

Topological spin: 9\ — 1, 9 £ — —l, 9 a = 9 = exp(|u/) ; 

Frobenius-Schur indicator: Xi = 1, x e = 1, x CT = x = (— l)^ 2-1 )/ 8 . 

Global dimension: £> 2 = f J2 U = 4- 

e x a = a, 



Fusion rules: e x e — 1, 

Associativity relations: 



ffXff = l + e. 



a o a oac a o a 

V / _X \ V _X_ 

o a a 



a a a 
X \V 



X 



V2 1 V2 



E O E £ O £ 



o £ a a £ o 





£ £ £ £ £ £ 

"V = \r 



e a a £ a a 



a o £ 



a £ £ O £ £ 

\ v 



£ £ a £ £ a 

V / 



a j 



a £ a a £ a 



1 1 

Braiding rules: 



E a a £ a a 



a a £ o o £ 



Definition of R x z y : 



y * 



y x 



it>f = -1, 



R £ f = R 





a 1 i 



R1 a = xexp(-fw), 
B?° = xexp(^z/). 



Topological ^-matrix: 



def 1 
xy J) 




Si 
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1 


1 






2 


2 






\ 
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V2 
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(S £ ) =exp(-f«/), 



Table 1: Algebraic properties of anyons in non-Abelian phases {y is odd). 
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fermion for v = — 1. A vortex emerging on the surface corresponds to a twist field a. The 
correlation functions for such fields are given by holomorphic (resp. antiholomorphic) conformal 
blocks for the Ising model. 

A partial bulk-edge correspondence can be established at a more elementary level. Let 
v — 1. A particularly important parameter of the edge theory is the conformal weight of the 
twist field, {h a ,h a ) = (^,0). The related bulk parameter is the topological spin of a vortex, 
Oct] it corresponds to the counterclockwise rotation by angle 2tt. Since the rotations in the bulk 
and on the surface must agree, the following equation holds for an arbitrary field a: 

9 a = e^a-fcO. (61) 

Thus 9 a = e l7r / 8 . Similarly, if v = —1, then 9 a = e~ i7r//8 . In the general case, there are 
|z/| = 2n + 1 species of the free fermion; the twist field also comes with some multiplicity!! 
One may argue that the conformal weight of the twist field is proportional to the number of 
fermionic species. Hence 

Q a = e lnu/8 . (62) 



8.2 Unpaired Majorana modes 

We now begin a rather lengthy derivation of the properties listed in Table [TJ dealing only with 
anyons in the bulk. First, we give a crude description of vortices by Majorana operators and 
find the fusion rules. The braiding rules and so-called associativity relations (also known as 
crossing symmetry, or 6j-symbols) are determined up to several free parameters; we discuss 
what combinations of these parameters have invariant meaning. Solving the so-called pen- 
tagon and hexagon equations, we reduce the continuum of possibilities down to eight consistent 
theories. The right theory is selected using Eq. (16"2"1) . 

The starting point is this: if v is odd, then each vortex carries an unpaired Majorana mode. 
A sufficiently rigorous proof of this statement can be found in Appendix O Here we give a 
rough explanation based on the quantum Hall analogy. 

It is known that the quantized Hall conductivity for noninteracting electrons (in units of 
e 2 /h) is equal to the Chern number of the projector onto the occupied electron states [44J. 
The essential difference from our case is that electrons are ordinary fermions, not Majorana 
fermions. However, we can concoct an analogue of an electron system from two copies of the 
Majorana system, which may be pictured as two layers. The Hamiltonian is as follows: 

i x i 

-^electron = 4 A 3 k ^ C 'j C 'k + C 'j C k) = ^2 iA jk a ] a k, where Oj = -(c'j + ic") . (63) 

j,k j,k 

(We forget about the original spin model for the purpose of this construction.) Note that 
the Hamiltonian possesses a global 0(2) symmetry that consists of orthogonal linear transfor- 
mations c'j i— > ac'j + j3c", c" i— > •yc'j + 5c". The rotational subgroup U(l) = S0(2) C 0(2) 
corresponds to the conservation of electric charge, whereas reflections (i.e., transformations 
with determinant —1) change the sign of the charge. 

3 Specifically, there are 2" copies of the twist field, which are transformed according to the fundamental 
representation of Spin(2n + 1). 
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Assuming that the Fermi energy is zero, the projector onto the occupied electron states 
coincides with the spectral projector for the Majorana system. Thus the Hall conductivity of 
the electron system is v. 

A vortex piercing both Majorana layers corresponds to half-integer magnetic flux. Such a 
vortex carries excessive charge q = u/2 + n, where n is an arbitrary integer. In particular, if 
v is odd, a state with q = 1/2 exist. By an 0(2) reflection, it is related to a q = —1/2 state 
with the same energy. The doublet of states with unit charge difference may be attributed to 
a zero-energy electron mode that can be empty or occupied. This mode can be represented by 
two Majorana modes. Since the layers are independent, we conclude that each layer contains 
one zero-energy Majorana mode. 

It has been previously shown that zero-energy Majorana modes exist in some exotic one- 
dimensional systems [64] as well as vortices in two-dimensional p-wave superconductors [48, 65J. 
A pair of Majorana modes at two vortices constitutes a full fermionic mode with a two- 
dimensional Fock space. The quantum state of such a pair is virtually inaccessible to mea- 
surements or perturbations as long as the vortices stay far apart from each other. A system 
of 2n vortices possesses a protected space of dimensionality 2 n (or 2 n , if we require that the 
whole system have even fermionic parity). 

8.3 Fusion and braiding rules 

In our model a vortex carrying an unpaired Majorana mode is just one of the superselection 
sectors, denoted by a. The other sectors are 1 (the vacuum) and e (a fermion). If two vortices 
fuse, they either annihilate completely or leave a fermion behind: a x a = 1 + e. The actual 
fusion outcome depends on the initial quantum state. Hence the protected space of the vortex 
pair has two basis vectors: \ipi a ) and \ip° a ). (The upper indices indicate the particle types 
before the fusion whereas the subscript indicates the resulting particle.) The complete set of 
fusion rules is as follows: 

e x e = 1, e x a = a, o x a = 1 + e, (64) 

plus trivial rules of the form 1 x x — x. Read backwards, these relations are understood as 
splitting rules: for example, an e-particle can split into two cr-particles. However, these rules 
do not capture more subtle aspects of fusion and splitting, which will be discussed later. 

Braiding rules for Majorana half-vortices in a spin-triplet superconductor have been derived 
by D. Ivanov [63]. Here we follow the main idea of Ivanov's work. We should, however, keep in 
mind one important difference between his setting and our model. A spin-triplet superconductor 
has a locally measurable vector order parameter, which contributes to vortex- vortex interaction 
and can interact with impurities. One vortex making a full turn around another may pick up 
a nonuniversal phase, hence the non-Abelian statistics is defined up to arbitrary phase factors. 
That is not the case for our model (or for spinless superconductors [IE]), so additional arguments 
are required to find the Abelian part of the vortex statistics. 

Once again, we use the fact that each vortex p carries an unpaired Majorana mode C p , 
which is a linear combination of the operators Cj on neighboring sites. The operators Cj do 
not commute with gauge transformations and therefore should be used with care. The gauge 
can be fixed in a neighborhood of each vortex, so constructing the linear combination is not 
a problem. However, the overall sign of C p does not have invariant meaning. This ambiguity 
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is avoided if we consider fermionic path operators (1161) . which are gauge-invariant. A suitable 
linear combination of elementary paths constitutes a path that begins or ends at a vortex. 

Let us choose some reference path l p connecting each vortex p = 1, 2 . . . to a reference point 
0. We will assume that the vortices lie on the horizontal axis whereas the reference point is 
located in the lower half-plane. Fixing the gauge along the paths, we may write 

W{l p ) = C p c . (65) 

Let us exchange vortices 1 and 2 by moving them counterclockwise. The exchange process 
is described by a unitary operator R acting on the physical Hilbert space. It also acts on 
operators by conjugation: X i— > RXR^. Clearly, 





RW(h)tf = W(l[) = W(l 2 ), 
RW(l 2 )R^ = W(l' 2 ) = -W(h) 



(66) 



(We have used the fact that W{1' 2 ) = —W(li). Indeed, the paths l 2 and h differ by a loop 
enclosing a vortex; transporting a fermion around the vortex gives rise to the minus sign.) 
Hence 

RdR^ = C 2 , / 7T \ 

where 9 is a phase factor. (We will see that 9 is actually the topological spin of a vortex; for 
now it is just an unknown parameter). 

At first sight, it is not clear whether the number 9 has an invariant meaning. Indeed, the 
operator that moves a vortex along a given path may be defined up to an arbitrary phase. 
However, the ambiguity can be avoided by a careful definition of the particle exchange process. 
The following argument is completely general; it is not based on adiabaticity or translational 
invariance. 

Let us consider a vortex path operator W a (l) that is composed of elementary steps, i.e., 
displacements of the vortex from hexagon to hexagon]^ Each elementary displacement is defined 
up to a phase, but the path operator satisfies these equations: 

W a (l) = W a (l)\ W a (hh) = W a {l 2 )W a (h). (68) 

Here I denotes the reverse path, whereas l\l 2 is the composite path (the vortex goes along li 
and then along l 2 ). Remarkably, two vortices can be exchanged in such a way that the arbitrary 
phases cancel each other, see Fig. [TUJ The vortices move along the indicated lines; each line 
segment is passed by a vortex in both directions. 

As mentioned above, the vortex pair has two states corresponding to the possible fusion 
outcomes, \ipi a ) and \ip° a ). They should be identified with the eigenvectors of CiC 2 , but we do 
not know which is which. We can only write: 

CiC 2 \W) = HVO, CiC 2 |Vr> = -HC 7 ), (69) 



4 The displacement may be realized as the action of a a on one of the two spins at the boundary between the 
hexagons (which changes the field configuration) followed by an operator of the form exp(^ FkmCkC m ) (which 
adjusts the fermionic subsystem). 
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Figure 10: A realization of the particle exchange in which all local contributions cancel, revealing 
the topological effect in a pure form. 



where a = ±1 is unknown. Thus the braiding operator R acts as follows: 

R\ij}™) = R™\il)% a ) : R™ = 6e~ iaw/4 , 

where , (70) 

R\ip™) = R™^™), R7 = 0e ian/A . 

Remark 8.1. One may wonder why formula (1671) contains the minus sign in front of Ci, but 
not in front of C 2 . This is, in fact, a consequence of the convention that the reference point 
lies in the lower half-plane. If we move it to the upper half-plane, the signs will change, a will 
turn into —a, and the signs in front of a in Eq. ( 1701 will change too. The numbers R° a , R 
and 9 are invariant though. 



8.4 Associativity relations 

The fusion rules only indicate what fusion and splitting events are possible. Nontrivial relations 
arise if we consider a sequence of such events. 

8.4.1 Relations based on the anticommutativity of fermionic operators. Let us 

consider the splitting of a cx-particle into e, a, and e. This can be done in two different ways, 
depending on the order the two e-particles are created. Let us suppose that the cx-particle is 
located at point 2 in the middle; the e-particles will appear at points 1 and 3 on the left and 
on the right, respectively. The particles (or their future locations) are connected by the paths 
li2 and l 23 . 

The er-particle is described by the associated Majorana mode C 2 . The e-particles can be 
created from the vacuum by operators a\ and a\. However, these operators are not physical 
by themselves. To actually split the first e-particle off the cx-particle, one needs to apply a 
fermionic path operator Vl^(Zi 2 ), which equals a\C 2 in a suitable gauge. The second e-particle 
is created by the operator W(l 2 3) = a\C 2 . Since W{li 2 )W{l 2 ?) = — W (l 2 z)W (h 2 ) , we obtain 
the following associativity relation between the two splitting processes: 



£ a £ £ a £ 




(71) 



a a 



Similarly, we can consider the splitting of one e-particle into three. For this, we use two 
conventionstj (i) moving an e-particle from place j to place k is described by the operator a\aj] 

5 Although these conventions are somewhat arbitrary, the result is invariant (see subsection 18. 4. 2|) . 
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(ii) creating a fermionic pair from the vacuum is described by 0!%, where j is located left of 
k. Thus we get: 



£ £ £ £ £ £ 




£ 



because (a\a\) (a^) = (a^Xajc^). ( m these graphs the "up" direction is time. The interme- 
diate state of the middle particle is shown in all cases, but lines corresponding to label 1, i.e., 
the vacuum superselection sector, are suppressed.) 

8.4.2 Some generalities. The analysis of other splitting processes (e.g., a — > aaa) is more 
complicated. Ideally, we would use the vortex path operator, but its exact form is unknown. 
Therefore we have to rely on more abstract arguments. Before we proceed, let us clarify some 
points pertaining to anyons in general. 

1. Speaking about anyons, we are interested in particle types (i.e., superselection sectors) 
and topology of the particle worldlines (braids, trees, etc.). An actual particle is also 
characterized by position and local degrees of freedom; for example, an e-particle may be 
"elementary" or consist of two adjacent vortices. We generally ignore such details though. 

2. In the study of splitting and fusion, we consider particles located on the horizontal axis. 
In this case, we only care about the order of the particles on the line, but not about their 
positions. All configurations with the same particle types and order can be identified with 
each other. Unlike in the two-dimensional case, cycling through several configurations has 
trivial effect on the quantum stated 

3. "Quantum state of a particle" is a rather subtle notion. It may be understood as a 
projector that enforces certain spin correlations in some neighborhood of the supposed 
particle location. (The neighborhood radius must be much larger than the correlation 
length). Such an object has no "overall phase". Superpositions of states from different 
superselection sectors cannot be constructed either. 

4. In spite of that, the state {i/jf 7 ) has a well-defined phase if we consider it relative to the 
vacuum. A similar argument holds for an arbitrary state \ip^ v ), which describes particles x 
and y with total "anyonic charge" z. For a general anyonic system, the x?/-pair may have 
several distinct states that belongs to the same superselection sector z; in other words, 
there may be several ways to split z into x and y. Such states form a finite-dimensional 
Hilbert space V^ y , which is called fusion spaced In our model, all such spaces have 
dimension one or zero. For example, |V ; i <J ) is a unit vector in the one-dimensional space 
V™ . 

6 This claim can be justified as follows. To move a particle z from one place to another, we apply a path 
operator W z (r, r'), which acts on the spins in some neighborhood of the interval [r, r']. These operators are 
defined in such a way that it does not matter whether the particle is moved at once or in several steps (cf. 
Eq. (|68|) ). When several particles are present, we do not let them pass through each other, which imposes a 
certain restriction on the sequence of operators applied. All such sequences are actually equivalent since the 
path operators for nonoverlapping intervals commute. 

7 Technically, it should be called "splitting space", whereas the fusion space is its dual, V* — (V^ y ) . 
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For a more rigorous definition of the space V* y , consider the splitting of z into x and y by an 
operator L that acts on the spins in some fixed finite region Q. (We call such operators local.) 
Let us also fix a quantum state \^> z ) that has a z-particle at a given place and no other particles 
in Q or within the correlation length from Q. Finally, we consider the set of local operators 
L for which the state L\ty z ) has an x-particle and an y-particle at the required places. By 
definition, these are all possible states of the x?/-pair that can be obtained from z. On the 
other hand, such states are in one-to-one correspondence with equivalence classes of operators 
L: we say that L and V are equivalent if L\§* z ) = L'\^ z ). For local operators, the equivalence 
relation does not depend on the choice of |\& z ). Thus we arrive at the following definition: 

V z y is the set of equivalence classes of local operators that split z into x and y. 

Each vertex in a splitting graph (as in Eqs. flTTj) . (1721) above) designates an equivalence class 

of local operators. For example, ^ = \ip £ a a ) € V £a ' . One may choose an arbitrary basis in 

the space V z y . In our case, this amounts to fixing the phase of the vectors IV^), \ r 4 ) ° a ), \^ £ f)i 
|^ £ ), and \ip\ e ). Different bases are related by a transformation of the form 

W? Y = «7W>, K v \ = i, (73) 

which generally affects relations like (17T1) and (1721) . However, these particular relations are 
invariant since they have the same set of basis vectors on both sides: IV^), \ip^ e ) in Eq. (1711) 
and \i){ £ ) in Eq. (1721. 

8.4.3 More relations. Let us write some other possible associativity relations in a general 
form: 

eoa eoa acre oae 



a e £ £ £ o 



a £ a a £ a 





a £ a o £ o 



* V" > A/ = r \X > (76) 



£ 



a 

1 1 



EOO £aa a a £ a a £ 



11 11 

The numbers ft, ft, 7i, 72, 5, A*i, A*2 are equal to 1 in magnitude. These and similar 
coefficients are called structural parameters. 

Using the transformation (1731) . we can eliminate some of the parameters, e.g., ft, ft, and 
7i. Indeed, they are transformed as follows: 
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We require that f3[ = f3' 2 = ^[ = 1 and solve for u x z v . Note that the solution is unique up to 
trivial variations not affecting the structural parameters, namely 

< y - CxCvC 1 ^. (79) 

Finally, let us consider the splitting of one cr-particle into three cr-particles. The space V° au 
is two-dimensional. A basis in this space can be chosen in two ways: 




Of course, one basis can be expressed in terms of the other: 

16) = aii|»7i) + aBi\Ve), \^ e ) = a lE \r]i) + a e£ \r] £ ) , where a xy = (r] x \Q. (81) 

The coefficients a xy and the other structural parameters will be found later. At this point, we 
can only tell that the matrix (a xy ) is unitary. 

To illustrate the physical meaning of the structural parameters, let us return to the gedanken 
experiment considered in the introduction. Suppose we create two pairs of cr-particles from the 
vacuum. Then we take one particle from each pair and fuse them. With probability |an| 2 they 
will annihilate, and with probability |a £ i| 2 they will fuse into an e-particle. (We will see that 
|an| 2 = |a el | 2 = l/2.) 



8.5 Algebraic consistency 

All the previous arguments were based on the result that vortices carry unpaired Majorana 
modes. Not surprisingly, this property alone is not sufficient to fully characterize the fusion 
and braiding of anyons. We now invoke some additional principles. The first one is consistency: 
successive fusion and braiding events must commute with each other in certain cases. A more 
careful statement of this requirement amounts to the formulation of an algebraic theory of 
anyons. 

Anyons may be described in the framework of topological quantum field theory (TQFT), 
which originates from Witten's paper on quantum Chern-Simons fields j5] and the work of 
Moore and Seiberg on conformal field theory [3]. Important mathematical studies in this area 
were done by Reshetikhin and Turaev [66] and Walker [67]. For our purposes, it suffices to use 
a construction called unitary modular category (UMC), which constitutes the algebraic core 
of TQFT [6B] ■ This construction will be outlined in Appendix [El Actually, the full theory 
is not necessary to understand the calculations below. On the contrary, one may use these 
calculations to motivate some of the UMC axioms. 

One of the axioms is known as the pentagon equation, see Fig. ITTa on page [69j It deals with 
the five ways to split a particle u into x, y, z, w, or five representation of the space V^ yzw . The 
arrows in the figure may be regarded as equality signs, therefore the diagram must commute. 
For example, consider the splitting process a — > eaaa via intermediate states q = 1 and p = a. 
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The upper path across the diagram looks like this: 



e o a a 



an 



eoaa eooct £ a a a eooa 

V 

+ a, 



V 





an \ ) - a £ i 

'a Vo 



whereas the lower path is as follows: 



e a a a 



£ a o a 



e a o a 





/^l l 




c a a a e a cj a 

/'I "2"- I; \ ^ + /il72«££ 





Thus /ii72ai e = an, /ii72« ee = — a e i- One can also start with the splitting graph that has 
intermediate state e instead of 1; this gives another pair of relations between the structural 
parameters. 

Overall, there are 23 spaces V* yzw with nontrivial upper indices. Thus we obtain 23 pen- 
tagon equations, 5 of which are satisfied automatically. The remaining 18 equations imply the 
following relations between the structural parameters: 



VI 



VI 



^l72«l£ 


= an, 




— ~ a eli 


Piau 


— a lei 




&££) 




— Oil, 


r^iOu 


= 7l«£l, 


8(3 x a eE 


= Ole, 




= 7i a £ £ ; 



v r : 



yU27i fl £i = on, 

Glides = -0-le, 
fadle — ~ a ee'i 

T/j, 2 a n = 7 2 an, 
S/j, 2 a El = 7 2 a ee , 
rfoau = an, 
5(3 2 a ££ = a B x\ 



(82) 



33) 



V" 



a XX ra le a eX 

o £ iOii + ™ £e a e i 
aiiOie + rai e a ee 



a e iai £ + ra 



V{ 



1 v e 1 



V{ 



-\raeea 
' e 

uecre prefects 



££ 
2 _ 



0. 

1; 



r = 7i72, 
S 2 = 7i72, 
5 = -r, 



T TO (J (J (J 
* £ 



a\x + = 0, 

a £l axx + 5a ee a £ i = f3 2} 
axxaxe + Sa l£ a ££ = fix, 

a £l a l£ + 5a 2 ££ = 0; 



V{ 
V{ 
VI 



££(7a \T££(J(J 
' £ 



00££ TJ(J(J££ 



£(J(T£ \r£<J(J£ 
1 V £ 



flx(3xl2 = 1, 

A*i^2 = A/3 2 - 



Equations (lezHHol) have only two solutions satisfying the convention that (3 X = P 2 
(3x = (32 = 71 = 72 = Hi = fJ>2 = r = 1, S = -1, 



an ai e 



X /l 1 

71 V 1 - 1 ' 



^4) 



7i 



where x = ±1. 



35) 



1: 



(86) 
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Remark 8.2. The number of independent equations can be reduced by using a symmetry 
between upper and lower indices: V xy = Vg x , where z denotes the antiparticle for z (in our 
case, e = e, a — a). However this symmetry generally involves nontrivial phase factors, so we 
find the brute-force calculation a safer approach. 

Let us now examine consistency between fusion and braiding. Braiding is fully characterized 
by the action of the counterclockwise rotation on fusion spaces: R xy : V xy — > V yx . Since in 
our case the spaces V xy and V yx are one-dimensional, the linear map R xy is given by a single 
matrix elemental 

%»lf> = i?m where R xy = ^\R xy W z y ) ■ (87) 

y x 

In graphic notation, R xy \ip xy ) = (J) , therefore Eq. ( 1H71) reads: 



y x y x 

^ = EC? i^J . (88) 

z z 

Nontrivial relations arise if we consider the action of braiding operators on the fusion space 
of three particles. These relations can be expressed by two commutative diagrams called the 
hexagon equations, see Fig. [20] on page [HTJ The arrows in those diagrams may be understood 
as equalities of vectors in the space V yzx . Let us consider the following example of the first 
equation (the first line below corresponds to the upper path across the hexagon, the second to 
the lower path): 



aoa a a a a a a aoa a a a oaa 




a a a a a a a 



Thus we get: 1 + R a a e = x^2(R™) 2 , 1 - R™ = x^R™ R™ . 

The full set of equations arising from the first hexagon is as follows: 

r i + i?f = W2W) 2 , 

v™ a . I i - R a a £ = xV2RrR7, K eee : (R? f = i; (89) 

{ i + i?^ = -xv / 2(i?r) 2 ) 



8 In Appendix [Ej we take the liberty to omit the hat from the notation R x z y . This should not cause confusion 
because we do not consider matrix elements there. 
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Vf™, Vf £a : RfRf = Rf, Vf a£ : (Rf) 2 = Rf , 

yeaa^yaea. R ae R aa = _ R acr^ yaae. _(^2 = ^ (gg) 

V; ff££ , V £a£ : Rf = -1, 1/ ££ff : -(/^ e ) 2 = 1- 

The second hexagon is obtained from the first by replacing R x z v with (i? -1 )^ = (Rf)' 1 . For 
example, the equation RfRf = Rf becomes (Rf)~ 1 (Rf)- 1 = (Rf)' 1 - 
Eliminating redundancies, we get this system of equations: 

3,(1 -L- R a£ ) 

Rf = -1, Rf = Rf = ±i, (Rf 7 ) 2 = 1 a Rf 7 = -RfRf, (91) 

where x = ±1. The system has eight solutions, which all fit Eq. (1701) . Specifically, the solutions 
have the form 

Rf = -1, Rf = Rf = -ia, Rf = 9e~ iaw/ \ Rf = 9e iaw/ \ (92) 
where the following combinations of 9, a, and x are possible: 



9 


e -7«7r/8 


„—5in/8 


g— 3i7r/8 


g-wr/ 8 


giTr/8 


g3i7r/8 


g5«7r/8 


e 7in/8 


a 


+ 1 


-1 


+ 1 


-1 


+ 1 


-1 


+ 1 


-1 


X 


+ 1 


-1 


-1 


+ 1 


+ 1 


-1 


-1 


+ 1 



(93) 



Thus the properties of anyons are defined by the value of 9, which satisfies 9 8 = — 1. 



8.6 Final details 

To choose the the right solution and to complete the description of non-Abelian anyons, we use 
Eq. (162]) in conjunction with some general results from Appendix [El 

The matrix element an = (?7i|£i) (cf. Eqs. flHUl) . (j5T|) . flSBl) ) is equal to x/v2- But according 
to Eq. (1187!) 011 P a g e EH °n — ^a/da, where d a is the quantum dimension of the cr-particle, 
and x a is the Frobenius-Schur indicator. Therefore 

d a = V2, x CT = x. (94) 

The topological spin of the cr-particle can be computed using Eq. (121 ip on page [82j 

9 a = df(Rf + Rf) =9. (95) 

This result also follows from Eq. (I215p . Matching it with the expression of 9 a in terms of the 
Chern number, we represent all eight cases in (1931) by these formulas: 

9 = e l ^ s , a = (-l)^" 1 )/ 2 , x = (-l)^ 2 - 1 )/ 8 . (96) 

Thus we have obtained almost all properties in Table [H The remaining structure, namely 
the topological S-matrix, can be found using general rules (see Appendix [E]). 
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Superselection sectors: 1 (vacuum), e (fermion), e, m (vortices). 



Quantum dimension: 


d\ = 1, 


d e = 1 


d e = 




Topological spin: 


0i = l, 


Oe = "I, 


#e = 


6 m = 6 = exp(f«/); 


Frobenius-Schur indicator: 


Xi = 1, 


X e = 1, 


X e = 


= x m = x = exp(fii/). 


Fusion rules: 










e x e — 1, e x e = m, 


e x m = e, 


e x e 


= m x 


m—1, e x m = e 


Case 1: u = 0,8 (mod 16); 


= ±1, X 


= 1. 







Associativity relations: all associativity relations are trivial. 
Braiding rules: 

i?f = -1, 

pee pram /) 

iTj — JT^ — (7, 



pee _ 
pern 



1. 



1, 
0, 



-ft m 



-1, 
-1, 



Case 2 : z/ = ±4 (mod 16); = ±i, x 
Nontrivial associativity relations: 



■1. 



e e e 



\Y 



m m m 



m m m 

\Y 



£ e £ 



£ e £ 



£ m £ £ m £ 



e £ e 





mem mem 





Braiding rules: 

Rf = -1, 

pee pmm /) 

ix^ — ix^ — (7. 



pee pee pern pme a 



Table 2: Properties of anyons for z/ = (mod 4). 
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Superselection sectors: 1 (vacuum), e (fermion), a, a (vortices). 

Quantum dimension: d\ = 1, d £ = 1 d a = d a = 1; 

Topological spin: #x = 1, £ = —1, # a = 9 a = 9 = exp(|zi>) 

Frobenius-Schur indicator: K\ = 1, x e = 1. 



±l±i. 



Fusion rules: 

a x e = a, a x e = a, e x e = 1, 
Nontrivial associativity relations: 



axa = axa = e, a x a = 1. 



£ a £ 



£ O £ 



£ a £ 





a £ a 

£ 



a £ a 





Braiding rules: 



i?f = -1, 

T)aa Daa n 

K £ — H E — V, 



Daa Daa /)— 1 

iX^ — iX^ — (7 



Table 3: Properties of anyons for v = 2 (mod 4). 

9 The sixteen-fold way 

Let us again consider the theory with Z 2 -vortices and free fermions whose spectrum is gapped 
and characterized by the Chern number v. The properties of anyons in this model depend on 
v mod 16. In the previous section we studied the case of odd v\ now we assume that v is even. 

For even v, the vortices do not carry unpaired Majorana modes (see Appendix O) , therefore 
a vortex cannot absorb a fermion while remaining in the same superselection sector. Thus, 
there are actually two types of vortices, which are transformed one to another by adding a 
fermion. Two vortices of the same type may either annihilate or fuse into a fermion. In the 
first case the vortices are denoted by e and m; they obey the fusion rules (T4T|) . In the second 
case, we denote the vortices by a and a; the fusion rules are as follows: 



axe 



axe 



e x e 



a x a = a x a 



a x a 



1. 



(97) 
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Both sets of rules are Abelian, i.e., they do not contain "+" on the right-hand side. 

The associativity relations can be found along the same lines as in the non-Abelian case. 
First, the fermions must obey the trivial relation (1721) . The argument used in the derivation of 
Eq. ( 17T1) remains valid, but we obtain two relations instead of one. Note that they are not in- 
variant individually; for example, the sign in both relations may be changed to "+" . Finally, we 
solve the pentagon equation. The fusion rules with e and m admit two nonequivalent solutions, 
one of which is trivial; see Table [2j Both solutions are invariant under the transformation 

ivn^-ivn, IO--IO. ivc>~-ivc> m 

(where \ipz V is the basis vector in the corresponding fusion space), as well as trivial transfor- 
mations 

\r z y )^UyC l \r z y )- (99) 

In the case of a and a, there is only one solution, which doesn't admit any nontrivial symmetry; 
see Table |3] 

The braiding rules are found by solving the hexagon equations in conjunction with the 
requirement that R{ £ = — 1 and that transporting a fermion around a vortex is described by 
the multiplication by —1. Then we compute the topological spin of the vortex and identify each 
solution with a particular value of v mod 16. The results are summarized in Tables [2] and [3j 
(We have eliminated redundant solutions that can be obtained by the transformation (198]) ). 

Let us mention a couple of interesting properties of these Abelian theories. For v = 8 
(mod 16), all three nontrivial particles (e, m, and e) are fermions. For v = ±4 (mod 16), the 
particles e and m are semions with trivial mutual braiding. (In comparison, the well-known 
Kalmeyer-Laughlin state [69J has one semionic species.) 

10 Odds and ends 

What follows are some open questions, as well as thoughts of how the present results can be 
extended. 

1. Duan, Demler, and Lukin [12] proposed an optical lattice implementation of the Hamil- 
tonian (T4]). It would be interesting to find a solid state realization as well. For example, 
the anisotropic exchange could be simulated by interaction of both lattice spins with a 
spin-1 atom coupled to a crystal field. 

2. The weak translational symmetry breaking in the Abelian phase has some interesting 
consequences. A particularly unusual phenomenon takes place when the lattice has a 
dislocation. A particle winding around the dislocation changes its type: e m. Since 
m = e x e, the fermionic parity appears not to be conserved. To restore the conservation 
law, we must assume that the dislocation carries an unpaired Majorana mode. Therefore, 
Abelian phases can also be used for the implementation of quantum memory. 

3. Chiral phases {y ^ 0) require that the time-reversal symmetry be broken. In the present 
model, this is achieved by applying a magnetic field. However, a spontaneous breaking of 
time-reversal symmetry occurs in the presence of odd cycles in the lattice. For example, 
one can replace each vertex of the honeycomb lattice by a triangle. In this case, a gapped 
v = ±1 phase is realized without external magnetic field [TP] . 
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4. The representation of a spin by four Majorana fermions might be useful for other models, 
even though it does not lead to an exact solution. In particular, one can consider a 
variational mean-field state in which Cj are decoupled from bj. It is unclear whether 
this type of states occur in Heisenberg antiferromagnets. On the triangular lattice, such 
a state has larger energy than the 120° Neel-type state. It would be interesting if the 
former could be stabilized by additional interactions and quantum fluctuations. 

5. Topological phases with free fermions coupled to an effective Z2 gauge field have been 
classified by v mod 16. However, this analysis does not include multilayer systems. One 
can argue that if the interaction between the layers is weak enough, topological particles 
cannot tunnel between the layers. Thus, each layer is described by one of the 16 theories 
studied above. In mathematical terms, we have the direct product of several unitary mod- 
ular categories (UMCs). Strongly interacting layers are roughly described by n fermionic 
species interacting with (Z2) n -vortices, but a complete classification of such phases is yet 
to be found. 

6. In addition to the particle classification, the chiral central charge c_ = v/2 is an important 
robust characteristic, though its topological meaning is not so clear. It appears that in 
multilayer systems the total value of v is shared between the layers, i.e., increasing v by 
16 in one layer while decreasing it in another does make a different topological phase. 
More generally, a topological phase is characterized by a UMC and a real number c_ 
satisfying the relation (11721) . To prove or disprove this statement, a mathematical notion 
of equivalence between topological phases is necessary. It may be based on local (or 
quasilocal) isomorphisms between operator algebras. 

7. A related question is whether the space-time boundary of an arbitrary topological phase 
can be described by a two-dimensional conformal field theory (CFT) and when two such 
theories have the same topological content. Two CFTs may be considered topologically 
equivalent if there is a consistent theory for a one-dimensional boundary between them. 
Conjecturally, this is the case if and only if both CFTs have the same chiral central charge 
c_ = c — c and correspond to the same UMC. 

8. Another topic to study is Bose-condensation of vortices, which occurs when a vortex 
energy becomes negative due to some parameter change. The condensation of e-particles 
in the v = phase is equivalent to the confinement of m- and e-particles. Thus, the 
topological order is destroyed. In the v = 16 system, this process produces a phase 
without anyons or fermions, but with the chiral central charge c_ = 8. Under special 
circumstances, the boundary of this phase is described by the E 8 CFT at level 1, though a 
generic (nonconformal) perturbation drives it into a state with 8 chiral bosons propagating 
with different velocities. 

9. The condensation of ee-pairs in two adjacent v = 1 layers leads to the binding of single- 
layer vortices into 00 pairs, which are equivalent to a or a in the v = 2 phase. Thus 
the direct product of two v = 1 theories is related to the v = 2 theory. It would be 
interesting to have a general mathematical characterization of such relations. The notion 
of lax tensor functor (see Remark IE. 261 on page II II) can be useful in this regard. 
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Appendix A: Numerical results on the stability of the 

vortex-free phase 



The goal of this study is to compare the energy of the vortex-free phase with that of other 
phases. The CcLS6 J x — Jy — J z — 1 llcLS been investigated most carefully. The energy of 
the vortex-free phase equals E m —1.5746 per unit cell (i.e., per two sites, or one hexagon). 
The actual computation was done for tori with periodic or antiperiodic boundary conditions 
in each direction. The vortex-free phase and other periodic phases with small period are 
computationally very simple, so dealing with large tori is not a problem. However, the energy 
calculation for an arbitrary vortex configuration requires finding the singular values of an N x iV 
matrix, where N is the number of unit cells. With our setup (MATLAB on a PC) we could 
handle matrices of size iV ^ 2500, which corresponds to systems of linear size L = ^ 50. 
Although these numbers are pretty large, finite-size effects are still appreciable due to the 
gapless nature of the spectrum (see Section HJ). 




3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48 
L 

a) tori with basis (Lni, Ln 2 ) 




3 6 9 12 15 18 21 24 27 30 33 36 39 42 45 48 
L 

b) tori with basis (Lni, Ln 2 + n x ) 



Figure 11: Finite size effects in the vortex-free phase. 



Let us first discuss the finite-size effects in the vortex-free phase. From now on, we consider 
the relative energy, i.e., we subtract NEq from actual results. The plot in Fig. [TTa shows the 
relative energy as a function of size for symmetric L x L tori. The oscillatory behavior with 
period 3 is due the fact that e(q) vanishes at q = q* (see (1341) ). One may argue that each period 
r of the torus contributes ~ | r | 1 cos(2(q„,, r)) to the total energy (there are infinitely many 
such terms since the periods form a lattice). Interestingly enough, these contributions almost 
cancel each other for periodic boundary conditions on the torus with basis (Lni, Ln 2 + nj), 
where ni and n 2 are the basis vectors of the lattice (see figure in Eq. (|32l) ). The corresponding 
plot is shown in Fig. [TTb . 

The energy of an isolated vortex is E VOItex ps 0.1536 above the ground stated The calculation 
was done for tori with basis (Ljxi, Ln 2 + n%) for L = 9, . . . , 32. (We actually put 4 vortices 
on the torus of twice this size because the number of vortices must be even). Then the results 
were extrapolated to L = oo by fitting the curve E(L) = E vortcx + a\L~ x + d2L~ 2 to the data, 
separately for L = 3k, L — 3k + 1, and L = 3k + 2 (see Fig. [T2l) . 

9 No rigorous precision analysis was attempted, but the error in this and the other figures is expected not to 
exceed 1-2 units of the least significant digit. 
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Figure 12: Extrapolation to infinite size: 4 equally spaced vortices on the torus with basis 
(2Ln l5 2(Ln 2 + ni)). 
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0.074 
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0.025 
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0.046 
0.092 
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3 
4 


0.072 
0.096 



Table 4: Numerical results for J x = J y = J z = 1. Vortices are shown in gray; for periodic 
phases the unit cell is indicated by a parallelogram. 
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The energy of a vortex pair is smaller than 2E vortex if the vortices are close to each other: 
by « 0.04 for nearest neighbors and by « 0.07 for next-nearest neighbors. (We didn't try 
to compensate the finite size effects, so the precision is poor). These numbers suggest that 
inter-vortex interaction is rather strong and could, in principle, result in some configurations 
having negative energy. However, our further calculations give evidence for the contrary. 

We have checked all periodic phases with unit cell containing 1, 2, 3 or 4 hexagons, see 
Table HI (As mentioned above, this computation requires much less computer resources than 
the study of separate vortices or vortex pairs). In all these cases the energy is positive and 
increases as more vortices are added. The smallest energy per vortex is achieved by phases 1 
and 5 — 0.067 and 0.078, respectively. 

All 14 phases have positive energies (relative to the vortex-free phase) for all nonzero values 

of J xi Jyi Jz- 



Appendix B: Edge modes in phases B u 

It is understood that the edge spectrum depends strongly on particular conditions at the edge. 
The calculations below are only meant to illustrate the universal feature of the spectrum - 
the existence of a chiral gapless mode. 

Let us suppose that the honeycomb lattice fills the lower half-plane and is cut as follows: 




Row consists of the Majorana operators M, which would be decoupled from the rest of 
the system if not for the magnetic field. Performing the Fourier transform in the horizontal 
direction, we compute effective couplings between the rows as functions of q x . 
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An sin y, 



(100) 



■7 



Let us first find edge modes ignoring the operators b*j, i.e., the first row and column in the 
above matrix. If k = 0, then for 2-7r/3 < q x < Att/3 the matrix has a null vector with elements 
ip(2j) = and ip(2j + 1) = (— 2 cos yY i which corresponds to a zero-energy state localized 
near the edge. If k is not zero but still small, we get this spectrum: 



e(q x ) w 12k sin q x 



q x e [2vr/3, 4vr/3]. 



(101) 
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It is shown in Fig. [T3"a . assuming that k > 0, i.e., v = +1. The point where the energy curve 
crosses zero, q x = ir, corresponds to a left-moving gapless mode. 

Now we take the operators b* into account. If k — 0, then a zero mode exists for q x e 
[— 2n/3, 2n/3], whereas for q x G [27r/3, 47r/3] there are two modes with energies 

e{q x ) w ±7^1 -4cos 2 f . (102) 

The spectrum for k > is shown schematically in Fig. fT3b . In this case a left-moving gapless 
mode occurs at q x = 0. 

In conclusion, let us give a simple (but not completely rigorous) proof of Eq. (|56|) . The idea 
is due to Laughlin |71j: we put the system on a cylinder and run magnetic flux through it. For 
simplicity, we consider a cylinder of the smallest radius, which is equal to the lattice period 
divided by 2tt. Thus the Hamiltonian is given by the matrix iA(q x ); the variable q x plays the 
role of the magnetic flux. 

Let q x vary from to 2ir. When the energy of an edge stat^l \ip) crosses zero, the spectral 
projector P(q x ) changes by \tp)(ip\. For an arbitrary operator Q the quantity Tr(QP(q x )) 
changes by (ip\Q\ip). Let 

rx i -\s / -\ / 1 near the edge, . , 

Q jk = 9(j)6 jk , M = [ Q f ar away from the edge (103) 

(the exact form of the function g is not important). Then (ip\Q\ip) pa 1. Since Tr(QP(q x )) is 
a periodic function of q x , the abrupt changes must be compensated by a continuous variation, 
which we call "adiabatic" . Thus 

^edge= / I(q x )dq x , where I(q x ) = - Tr(QP), P = ( ^ ) . (104) 

J \dQx/ adiabatic 

10 We are using the first quantization formalism, therefore "states" are just superpositions of lattice points. 
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The adiabatic evolution of the spectral projector may be represented by the Heisenberg 
equation with a suitable Hamiltonian: 



P = i[H, P], where H = i[P, P]. (105) 

Indeed, P 2 = P, hence PP+PP = P and PPP = 0, implying that -[[P, P],P] = P. Therefore 

= -i Tr(QHP - QPH) = -i Ti{HPQ - HQP) = -iTr(H[P,Q]). (106) 

(These transformations are valid because each trace is represented by a finite sum.) The last 
expression in Eq. ( 1106[) may be calculated using the spectral projector for the bulk. Indeed, if 
the function g were constant, the commutator [P, Q] would vanish; thus the main contribution 
comes from the region where g changes from 1 to 0. 

In the final part of the proof, we take into account the periodicity of the bulk spectral 
projector, representing the site index j as (s, A). As is usual, lattice cells are indexed by s; 
it is an integer that increases in the positive y-direction. We may assume that the function g 
depends only on s and satisfies lim s _ >00 p(s) = 1, linis^-oo g(s) = 0. Then 



{[P,Q]) 



(g(s)-g(t))P t ,, 



sX ! 



^2(g(t + r)-g(t)) 



Using these auxiliary identities, we pass to the momentum representation: 
Tt(H[P,Q])= H sX ^P t ^ x {g(s)-g(t)) = ^ H QX ^ rfl P ^ rX ^(<?(t + r) - g(t)) 



s,X,t,fj, 



r,X,fj, 



J V 3qJ 2* 
Substituting the result into (11061) and (1104[) . we get 



HoX- rf _iPof_i,rX r 



Tr 



P. 



dP] dP\ dq y 
dq x \ dq y J 2n 



I(q x ) dq x 



—i 
2^ 



Tr 



P- 



OP 

dq x 



dP 
dq y 



dq x dq y 



v. 



Appendix C: Quasidiagonal matrices 

The goal of this appendix is to provide a formal argument for the existence of unpaired Major ana 
modes on vortices for odd values of the Chern number v. However, the developed formalism 
may be interesting on its own right. It suggests a rather efficient approach to problems like 
quantization of Hall conductivity in disordered systems without the use of excessively heavy 
tools. Some of results (in particular, the ones concerning the flow of a matrix and the Chern 
number) are actually poor man's variants of known mathematical theorems^ A powerful 
theory, called noncommutative geometry [45] was used by Bellissard at al [1] to prove rigorously 
that the Hall conductivity is quantized provided the electron are localized (see [721 173] for 

n Note for experts: we are effectively trying to build a K-theory on a manifold by considering functions that 
are not continuous, but rather constant on cells that are dual to simplices. This seems to be an awkward 
approach, but its possible advantage is the relation to a second-quantized case, see Appendix [Pi 
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n n n 

a) b) c) 

Figure 14: Cutting a Majorana chain: A4 V (B) = 1 in cases (a) and (c) whereas A4 V (B) = — 1 
in case (b). The dashed lines in (c) show a possible way to reconnect the broken pairs. 

another proof). We will use similar ideas, but focus mainly on providing the intuition rather 
than mathematical rigor. 

A physical example of a quasidiagonal matrix is an electron hopping matrix T = (tjk) on 
a (i-dimensional lattice. It is a Hermitian matrix with the property that is bounded in 
magnitude and vanishes if the distance \j — k\ between the sites j and k is greater that some 
constant L. Furthermore, if T has a spectral gap, i.e., if the eigenvalues are bounded away from 
zero, then the corresponding spectral projector P = §(l — sgnT) is also quasidiagonal. More 
specifically, the matrix elements Pjk decay exponentially with distance f^l 

In general, a quasidiagonal matrix is a lattice-indexed matrix A = (Ajk) with sufficiently 
rapidly decaying off-diagonal elements. Technically, one requires that 

\ A jk\ < c \j ~ k\~ a , a> d, 

where c and a are some constants, and d is the dimension of the space. Note that "lattice" 
is simply a way to impose coarse M. d geometry at large distances. We may think about the 
problem in these terms: matrices are operators acting in some Hilbert space, and lattice points 
are basis vectors. But the choice of the basis need not be fixed. One may safely replace the 
basis vector corresponding to a given lattice point by a linear combination of nearby points. 
One may also use some kind of coarse-graining, replacing the basis by a decomposition into 
orthogonal subspaces corresponding to groups of points, or regions in M. d . 

Let us outline the main results. We first consider quasidiagonal unitary matrices in one 
dimension and define an integral topological characteristic called flow. Then we study projection 
matrices in two dimensions. The Chern number v{P) of a quasidiagonal projection matrix P is 
expressed directly in terms of the matrix elements Pjk, see Eqs. ( 11241) and (11221) . This definition 
does not rely on translational invariance, but in the translationally invariant case we reproduce 
Eq. ([MD. 

After those preliminaries, we switch to questions related to Majorana fermions (all necessary 
background is given in Section[3]). From the mathematical point of view, we study quasidiagonal 
real skew-symmetric matrices B satisfying the condition B 2 = — 1. In one dimension, such a 
matrix is assigned a cutting obstruction Ai n {B) = ±1 with respect to an arbitrary cut rj dividing 
the chain into two parts. Let us explain the physical meaning of this number. The matrix B 
defines a pairing of Majorana modes (in the exact sense stated right after Eq. fT27|) . but we 
will use a cartoon description for illustration). If we cut the chain, the resulting pieces may 
carry unpaired boundary modes [64] . This happens when an odd number of pairs is broken, as 



12 Indeed, the function fix) = sgna; can be approximated by a sequence of polynomials p n {x) of degree n — ► oo 
which converges exponentially on the spectrum of T. Therefore \\P — p n (T)\\ < ab n for some b < 1 . On the 
other hand, the matrix elements (p n (T)) . fc vanish if \j — k\ > nL. 
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a) 



b) 



Figure 15: A map from two to one dimension: annular regions around the vortex (a) are 
collapsed into points on a ray (b). 

in Fig. [14b; otherwise one can modify the matrix near the cut so that to avoid broken pairs, 
see Fig. [14b. More formally, let us consider quasidiagonal real skew-symmetric matrices B' 
that agree with B at infinity but have no nonzero elements across the cut. We prove that if 
A4 V (B) = —1, then no such B' exists that satisfies the equation (B') 2 = —1. Note that while 
M. rj (B) depends on the cut, the relative Majorana number of two matrices on the same chain, 
M(A, B) = M v {A)Mr,{B) is invariant. 

The concept of Majorana number is also applicable in two dimensions. Let B 2 = — 1, 
and let us construct another matrix B' by inserting a Z2-vortex at the origin. More exactly, 
B'- k = ±Bjk, where the minus sign occurs when the link (j, k) crosses some fixed ray r (an 
analogue of the Dirac string). Note that the condition (B') 2 = — 1 is only true asymptotically 
(i.e., far away from the origin); it may or may not be possible to satisfy this equation by altering 
the matrix elements near the vortex. In fact, the presence of an unpaired Majorana mode at 
the vortex is defined as locally unrepairable failure of the equation (B') 2 = — 1. Such a mode is 
detected by an absolute Majorana number Ai(B'). For its construction, let us regard B and 
B' as one- dimensional by keeping track of the radial direction only: we divide the plane into 
concentric rings and map all sites in each ring to a single location on a ray (see Fig. [T5"|) . To 
simplify the calculation, we cut out the interior of a sufficiently large circle rj; let \t)\ be the 
number of sites removed. The presence or absence of an unpaired mode in this case is given 
by A4 V (B'), which is defined using some annular neighborhood of rj. Extrapolating to \rj\ = 0, 
we get M(B') = (—1)^M V (B') = M(B,B'). The last number only depends on the matrix 
elements of B in some neighborhood of the intersection point between rj and r. Finally, we 
show that Ai(B, B') = (—l) u( - p \ where v(P) is the Chern number associated with the projector 
P = |(1 — iB). Thus, vortices carry unpaired Majorana modes if and only if u(P) is odd. 

C.l The flow of a unitary matrix 

C.l.l Definition. Let us consider an arbitrary (possibly infinite) unitary matrix U = (JJjk)- 
We refer to the values of j and k as "sites" and define a "current" flowing from k to j: 





(108) 



j 
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Suppose that the sites are positioned on a line (more specifically, enumerated by integers), 
and that the off-diagonal matrix elements of U decay fast enough. We will see later that the 
condition 

\U jk \ < c\j - k\~ a , where a > 1, (109) 

is sufficient for our purposes. For now, let us simply assume that Ujk vanishes if \j — k\ is greater 
than some constant L. Then it is obvious that the total current through a "cross section" 77, 



j^rj k<r) 

does not depend on the choice of 77. The number J-(U) is called the flow of U. For example, 
the matrix with entries Ujk = 5j,k+s has flow s. 

C.1.2 Integrality of the flow. A nontrivial property is that the flow is quantized, i.e., 
has an integer value. To prove this statement, let us introduce the projector onto the sites M 
through N: 

u [m,n] _ f 1 ^ M < j = k < N, 
J [ otherwise. 

The projector onto an infinite interval is defined similarly. Let II = n' ' 00 ^; then 

F(U) = Tr(U*nU(l-n))-Tr(U*(l-IL)UIL) = Tr^UU - U). (112) 

(Caution: In the last expression, the order of operations is important: we first compute the 
matrix A = WHU — II, and then its trace, V. Am- ft i s n °t possible to use the cyclic property 
of the trace, Tr AB = Tr BA, since it is only valid if one of the matrices has a finite number of 
nonzero elements or, more generally, if the sum involved in the calculation of the trace converges 
absolutely.) 

Note that the matrices II and WHU are orthogonal projectors. If they were finite, their 
traces would be integers, and the difference would also be an integer. In the infinite case, one 
may refer to the general notion of the relative index of two projectors [731 13] • However, we 
will proceed in a more pedestrian fashion and simply truncate the matrices. 

Let A = WHU — II. It is clear that the matrix element Ajk vanishes if \j\ or \k\ is greater 
than L, so A is not changed by the truncation to the interval [-L, L\. Hence 

T{U) = Tr(C/tnC/) truilc - Tr n trunc , where A tIWC = n^Unl-^. (113) 

If A is an orthogonal projector and commutes with Il[ _L ' L l, then A trunc is also an orthogonal 
projector. Obviously, II and A commute with Ii^~ L,L \ and so does the matrix WHU = II + A. 
Thus both n trunc and [U^HU) tranc are orthogonal projectors. It follows that J-(U) is an integer. 

In equation (I112p . the projector II may be replaced by another operator with the same 
asymptotics at infinity. (In other words, we can "blur" the boundary between the left and the 
right half-line.) Indeed, the expression Tr(WAU — A) vanishes if the cyclic property of the 
trace is true, i.e., if the elements of A decay at infinity fast enough. Adding such a matrix A 
to II will not change the result. Thus 

(5 .i for 7 k ^ ~~ \~ oo 
n , , ' (H4) 

tor j,k — > — oo. 
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(Here we also assume that the matrix Q is quasidiagonal, i.e, \Qjk\ is bounded by a rapidly 
decaying function of \j — k\, cf. (11091) .) 

C.1.3 Translationally-invariant case. The flow can be easily calculated if U possesses a 
translational symmetry. Let us group the sites by unit cells and index them as (s, A), where s 
is the number of the cell, and A refers to a position type inside the cell. We assume that U s \^ 
depends only on t — s, A and fi. In this case, we can express J-{U) in terms of the trace per 
unit cell (denoted by tr) and a position operator X: 

T{U) = tr(t/t[x, U\) , where X sX ^ = s 5 st 8 x „. (115) 

(For a proof, replace the factor s in the definition of X by the function g(s) such that g(s) = s 
for \s\ ^ N and g(s) = iVsgns for \s\ > N, where N is large; then use Eq. (11141) ). 

In the momentum representation, the operator U becomes U^q) = Ylt e ^ q '*^oA,t/i- We 
may also use these simple rules: 

[X,A] - i^, trA^ T^TrA (116) 
dq J-* 2 '* 



Thus 



H U) = I />(^f ) ^ = -± rjMetUiq)). 



(117) 



C.2 General setting 



The result about the flow quantization can be extended to matrices satisfying condition ( 1109}) . 
which guarantees that the sum in Eq. (II 101) converges absolutely. 

Another generalization is coarse-graining, which means to allow multiple states per site. The 
sites are now described by orthogonal subspaces of the total Hilbert space. This generalization 
is useful if want to treat a two-dimensional system as one- dimensional. 

Finally, we may try to apply the notion of flow to finite systems. Let us assume that the 
unitarity condition, £\ U^Ujk = ^ • U^U^ = 5ik, is approximate and holds only if j and k 
belong to some interval. In this case, one needs to restrict the sum (II 101) to the same interval, 
and the result will not be an exact integer. But one can actually give an upper bound for the 
deviation from the closest integer. 

Setting the last generalization aside, let us put the integrality theorem into a rigorous form. 
(The reader may safely skip this formalism and proceed to the next subsection.) 

Theorem C.l. Let U be a unitary operator acting in the Hilbert space 7i = ©"l.^Tij and 
represented by a matrix whose entries Ujk are linear maps from TCk to Hj. Suppose that the 
off-diagonal entries are Hilbert- Schmidt operators satisfying the condition 

HMhs ^c\j-k\- a (118) 

where c and a > 1 are some constants, and \\ ■ ||hs is the Hilbert- Schmidt norm. Then the sum 
in the expression for the flow, 

= EE( Tr (^^) - ( 119 ) 

i^O A:<0 

converges absolutely and has an integer value. 
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Proof sketch. Let 7Y_ = ® J<0 7ij and 7i + = ® J>0 7~(?- An operator A is said to be quasi- 
trace-class relative to the decomposition Ji = Ti- ©7^+ if the matrix elements A ++ and A are 

trace class, whereas A + _ and A^ + are Hilbert-Schmidt. By definition, Tr A = Tr A ++ + Tr A 

An operator B is said to be quasi-bounded if B ++ , B are bounded and B |_ are 

Hilbert-Schmidt. Both types of operators form Banach spaces with respect to suitable norms; 
quasi-trace-class operators form an ideal in the algebra of quasi-bounded operators. 

One can show that the operator A = WHU — II is quasi-trace-class. Also, the operators 
P(L) = U^ L ' L ^WUUIl^ L ' L ^ are almost projectors, in the sense that lim L ^ 00 (P(L) 2 -P(L)) = 
with respect to the quasi-trace norm. Hence P(L) can be approximated by a projector, which 
has an integer trace. □ 

We omit the details and proceed to higher dimensions. Generally, our style will not be very 
rigorous, but one can hopefully elaborate the results using the following definition. 

Definition C.2. Let a positive integer d (the dimension) and a real number a > d be fixed. 
A matrix (Ajk : j,k G Z d ) with operator entries is called quasidiagonal if all its off-diagonal 
entries are Hilbert-Schmidt, and there are some constants c and d such that 

IIAtfll < c ' II4>*IIhs < c\j - k\~ a for j + k. (120) 

Note that if A is unitary, then \\Ajj\\ ^ 1 automatically. Quasidiagonal matrices form a 
Banach algebra with the norm 

P|| qd = sup ||A,j|| +7sup|j-A;| a ||A ifc ||Hs, (121) 
where 7 is a sufficiently large constant. 



C.3 Chern number 

C.3.1 General definition of v(P). Let P be an orthogonal projector represented by a 
quasidiagonal matrix in two dimensions. (For simplicity, we assume that the matrix elements 
Pjk are scalars, though they may be operators as well.) For each triple (j, k, I) we define a 
"2-current"0 

hjki = 127ri(PjkPkiPij — PjiPikPkj) ■ (122) 

It is clear that hjki is antisymmetric in all three indices. Since P is Hermitian, hjki is a real 
number. Moreover, since P 2 = P, 

(dh) kl d ^j2 h ^ = °- ( 123 ) 

Let us partition the plane into three sectors and define a quantity f(P), which will be shown 
to generalize the notion of the Chern number: 



v{p) = h(A, b, c) = J2J2H h i»- ( 124 ) 

jeA keB leo 



13 Another name of this object is "simplicial 2-chain". More exactly, the 2-chain is the formal sum 
Tlij<k<l hjkiAjki, where Ajki is a combinatorial simplex. 
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This sum converges absolutely. On the other hand, its value does not change if one reassigns any 
site from one sector to another (due to Eq. fl 1 2 3 [) ) . Therefore v{P) is a topological invariant: it 
is constant, provided that A, B, C are arranged in the counterclockwise order. We will see that 
v(P) is actually an integer, and its value agrees with Eq. ( |54l) in the translationally invariant 
case. 

C.3.2 Some properties and the translationally invariant case. Let us rewrite Eq. (11241) 
in an operator form using the relation A + B + C = 1, where the symbols A, B, C designate 
the projectors on the corresponding sectors: 

u{P) = 12m(Tr(APBPCP) -Tr(APCPBP)) 

= Ani Tr(APBPCP + B PC PAP + C PAP BP - APCPBP - CPBPAP - BPAPCP) 

= Ani Tr (PAPBP - PBPAP) = Am Tr [PAP, PBP] = v{P, A, B). (125) 

Note that the value of z/(P, A, B) will not change if we add arbitrary finite matrices to A and 
B. Moreover, the commutator [PAP, PBP] is nonzero only in the region where the supports 
of A, B, and C = 1 — A — B touch each other, hence we can make arbitrary changes away 
from the triple contact point. Thus the operators A, B, C do not have to be projectors or 
even commute with each other; it is only important that they are supported by regions with 
the particular topological configuration (possibly with some overlap along the boundaries and 
triple overlap at the center). 

This freedom in the choice of A and B is very useful. As one application, we show that v(P) 
is additive. Let Pi and Pi be projectors onto orthogonal subspaces (i.e., P1P2 = P2P1 = 0) and 
let P 3 = 1 - Pi - P 2 . We may replace A, B, and C = 1 - A - B with 

333 

a' = J2 p kAP k , B ' = Y1 PkBPk > C ' = J2 PkCP ^ 

k=l k=l k=l 

which differ from A, B, and C only at the boundaries of the corresponding regions. Then 
(Pi + P 2 )A'(P 1 + P 2 )B'(P 1 + P 2 ) = P 1 AP 1 BP 1 + P 2 AP 2 BP 2 , 
u(P 1 + P 2 ,A,B) = u(P 1 + P 2 ,A',B') = u(P 1 ,A,B)+u(P 2 ,A,B). 

Thus 

v{P x + P 2 ) = v{P x ) + v{P 2 ), 1/(1 - P) = -v{P). (126) 

Now, we calculate z/(P, A, B) for a different topological topological configuration. Let = 
iT 1 ) + iT 4 ) and = IT 2 ) + IT 1 ) be the projectors onto the right and the upper half-plane, 
respectively (see the picture below). Then 

....J... u{P, U {x \U (y) ) 

MiUdl x = u(P, u^\u^) + u(p, u^\u^) + u(P, nw, n< 2 >) + u(p, n< 4 >, n«) (127) 
::: 3 ::::1:i = v(P) + + + v{P) = 2u(P). 

It follows that 

u(P) = 2vriTr[Pn (a;) P, Pn (y) P] = 2vri Tr(^P[[n (a;) , P], [Tl {y \ P}] + P[n (x) , n fe) ]p) , (128) 
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where the last term vanishes because 11^ and 11^ commute. Of course, H^ x \ H w may be 
replaced by any topologically equivalent pair of operators, i.e., we may deform or blur the 
boundaries of the corresponding half-planes. 

In the translationally invariant case, the right-hand side of Eq. ( 1128}) can be calculated by 
widening the fuzzy boundaries so that they turn into linear functions. Thus 

u(P) = 27utr(p[[X,P],[r,P]]), (129) 

where X and Y are the operators of x- and ^-coordinate (resp.), and tr is the trace per unit 
cell. Passing to the momentum representation and using a two-dimensional analogue of the 
correspondence (I116p . we recover Eq. (1541) . 

C.3.3 The integrality of v{P). We have seen v{P) has topological nature; let us show 
that it is an integer. Let 

q{x) = pn^P, Q(») = Pn fe) P, U {x) = exp(2mQ {y) ), U {y) = exp(-2niQ {x) ) . (130) 

The operator coincides with the identity matrix away from the x-axis, and is equal 
to the identity away from the |/-axis. We may proceed in two different ways. 

First proof. Let us regard the two-dimensional lattice as one- dimensional by identifying all 
sites in each vertical row. Then U^ x ' satisfies the hypothesis of Theorem IC.ll hence it has an 
integer flow. But the flow can be expressed using Eq. (P33J with Q = Q(*>. Thus 



T{U^) = Ti{e~ 2mQ(v) Q^e 2mQ{v) -Q^) = jf \r (J- ( e -W»'QW e W»)jj 
= 27riTr[Q ix \Q (y) ] = u(P). 



(131) 



Note an analogy with the Laughlin argument: the integration over ip corresponds to the adia- 
batic insertion of a magnetic flux quantum. □ 

Second proof (sketch). Let us employ the notion of Fredholm determinant. It is a general- 
ization of the determinant to infinite matrices that are close to the identity, specifically of the 
form 1 + K, where K has a well-defined traced! 



det(l+iT) d = exp(Trln(l + K)) = l + TiK+ (~ Tr K 2 + ^(Tr Kf^j + .... 



(132) 



The exact meaning of this definition is as follows: we first obtain a formal expression as a power 
series in K and then evaluate each term. (If K is a matrix of size n, then all terms of degree 
higher than n vanish.) We claim that 

det(u {x) U {y) {U {x) Y\u {y) y l ) =!■ ( 133 ) 

Indeed, the identity det (ABA~ 1 B~ 1 ) = 1 is true at the formal power series level, assuming 
that the cyclic property of the trace holds for any product of A — 1, A -1 — 1, B — 1, P -1 — 1 



14 



It is important that any matrix product that occurs in the problem has a trace satisfying the cyclic property. 



56 



containing at least one factor with A and at least one factor with B. On the other hand, one 
can show that 



(») iw-OOO 



det e l ^ Q 



)(«) 



exp(^Tr[gW,Q^]). 



(134) 
□ 



C.4 Majorana numbers 

In this section we study real skew-symmetric matrices satisfying the equation B 2 = —1. Re- 
call that such a matrix defines the ground state of a Majorana system with a quadratic 
Hamiltonian (cf. Eqs. ( 1261) and ( 1271) ). First, we consider finite or general infinite matrices (i.e., 
work in dimension zero), then proceed to quasidiagonal matrices in dimension one, and finally 
to dimension two. 

C.4.1 Dimension zero: fermionic parity and the PfafRan. A few more words about 
the Majorana system, then we will turn to matrices. Let us group the sites (or the corresponding 
Majorana operators) into pairs (2k — 1, 2k) to form "full" fermionic modes with the occupation 
number operators a\ak = |(1 + ic2k-iC2k)- Note that the operator —ic 2 k-iC2k has eigenvalue 1 if 
the mode is empty and —1 if it is occupied. Thus the total fermionic parity is characterized by 
the operator Ylk(~^k-iC2k) ■ The fermionic parity of the ground state is equal to the Pfafflan 
of the matrix B: 

\[{-tC2k-iC2k)m = (Pf B)|tf>, PiB = ±1. (135) 

k 

Let us recall some standard definitions. The Pfafflan of a skew-symmetric matrix is ex- 
pressed as a sum over all partitions of the set {1, . . . , 2N} into pairs, or over all elements of the 
permutation group S2N- 

Pf A = _. T „. ^ Sgn(r) A r (i) )T ( 2 ) • • • A r (2AT-l),r(2JV) ■ (136) 
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For example, 
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Pf 



/ a 12 
— a 12 

—0-V3 — Ct23 

y— <2i4 — 0-24 



«13 

a-23 


— a.34 



0-2A 
034 



^12^34 + Ol4«23 — 0-130-24- 



The Pfafflan satisfies the equations 
(PfA) 2 = det A, 



Yi(WAW T ) = Pf(A) det(W). 



(137) 



We will also need a generalization of the Pfafflan to infinite matrices, which is called Fredholm 
Pfafflan [TH]. Let A and B be real skew-symmetric, B invertible, and A — B have a well-defined 
trace. Then 



Pf (A, B) = y/det(l + K) = exp Q Tr ln(l + K)\ 

: (TrK) 



1 + -TrK + 
2 



1 'IV A - • ' - 

4 



(138) 



+ 



where K = AB 



1. 
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This definition is understood as that of the Fredholm determinant: we first obtain a formal 
expression as a power series in K and then evaluate each term. (For finite matrices of size 2N 
we get Pf(A, B) = (Pf A) (Pi -B)" 1 ; note that the terms of degree higher than iV vanish). The 



Fredholm Pfaffian has the following properties: 

Pi{B,B) = l, Pi(A,B)Pi(B,C)=Pi(A,C), (139) 

Pi(A, Bf = de^AB" 1 ), Pi(WAW T , V BV T ) = Pi(A, B) det(y _1 W); (140) 

HA 2 = B 2 = -1, then Pi(A, B) — ±1. (141) 



C.4.2 Dimension one: cutting obstruction and the relative Majorana number. 

Now, consider a matrix that describes the ground state of Majorana fermions on an infinite 
chain. From a formal point of view, it suffices to say that B = (Bjk : j,k G Z) is quasidiagonal, 
real skew-symmetric, and that B 2 = —1. Suppose that we cut the chain between sites rj — 1 
and r\. Is it possible to find a new matrix B' that contains no cross elements (i.e., has the form 

B'_ © B' ++ ), satisfies the condition (B') 2 = —1, and coincides with B far away from the cut? 

The following quantity represents an obstruction for the construction of such a matrix: 

M V (B) = Pi{V(q)BV(q), B) = ±1, where V(q) = 1 - 2n [ "' oo) . (142) 

Recall that n^' 00 ) denotes the projector onto sites i] and higher; (V (r])BV (rf)) . k = ±Bjk, where 
the minus sign occurs when j and k lie on different sides of the cut. 

We claim that the conditions on B' cannot be met if M. V {B) — — 1. Indeed, 

M,{B) M V {B)- 1 = Pi (V( 7 )W( 7 ), V(r,)BV(r,j) 

= det(y(ry)- 1 y( 7 )) = (-1)^. 
The same is true for B', hence the relative Majorana number, 

M(B,B') = M 1 (B)M 1 (B') (144) 

does not depend on 7. However, Ai 7 (B') = Ai 1 (B) if | 7 — 77 1 is large enough (since B and B' 

agree far away from the cut), while A4 V (B') = 1 (because B' = B'_ © B' ++ ) and Ai v (B) = —1 

(by assumption). We have arrived at a contradiction. 

From the physical point of view, the impossibility to construct a suitable matrix B' indicates 
the presence of unpaired Majorana modes on both sides of the cut. This situation may be 
described by matrices B'_ and B' ++ having one zero eigenvalue. 

The cutting obstruction and the relative Majorana number can also be expressed as Fred- 
holm determinants. To this end, consider the quasidiagonal matrix 

W = V(0) exp(TrX), where X = ^(ll {x) B + BU^, U {x) = U [0 '°°\ (145) 

Clearly, X is real skew-symmetric and commutes with B. Hence W is orthogonal, and WBW T = 
V(0)BV(0). On the other hand, X is equal to B at +00 and vanishes at —00, therefore W is 
equal to 1 at ±00. One can actually show that W — 1 has a well-defined trace. It follows that 

M (B) = Pf (WBW T , B) = detW. (146) 

Similarly, one can define X', W' as functions of B' . Thus 

M(B, B') = det(W- 1 W) = det(e~ wX e nX '). (147) 
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C.4.3 Applications to two dimensions. Eq. (11471) may be used as a more general defini- 
tion of the relative Majorana number, which works even if J\4 V (B) and Ad^B') do not exist. For 
example, let B be a quasidiagonal matrix in two dimensions and B' — (1 — 211^)5(1 — 2n^). 
In other words, B'- k = ±Bjk, where the minus sign occurs when j and k lie on different sides 
of the x-axis. If we regard the system as one- dimensional (collapsing it in the y direction), 
the Majorana number M(B, B') is well-defined and can be expressed in terms of the projector 
P = |(1 — iB). The calculation follows. 

First, we obtain an expression for the operator X: 



x = ifpu^p - (1 - P)n w (l - P)Y 



Then we use the relation 1 — 211^ = exp(±i7rll^): 



M(B, B') = det(e^ x e- mmy 'e 
Finally, we compute the trace of the commutator: 



exp(m 2 Tr[X,n {s/) ]). 



Tr[X,IT (y) ] = Tr((P 2 + (l-P) 2 )[X,n (3/) ]) = Tr (p [X, U {y) ] P + (1 - P) [X, U iv) ] (1 - P) 



Pii^p, pu^p 



— iTr 



;i - p)n w (i - p), (i - p)n (j/) (i - p) 



Uv{P) - 1/(1 - P)) = -u(P), 
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because v{\ — P) = —v{P). Thus 



M(B,B') = (-l)"( p l 



(148) 



Equation (I148j) can be applied to the geometry shown in Fig. [151 In this case, B and B' 
correspond to two different Z 2 -field configurations: the first is regarded as vortex-free and the 
second has a vortex at the origin. The specific assumptions are as follows: (i) B and B' are 
real skew-symmetric quasidiagonal matrices in two dimensions; (ii) B and B' coincide, except 
that B'- k = —Bjk for those links which cross the ray r; (iii) B 2 = —1 (therefore {B') 2 = —1 far 
away from the vortex). Note that Ai v (B) and M^B') may be defined for a sufficiently large 
loop 7] enclosing the vortex. Let 1 77 1 be the number of sites inside the loop. Then the absolute 

Majorana numbers M(B) = f (-l) |r?l M V {B) and M(B') d = (-l) |r?l M V {B') do not depend on 1] 
and indicate the presence of unpaired Majorana modes in B and B', respectively. However, an 
unpaired mode cannot exist in B since B 2 = —1; therefore Ai(B,B') = Ai(B'). We conclude 
that 

M(B') = (-ir (p) . (149) 



Appendix D: Some remarks on the chiral central charge 

This appendix is an attempt to understand the physical and mathematical meaning of the 
chiral central charge beyond the CFT framework. Recall that chiral central charge is just the 
coefficient c_ in the edge energy current formula floTj) . This definition will be refined in Sec. ID. II 
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It turns out that the edge energy current is a property of the bulk ground state. The edge current 
is related to a bulk 2-current. Although the theory is in its embryonic stage, it looks like a 
second-quantized version of the Chern number for quasidiagonal matrices (cf. Sec. IC.3j) . 

In Sec. ID. 21 we discuss modular transformations of the partition function on a space-time 
torus. In particular, the phase factor e - 2mc -/ 2A appears in the description of the Dehn twist 
along a time circle. However, other modular transformations are difficult to define because 
space and time are not physically equivalent. 



D.l The edge energy current and a bulk 2-current 

D.l.l Energy current in the Hamiltonian formalism. Let us represent the Hamilto- 
nian as a sum of local terms: 

H = J2 H i> ( 15 °) 

3 

where Hj is a Hermitian operator acting only on spins in some neighborhood of the point j. 
Note that the decomposition into local terms is not unique. We also assume that the ground 
state \^f) is separated from excited states by an energy gap; thus equal-time spin correlators 
decay exponentially with distance p2]. We may slightly extend the notion of locality so that 
Hj acts on all spins but 

|| [H jt ofl || ^ u(\j-k\), (151) 

where u is some function with fast decay at infinity (faster than any power). Then the equal- 
time correlators also decay faster than any power, provided the spectrum is gapped. 

Microscopic energy current can be defined for any system with Hamiltonian (11501) at finite 
temperature. It is convenient to fix the temperature at T = 1 and vary the Hamiltonian instead. 
The thermal average of an operator X is defined in the standard way: 

(X)=Tr(pX), where p = Z^e' 11 , Z = Tre~ H . (152) 

According to the Heisenberg equation, 



3 



dt 



-i[Hj,H] = y^/jfc, where f jh = -i[Hj,H k ]. (153) 



Thus the operator fj k characterizes the energy current from k to j. We denote its thermal 
average by f jk : 

f jk = f(H j ,H k ), where f(A, B) = f (—i[A, B]). (154) 

(Note that this definition depends on the decomposition of the Hamiltonian into local terms.) 
Of course, the energy current is conserved: 

(df) k = J2fik = °- ( 155 ) 

3 

Indeed, 

(df) k = f(H,H k ) = Tr(-ip[H,H k ]) = Tr(-i[p,H]H k ) = 0, 
since p and H commute. 
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Let us try to apply the general formula (11541) to a specific geometry. Let if be the Hamil- 
tonian of an infinite two-dimensional system with a ground state |^) and gapped excitations. 
To simulate an edge, we introduce another Hamiltonian, 



Ef if y(j) — > +00, 

Hj, Hi =1 (oo) (156) 



where y(j) is the y-coordinate of the point j. Thus the lower part of the system is "bulk 
material" characterized by the state |^ r )(^ r |, whereas the upper part is "empty space" (actually, 
a set of noninteracting spins in the maximally mixed state). By definition, the chiral central 
charge is 12 jn times the total energy current in the negative x-direction. 

D.1.2 Representing the ground state by local constraints. At this point, we have 
encountered a technical difficulty: the energy current does not necessarily vanish in the bulk. 
The problem can be avoided if we describe the ground state by a set of local constraints. (A 
local constraint is a local operator that annihilates the ground state.) Hamiltonians that are 
explicitly written as sums of local constraints, e.g., the Rokhsar-Kivelson model [76j [19], have 
provided a lot of insight into properties of quantum many-body systems. 

Proposition D.l. Any gapped local Hamiltonian H = Hj whose ground state \^) has zero 
energy can be represented as a sum of Hermitian local operators Hj such that Hj\ty) =0. 

Proof sketch. Without loss of generality we may assume that (^\Hj\ty) = 0. (If not, change 
Hj by a constant.) Let all excited states have energy greater or equal to A. Choose a smooth 
function w(e) such that 

w(-e) = w(e)*, w(0) = 1, w{e) = for |e| ^ A. 

Its Fourier transform, w(t) = f w(e) e~ £t || decays faster than any power of t though more 
slowly than e -7 '*' (because w is smooth but not analytic). We set 

POO 

Hj = / e iHt H 3 e- iHt w{t)dt. (157) 



00 



It is easy to check that Ylj Hj = H and Hj\ty) = 0. A locality condition of the form (11511) can be 
established using the bound of Lieb and Robinson [77] on correlation propagation (cf. [IS]). □ 

For our purposes, we can actually use an arbitrary gapped Hamiltonian X = £\ Xj made 
up of local constraints for the given state It turns out that the chiral central charge depends 
only on the ground state. Indeed, one can interpolate between different sets of local constraints 
for the same state. Consider the situation where the constrains are Xj in close proximity of 
the edge, gradually changing to X'- deeper in the bulk. No energy current is associated with 
the transition between X and X' . 

As an aside, we conjecture that there is a set of local constraints Yj such that the Hamiltonian 
Y = J2j YjYj i s dipped. Thus excitations can be efficiently detected locally, by coupling each 
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Yj to an external indicator: -^detector = J2j{ Y j ® +Yj ® A stronger version 

of this conjecture asserts that the original Hamiltonian H can be represented as JE Y } Y v where 
Yj are local constraints. This is true for the quadratic fermionic Hamiltonian ( fl8|) . which needs 
to be offset by its smallest eigenvalue (1251) so that the ground state energy becomes zero. It 
is easy to see that the matrix \iA\ is symmetric and that JE fc (|z\A|) - k CjC k is equal to Tr \iA\ 
times the identity operator. Hence 

H = ^2(iA jk + (\iA\).^j Cj c k = J2 D ikCjC k , where D = ^(iA+\iA\) ^ 

j,k j,k 

(meaning that D is positive semidefinite) . The matrix \/D is quasidiagonal, provided A is 
quasidiagonal and has a spectral gap around zero. Thus we can represent the Hamiltonian in 
the desired form: 

H = y^ Y l Y m, where E m = V(Vd) c k . (158) 

m k 

D.1.3 Bulk 2-current. In the disk geometry (i.e., when interaction is strong near the 
origin and vanishes at infinity), the chiral central charge is 12/7T times the energy current in 
the counterclockwise direction. Let us divide the disk into three sectors, A, B, and C in the 
counterclockwise order (see picture in Eq. (11241) . or in Eq. (11601) below). In this setting, the 
energy current is EfceB E/eA 

Let us try to represent the chiral central charge c_ in a form that would explicitly depend 
on the bulk Hamiltonian rather than the edge. To this end, we construct a 2-current h such 
that 

f = dh, where (dh) kl = h jkl . (159) 

i 

(A 2-current is an antisymmetric function of three lattice sites that decays fast enough as the 
distance between any two sites increases.) Then in the three-sector geometry we have: 

EE a. = EEfe'v^ = -EEE"* 

keB /eA keB /eA \ j / jeA k£B lec 

Thus, 



c = --£EE*i«- (160) 

j'eA keB leC 

The last sum is dominated by a neighborhood of the triple contact point between A, B, and C . 

Unfortunately, there is no canonical expression for hjki- Instead, we can define a canonical 
object g which is a 2-current on the lattice and a 1-form on the space of local Hamiltonians. It 
satisfies the condition 

dg = df. (161) 

A suitable 2-current h is obtained by integrating g over an arbitrary path H(/3) : /3 £ [0, oo) 
in the space of spin Hamiltonians. Here we assume that H(/3) is a sum of local terms Hj(/3), 
the corresponding thermal state has short-range correlators, H(0) = 0, and H(f3) ~ (3H^ as 
f3 — > oo. Thus we implicitly accept the following conjecture. 
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Conjecture D.2. The ground state of a gapped local Hamiltonian on a two-dimensional lattice 
can be transformed into the maximally mixed state without a phase transition. 

(In this formulation, (3 changes from oo to 0.) Note that usual symmetry-related phase tran- 
sitions are not a problem because H(/3) is not required to be symmetric. Rather, we should 
worry about a transition from a topological phase to a non-topological one. However, in two 
dimensions the long-range topological order in the ground state is destroyed at any finite tem- 
perature. Indeed, all topological excitations are point-like and not confined, therefore they 
have a finite density at any non-zero temperature. This argument does not apply to three- 
dimensional systems, where string-like excitations are possible. 

Let us give an explicit formula for g. For any two operators A and B, we define the 
Matsubara time-ordered average, as well as its "truncated" version 1^1 

(A(t)B(O)) = Z- 1 Ti{e- (1 - T)H Ae- TH B), where Z = Tie~ H ) < r < 1; (162) 

((A(r)B(O))) = (A(r)B(O)) - (A)(B). (163) 

Next, we define a function of three operators: 

fi(A,B,C) = i [ ((A(r)[B,C](0)))dr. (164) 
Jo 

Finally, we assume that the Hamiltonian H and the operators A, B, C depend on some pa- 
rameters, so that we can differentiate them. Now we can define g: 

g(A, B, C) = fi(dA, B, C) + fi{dB, C, A) + y,{dC, A, B), g jkl = g(H v H k , Hi). (165) 

Let us verify Eq. (|IoT]) . Note that fx(X,Y,H) = -i([X,Y)) because (X(t) [Y,H)(0)) = 
±{X(t)Y(0)) and (\Y,H]) = 0. We proceed as follows: 

dg(B,C) = g(H,B,C) = fi(dH, B, C) + fi(dB, C, H) + fi(dC, H, B) 

= (-iZ- 1 Tr (de- H [B,C\) - i(dH)([B, C])) - i([dB, C\) + i([dC, B}) 

= d{-iZ- 1 Ti{e- H [B,C})^ = df{B,C). 
D.2 Modular transformations beyond CFT 

The modular invariance [22] is usually formulated in the CFT framework. In this section, we 
demonstrate some modular properties of the partition function of edge excitations, not assuming 
conformal, rotational, or any other symmetry. First, let us generalize the transformation of the 
vacuum charactei0 under the Dehn twist: 

Xi{w + 1) = e- 2 ™-/ 24 Xl (w), (166) 

15 This notation is reminiscent of the formula A(r) — e rH Ae~ rH , but the latter may be problematic to use 
because e rH w exp(T/3iJ^ 00 - ) ) diverges as /3 — > oo. Thus the expression (A(r)i?(0)) should be viewed as a whole. 

16 Some care should be taken to make a proper connection between the physics of edge modes and the the 
CFT formalism. If all edge modes have the same chirality, characters correspond to representations of some 
chiral algebra (as is usual in CFT). In general, the vacuum character xi is defined as a sum over all fields with 
trivial monodromy whereas \a includes all fields that occur in the two-dimensional theory when an anyon of 
type a emerges on the surface (see Fig. [9Jd on page [29]) . 
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where w is the modulus of the torus and the subscript 1 refers to the vacuum sector. 

In fact, the derivation of Eq. (11661) does not require much more than the knowledge of the 
energy- momentum tensor, or just its rx-component, which is exactly the energy current /. 
(We use imaginary time, r = it.) In the thermodinamic limit, the coordinate transformational 
(x,t) — > (x, r + £(x)) induces the action change 5S = il J ^dxdr. If the imaginary time 
is cyclic, then twisting it by a full circle results in the action change 5S = and the 

multiplication of the partition function by e~ s = e - 27nc -/ 24 . We now describe a particular 
setting in which this result is applicable. 

Let us consider the partition function Z = Tr e~~ H l T of a spin system on a disk. As is usual, 
we think of the evolution over the time period 1/T, making the edge into a torus and the disk 
itself into the three-dimensional manifold M = D 2 x S 1 . Let us suppose that the temperature 
T is much smaller than the energy gap in the bulk, so that Z is mostly determined by edge 
excitations. In fact, the partition function depends on full detail of the edge Hamiltonian as 
well as some conditions in the bulk. For example, we may place an anyon of type a somewhere 
in the disk, which slightly changes the edge excitation spectrum. In other words, we insert a 
particle worldline I into the manifold M, or act by the time-like Wilson loop operator W a (£). 
Thus a partition function Z a is defined. If all edge modes propagate with the same velocity v 
and can be described by a CFT, then Z a = x (— 1/iu), where w = iLT/v. In the limit w — > oo, 

Z a = Xa(-l/w) = ^2S ab Xb(w) « S a iXi(w), (167) 

b 

because all characters Xb{w) for b ^ 1 are exponentially smaller than the vacuum character 
Xi(w). 

The following argument does not depend on CFT. As already mentioned, twisting the torus 
by £ along a time circle leads to the the action change by iI£/T, or the multiplication of the 
partition function by e~ ll ^l T . (The sign in the exponent depends on the direction of the twist; 
we choose it to be consistent with the CFT formula.) For the full twist, £ = 1/T, we get this 
equation: 

Z' a « e~ 27Tic - /24 Z a . (168) 

Note that the finite-size quantization of edge modes might affect the energy current. We 
therefore assume that the disk circumference L is much larger than the edge correlation length 
^edge ~ V /T, where v is the maximum group velocity of any excitation in the systemic Small 
corrections proportional to exp(— L// e d g e) are expected. Thus, Eq. (I168P and its estimated 
precision are in perfect agreement with the CFT formulas (I166P and (I167p . 
With some more work, one can also derive an analogue of this CFT result: 

Xb(w + 1) = e- 2 ™-/ 2 % X bW, (169) 

17 For lattice models, one should consider the transformation (j,r) — > (J, t + It should be possible to 
give an exact quantum-mechanical interpretation of such a reparametrization by coupling the discrete derivative 
Vjk = £j — £fc to the energy current operator fjk = —i[Hj, H^], where Hj, are local terms in the Hamiltonian. 
Such coupling is unambiguously defined only in the first order in rjjk , but that does not affect the thermodynamic 
limit. 

18 A rigorous and completely general upper bound for the group velocity was obtained by Lieb and Robin- 
son |77j . It is consistent with the rough estimate v < rJ, where J is the inter-spin interaction strength and r 
is its range. 
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where Of, is the topological spin of b. In the non-conformal setting, the role of the characters 
Xb{w) is played by linear combinations of the partition functions: 

Z b = J2 S baZa, (170) 

a 

where s a b are the entries of a topological S-matrix defined in terms of anyonic braiding, see 
Eq. (12231) on page [SU (We distinguish it form the modular matrix S, though they actually 
coincide when the latter is defined.) A calculation based on the insertion of a space-like Wilson 
loop Wb(£) shows that 

Z' b « e- 2mc ~ /2 %Z b . (171) 

This equation is indeed similar to (11691) with one important difference: unlike Xa{w) and 
Xa{— l/ w )> the numbers Z a and Z a are not obtained by evaluating the same function. An- 
other problem is that the Dehn twist can only be performed along a time circle. It is therefore 
not clear whether the full modular invariance can be established without conformal symmetry, 
or at least a 90°-rotation symmetry of the space-time manifold. 



Appendix E: Algebraic theory of anyons 

This appendix is an attempt to present an existing but difficult and somewhat obscure the- 
ory in an accessible form, especially for the reader without extensive field theory knowledge. 
(But we do assume some mathematical background and the willingness to follow abstract ar- 
guments.) The presently available resources may be divided into three categories: the original 
field-theoretic papers where the relevant mathematical structure was discovered [4j [5], local 
field theory theory papers [HI [7] (see Haag's book [7SJ for a general reference on this subject), 
and purely mathematical expositions [HSJ [79], |80j. Unfortunately, all these texts have different 
focus and/or are difficult to read. A nice elementary introduction can be found in Preskill's 
lecture notes on quantum computation [ST], but more detail is needed for our purposes. 

The theory described below is applicable to two-dimensional many-body systems with short- 
range interactions and an energy gap. Using some intuition about local excitations in such 
systems, we characterize the properties of anyons by a set of axioms and derive some corollaries. 
This approach is somewhat similar to what has been done rigorously in local field theory. 
However, we avoid many difficulties by keeping the discussion at the physical level until all 
essential properties are cast into a finite algebraic form, at which point we enter the realm of 
mathematics. 

More specifically, we claim that finite-energy excitations are classified by supers election sec- 
tors: each sector consists of states that can be transformed one to another by local operations. 
(We do not care about rigorous definition here.) It is assumed that the full classification can 
be established by considering a small neighborhood of an arbitrary point (of size compared to 
the correlation length). Another assumption is that local excitations can be moved from one 
place to another by applying an operator along an arbitrary path. Note that the quasiparticle 
transport need not be adiabatic, nor we require translational symmetry. Topological proper- 
ties, such as braiding rules, have invariant meaning independent of possible disorder, cf. Fig. [10] 
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Figure 16: Some "obvious" identities. 



on page EU Using some additional arguments, we will arrive at a theory that may be called 
unitary braided fusion category (UBFC) (= unitary ribbon category) in abstract language o 

With this approach, we get a weaker set of axioms compared to an analogous algebraic 
structure in conformal field theory. Therefore some "obvious" properties require a proof. (A 
similar situation occurs in local field theory, but our exposition is different in some details and 
hopefully simpler.) For nonrelativistic quantum- mechanical models, the space and time are 
not equivalent even if we use the imaginary time formalism. Thus the notion of an antiparticle 
being a particle propagating back in time must be applied with caution. The identity in Fig. [T6a 
illustrates a standard use of this notion. (Shown in the figure are space-time diagrams; time 
goes up.) On the left-hand side of the identity a particle a annihilates with its antiparticle a that 
was created as part of an da pair. If space and time were related by a rotational symmetry, then 
the zigzag on the particle worldline could simply be removed. In our formalism, this equation is 
satisfied by a suitable choice of normalization factors for the creation and annihilation operators 
of particle-antiparticle pairs. However, the relation in Fig. [T6b is a nontrivial theorem; it follows 
from the positivity of Hermitian inner product, see Section TE.21 (In the categorical language, 
this result reads: "Any unitary fusion category admits a pivotal structure".) 

Having formulated fusion and braiding axioms, we will touch upon another question: What 
UBFCs can be realized by a Hamiltonian with given symmetry properties? The simplest case 
is where there is no built-in symmetry or associated conservation laws. That is, we assume that 
the system consists of spins (or other bosonic degrees of freedom) and the Hamiltonian is in 
generic position. (Fermions are excluded because their number is always conserved modulo 2.) 
Thus, any effective conservation law for quasiparticles must have topological nature. In Sec. IE. 51 
we formulate this condition as a certain nondegeneracy property of the braiding rules. Theories 
satisfying this axiom are called unitary modular categories (UMC)ffil 

Unfortunately, UMC does not capture one important robust characteristic of the physical 
system, namely the chiral central charge c_ . However, there is a beautiful relation [3, [5T] that 
fixes c_ mod 8: 



19 Note for experts (the general reader shouldn't worry about this): a fusion category is a /c-linear semisim- 
ple rigid tensor category with finite-dimensional morphism spaces and finitely many simple objects such that 
End(l) = k, see 30 . We take k — C. A unitary fusion category is automatically pivotal, spherical, and 
nondegenerate. In the presence of braiding, these properties imply the existence of a ribbon twist. 

20 The term "modular" refers to the possibility to define a Hilbert space associated with a torus, on which the 
modular group acts. More generally, one can construct a so-called modular functor and a topological quantum 
field theory [67l [68l ES] (but we do not cover these topics) . 




where 




(172) 



a 
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d a is the quantum dimension of the superselection sector a (see Sec. |E.2p and 8 a is the topological 
spin (see Sec. IE.3j) . The left-hand side of Eq. ( I172p is denoted by 0; its algebraic meaning will 



be explained in Sec. IE.51 Roughly, 6 = S^^-TSTST, where S is a so-called topological S-matrix 
and T is the diagonal matrix with entries 8 a . We will see that is a root of unity, hence 
Eq. (11721) implies that c_ is rational. 

The standard proof of Eq. (11721) is based on the assumption that the system admits con- 
formally invariant edge modes under suitable boundary conditions. Then one can employ the 
modular invariance of the partition function on a space-time torus, which is established in the 
CFT framework [62] . It is not known whether the conformality hypothesis is true in general. 
One may, however, hope to find a proof that would not depend on conformal symmetry. In 
Sec. ID. 21 we made partial progress toward this goal. Specifically, we showed that the matrix 
T = e~ 2mc ~ l 2A T corresponds to the Dehn twist along a time circle. While S can also be re- 
garded as a modular transformation in a suitable mathematical theory, it is not clear whether 
that theory is applicable in the same context as the expression for T . 

In the remaining part of the appendix we do not discuss the chiral central charge. Our 
primary goal is to define a unitary braided fusion category (UBFC) and its special case, unitary 
modular category (UMC). The definition is rather long, so we break it into several parts. Each 
part contains some physical motivation, formal axioms in terms of basic data and equations, 
as well as some corollaries. Note that our formulation of the axioms does not involve the 
mathematical notion of category, so we call the corresponding objects "theories" (e.g., "fusion 
theory" instead of "fusion category"). Only later do we introduce the powerful but heavy 
language of categories and functors (also known as "abstract nonsense"). It is surely overkill 
for most of the problems considered in this paper, but the concept of tensor functor may 
actually be useful in the study of phase transitions in anyonic systems, e.g., Bose-condensation 
of spinless particle and weak symmetry breaking. (The last topic is considered in Appendix IFl) 

E.l Fusion theory 

In this section we consider anyons on a line (though the physical system is two-dimensional). 
The motivation for the definitions has been provided in Sec. 18.41 (esp. in 18.4.21) . Each particle is 
characterized by a superselection sector label (describing its "anyonic charge"). The position 
and local degrees of freedom may be ignored, but the order of particles is important. The main 
element of the the fusion theory is the space V c ab — the state space of particles a and b restricted 
to have total anyonic charge c. More exactly, vectors in this space correspond to different ways 
of splitting c into a and b, or equivalence classes of local operators that effect the splitting. 

E.l.l Splitting and fusion operators. While V^ b may be called "splitting space", the 
corresponding fusion space is denoted by V£ b . If ip e Vf' is a splitting operator, then e V^ b 
is a fusion operator: 

a b c 
splitting: , fusion: , (173) 



where ip and ^ label the corresponding vertices. One may also use Dirac's notation, i.e., write 
instead of ip and instead of ip\ By definition, the Hermitian inner product corresponds 
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to the operator multiplication. Specifically, if £,77 G V^ b } then 

c c 

i.e., r/t£ = (r/|£)id c . (174) 

c c 

Similarly, one can define splitting and fusion of multiple particles. The most general process 
is one transforming particles b%, . . . ,b m into a 1; . . . , a&. It can be performed in two steps: first, 
fusing all the particles into one, and then splitting it as required. Therefore 

C:C = 0C-" flt ®^.., m - (175) 

c 

For example, the identity operator acting on particles a and b can be decomposed as follows: 




(v\0 



a b a b 

v 

= EE % . ( 176 ) 

a b a b 

where j stands for both ipj and i/jj — basis vectors in Vf' and V£ b , respectively. 

The description of splitting and fusion will be complete when we understand relations be- 
tween different spaces l/ ai "'° fc . That is the topic of the following subsection. 

E.1.2 Basic data and equations. The basic data include a set of superselection sectors, 
fusion rules, and associativity relations. The latter satisfy so-called pentagon and triangle 
equations. An additional condition (requiring the nondegeneracy of certain operators) will be 
formulated in Sec. IE. 21 

Superselection sectors (also called "particle types", or "labels"): They form a finite set M. 

Fusion rules: For any combination of a, b, c G M there is a fixed finite-dimensional Hilbert 
space V^ b . The numbers N^ b = dim V£ b = dimV^ 6 are called fusion multiplicities. The choice 
of a special element 1 G M (the vacuum sector) is also considered part of the fusion rules. 

Associativity relations: For each a,b,c,d G M there is a canonical unitary isomorphism 
between two Hilbert spaces: 



a b c a b c 

K bc ■■ V e ab (8 V: c - V a J ® Vf. (177) 

u u 

Note that both spaces in question are simply different representation of the same physical space 
yabc — faofc Q f the particles a, b, and c restricted to have total anyonic charge u. Therefore one 
may regard the isomorphism F£ bc as equality. This view is formally justified by the pentagon 
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equation and MacLane's coherence theorem, which will be formulated later. At a cruder, 
combinatorial level associativity may be expressed as follows: J2 e ^tb-^ec = ^af-^L- 

From the physical perspective, splitting off (or fusing with) the vacuum sector is trivial. In 
any splitting diagram all lines labeled by 1 may be simply erased. To guarantee the consistency 
of this procedure, we require that the spaces V® 1 and V^ a are not only one-dimensional, but 
canonically isomorphic to C. In other words, these spaces are endowed with fixed unit vectors 
|a a ) G V® 1 and \(3 a ) G V} a . The canonical isomorphisms are given by the formulas 

a a :C^V a al : z^z\a a ), (3 a : C -> V} a : z»z\p a ), (178) 

where z6C These isomorphisms must satisfy so-called triangle equations (see below). 

Pentagon equation: The graphic representation of this axiom is shown in Fig. [T7a . Its 
exact meaning is this: for any x, y, z,w,u G M the following diagram commutes: 




(179) 



Arrow labels are abbreviated: for example, the arrow on the left (labeled by F^ yz ) actually 
designates the map ^ F* yz £g> id v v™, where id v v™ is the unit operator on V£ w . 

According to MacLane's coherence theorem [82J, any sequence of F-moves between two 
given trees results in the same isomorphism between the corresponding spaces: the equality of 
all such isomorphisms follows from the pentagon equation. (Actually, the theorem is a bit more 
general and also includes a- and /3-moves, but we ignore them for the moment.) This result 
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a) b) 

Figure 18: The graph of F-moves between binary planar trees with n = 5 leaves (a) and the 
procedure used to assign coordinates x\, . . . , x n -\ to each tree (b). 



may be regarded as a combinatorial statement. Indeed, let us consider the graph T n whose 
vertices are binary planar trees with n leaves and whose edges correspond to F-moves. (The 
graph is shown in Fig. fT8k.) Then the theorem says that any cycle in the graph can be filled 
with pentagons and quadrilaterals. The latter correspond to the obvious fact that independent 
F-moves commute; such moves can be nested as in Fig. [TTb or occur in disjoint branches. 
MacLane actually shows that any two directed paths between the same pair of vertices can be 
transformed one to the other using the pentagon and quadrilateral equations; the proof is by 
induction on a suitable parameter. 

Let us also mention a beautiful geometric proof of the coherence theorem: the graph T n 
together with the set of 2-cells (i.e., the pentagons and quadrilaterals) can be realized by the 
vertices, edges, and 2-faces of some convex polytope K n of dimension n — 2. This polytope 
appeared in the work of Stasheff [83J and now bears his name; it is also called associahedron. 
More exactly, this term refers to a certain combinatorial type of a polytope while its exact 
shape may vary. The faces of K n in all dimensions are associated with general (not necessarily 
binary) planar trees. In particular, the edges correspond to trees with one degree-3 vertex, the 
quadrilaterals correspond to trees with two degree-3 vertices, and the pentagons correspond to 
trees with one degree-4 vertex. Described in Fig. [T8b is a concrete geometric realization of K n 
(cf. Ref. [M]). It is fairly easy to show that the convex hull of the points (xi, . . . , £ n -i) £ M n_1 
obtained this way is characterized by one equation and n ( n-1 ) — \ inequalities: 

n-1 l-l 

y*]xj = r] ln , > Vki (l^k<l^n), where r] U = ^ Vuv (180) 

j=l j=k k^u<v^l 



70 



X w 



X w 



X w 




a) b) c) 

Figure 19: The fundamental triangle equation (a) and its corollaries (b and c). 



Triangle equations: The three equations in Fig. [T9] guarantee that adding or removing 
trivial lines commutes with F-moves. Let us write the first equation in a more conventional 
form: 




commutes, i.e. 



yxw y lw 



pxlw 



(\a x ) ® = ® \/3 w ) 



for any 
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Lemma E.l. a\ = [3\. 

Proof. Let x = u and w — 1. Then first and the third triangle equations coincide, except that 
the right arrow is (3\ in Fig. [T9~a and cti in Fig. [T9~b . But the right arrow is the composition of 
the other two, therefore it is the same in both cases. □ 

Lemma E.2 (cf. (HDj, Lemma XI.2.2). The first triangle equation (together with the pentagon 
and quadrilaterals) implies the second and the third. 

The statement of the lemma refers to formal tree calculus. Moves between trees are bidirectional 
(therefore we may represent them by lines rather then arrows). The quadrilaterals express the 
commutativity of disjoint F-, a-, and /3-moves. 

Proof. Let us prove the equation in Fig. [19b. 

First, we join an additional trivial branch to 
the trunk of each tree and show that the result- 
ing equation is equivalent to the original one. In 
the diagram on the right, the old and the new 
equation constitute the top and the bottom of a 
triangular prism, respectively. The sides of the 
prism are commutative quadrilateral, therefore 
the top triangle commutes if and only if the bot- 
tom one does. 
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We now consider the bottom of the prism, 
which may be identified with triangle 4-7-5 in 
the new diagram. In this part of the argument 
it is not important that the leftmost branch of 
each tree is trivial, so we represent it by a solid 
line. The outline of the diagram corresponds 
to the pentagon equation. Triangles 3-1-6 and 
4-2-5 are instances of the equation in Fig. [T§a . 
which is true by hypothesis. The quadrilater- 
als 1-2-4-3 and 6-3-4-7 also commute. Thus the 
required triangle 4-7-5 commutes as well. 

The equation in Fig. fl~9b is proved analogously 

E.1.3 Examples. Several concrete theories have been presented in the main text, see tables 
on pages [30], IUJ |42j Let us now describe some general constructions leading to infinite series 
of examples. 

1. Let the label set M be a finite group and the fusion rules correspond to the group 
multiplication. That is, the space V^ ,y is one-dimensional (with a basis vector ijj^ y ) 
if z = xy, otherwise V^' y = 0. The associativity relations are trivial. Note that if the 
group M is non-Abelian, this fusion theory does not admit braiding. 

2. The associativity constraints in the above example can be deformed as follows: 





f(x,y,z 




(182) 



.vv- 



where the vertices of the trees correspond to basis vectors and f(x,y,z) are some phase 
factors. In this case the pentagon equation reads: 

f(xy,z,w) f(x,y,zw) = f(x,y,z)f(x,yz,w)f(y,z,w). (183) 

This is a cocycle equation, i.e., / is a 3-cocycle on M with values in U(l). One can show 
that the triangle equations do not put any additional restriction on / but rather define 
a and (5: 

a x = C f(x, 1, 1)~ A« = C /(I) I; w )^w W j (C is an arbitrary phase) . (184) 

The basis change ip% ,y — > u(x,y)^*y will result in a new cocycle g such that 

f(x,y,zY l g(x,y,z) = u(x,y)u(x,yz)~ 1 u(xy,z)u{y,zy 1 . (185) 

The right-hand side of this equation is the coboundary of the 2-cochain u, i.e., g and / 
belong to the same cohomology class /. Thus, the associativity relations are classified by 
# 3 (M,U(1)). 
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For a concrete example, let M = Z 2 = It is known that if 3 (Z 2 , U(l)) = Z 2 , i.e., 

Eq. (11831) has one trivial and one nontrivial solution. The nontrivial solution is given by 
f(a,a,a) = — 1, the other seven values being equal to 1. This particular fusion theory 
admits two braidings: the counterclockwise exchange of two a-particles is characterized 
by either +i or —i. Such particles are called semions. Recall that Case 2 in the table on 
page HH represents two independent species of semions. 

3. Let M be the set of irreducible representations of a finite group G and let Vf' = 
HomqG] (c, a ® b) . In other words, elements of V c ab are intertwiners between c and a <S> b, 
i.e., linear maps that commute with the group action. The associativity relations are 
given by 6j'-symbols. This theory admits trivial braiding. 

4. Examples 3 and 1 (with M — G) can be combined in what is called representation theory 
of Drinf eld's quantum double [85, 80J. Models with anyons described by this theory 
were proposed in [86] (gauge-symmetric Lagrangian) and in [24] (lattice Hamiltonian not 
constrained by an external symmetry) . This construction can be deformed by an arbitrary 
cohomology class fe H 3 (G, U(l)), see Refs. [83188]. 

5. A very interesting set of fusion theories is based on a Temperley-Lieb category, see e.g. 
Chapter XII in book [HE]- These theories are also known as "representations of quantum 
SU(2)". 

E.1.4 Calculations with planar graphs. The rules we have described allow us to work 
not only with trees but with arbitrary oriented planar graphs. In particular, loops can be 
removed using Eq. (11741) . If we encounter a subgraph like shown in Eq. (j!74p but with different 
labels at the top and at the bottom (say, d and c"), then its value is zero. Let us illustrate 
that by a concrete example, using the fusion theory from Table [Hon page 1301 for h = +1: 



a a £ ooe a a £ a a ^ 




11 11 



In general, we consider planar graphs that satisfy the following conditions. 

1. Vertices are associated with fusion/splitting spaces V? l '"' l ' Ua ', the indices are represented 
by prongs. There may also be external labels (terminals) at the bottom and the top of 
the graph. Bottom terminals and upper indices are called sources, lower indices and top 
terminals are called targets. 

2. Each edge connects a source to an identically labeled target. (In pictures we use a single 
label for the whole edge.) 

3. The edges are transversal to the horizontal direction and oriented upwards. (We do not 
show the orientation in drawing because it is obvious.) 
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To each vertex we also assign an element of the corresponding space V£ u "i' ". Using asso- 
ciativity relations and Eqs. (1174[) . (11761) . we can compute the value of the graph, X G V^ 1 '"'^, 
where ax, . . . , a* and bx, ■ ■ ■ ,b m are the external labels. Note that we may freely add or remove 
edges labeled with 1, thanks to the triangle equations. 

A general calculation strategy is based on the idea that X can be represented by a set of 
maps X c : V^ 1 ' '"' bm — > V c ai ' "' afc for each label c; this representation is closely related to the 
decomposition of identity (11761) . Thus we apply Eq. (1 1 76 [) first, and then start shrinking and 
removing loops as illustrated above. For example: 



a a 

e. 






a a a o 



V 

1 1 

A 



a a a a 



We will see that for theories with particle-antiparticle duality, condition [3] can be dispensed 
with. The formalism will be revised so that edges can bend, changing their direction from 
upward to downward and back, and vertices can rotate by 360° as in Fig. [T()b . 

E.2 Particle-antiparticle duality 

In this section we complement the fusion theory by the following condition. 

Duality axiom. For each label a there is some label a and vectors |£) G V® a , \q) G V^ a 
such that 

<a a ®r / |Fr a |e®/3a>^0. 

Note that the matrix element in question corresponds to a physical process in which the a- 
particle from an aa-pair annihilates with a different copy of a: 




[a a ®7l\*? a \Z®Pa) 



;i86) 



(The dotted lines and the associated labels, i.e., 1, a a , j3 a may be ignored.) 

Let (a a <S> v\Fa aa \£, ® Pa) = u, assuming that |£) and \rj) are unit vectors. We will see that 
the spaces V^ a and V^ a are one- dimensional, therefore \u\ is uniquely defined. The number 
d a = | | " 1 , called quantum dimension, plays an important role in the theory. In particular, 
d a — 1 for Abelian particles (i.e., such that a x a — 1), otherwise d a > 1. If a = a, then we may 
set |£) = \rj) so that u itself has an invariant meaning. Specifically, 




d a 



187) 
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where x a is a uniquely defined phase, which is actually equal to ±1 (see below). This number 
is called Frobenius-Schur indicator; in the present context it was introduced by Fredenhagen, 
Rehren, and Schroer [89] . 

E.2.1 Physical motivation for the duality axiom. The existence of antiparticles follows 
from a locality principle mentioned at the beginning of this appendix: a quasiparticle can be 
moved from one place to another by applying an operator that acts on spins along a path 
connecting the given points. (We are still considering anyons on a line, therefore the path is 
unique.) The action on different spins can be performed at once or in any particular order. For 
example, to move a particle from point 3 to point 1 on its left, it may be natural to move it 
first to some middle point 2, and then to the final destination. But it is also possible to act on 
the interval [1,2] before [2,3]. The intermediate state must be a particle-antiparticle pair. 

For a slightly more rigorous argument, let X be an operator that moves a particle of type a 
from 3 to 1, acting on spins along the interval. We can represent it as follows: 



X 



1 



/] Y k Z k , 



(188) 



where Y k and Z k act on disjoint sets of spins in some regions around [1, 2] and [2, 3], respectively. 
Let n?! (s = 1,2,3) be the projector onto states having particle x at point s; such operators 
can be realized locally (to act in the gray circles in the above picture) if we only allow states 
that are not excited away from the given three points. Finally, we define new versions of X, 
Y k , and Z k : 

x' = n^n^n^ xn^ripjii^ = y] z' k Y k \ 

k (189) 



Z' 



n^np] z k n" 3 ], 



Y! 
1 k 



nfi] Y k n^rip]. 



The operator X' still moves a from 2 to 3, but Y k ' and Z' k overlap geometrically and therefore 
no longer commute. It is clear that each operator Y k creates an a-particle at point 1 and some 
particle at point 2, whereas Z' k removes the second particle as well as the original particle at 
point 3. Thus each product Z' k Y k ' effects the process shown in Eq. f)186p . up to an overall factor. 
This factor is nonzero for at least one value of k. 



E.2.2 Some properties and normalization conventions. 

Lemma E.3. For each label a the corresponding label a is unique. Moreover, N^ b 

dcf 



5, 



ba ■ 



Proof. Recall that N^, 



dimV^ = dimF^ 6 . We first show that this number is equal to one 
if b = a and zero otherwise. Let £, 77 be as in the duality axiom, and let u be the corresponding 
matrix element. For an arbitrary element ip G V® we have: 




some 
operator 



7id s for some 7 G C if b = a 
if b ^ a. 
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Therefore 




u X 7 £ if b 



1p = U as a/ a\ b, = 

by duality axiom \ / \ / duo to the above Otherwise. 



To prove that N^ a = <5 te , we do a similar calculation with tp' £ V[ a , putting (V'O* on the left 
of£. □ 

Corollary E.4. 1 = 1 and a = a. 

Let us now set new normalization conventions. In the definition of d a we assumed that |£) 
and \r]) are unit vectors, i.e., = r^i] = idi. However, it is more convenient to multiply 
£ and 77 by \fd~ a (and some suitable phase factors) so that the zigzag on the left-hand side 
of equation ( 11861) could be simply removed. The so normalized operators for the creation of 
particle-antiparticle pairs will be represented as smooth "cups" with a triangle mark at the 
bottom. The adjoint operators are "caps" obtained by flipping the pictures about a horizontal 
line@ 

a = (uY. a = fuY (190) 



In all four cases, the triangle points from a to a as we follow the line. The four operators are 
defined up to a single phase and satisfy the following relations: 

d a = a(^)«, (191) 



a, a I a 



id a , S ''^\ 1> ^ = a/A> \^J = ida ' a K^f^ = ( 192 ' ) 

Eq. f)19ip and the first equality in (11921) are true by definition. The second equality can be proved 
analogously to Lemma lE"3l the last two equalities are obtained by passing to adjoint operators. 
Roughly, the normalization conditions may be described as follows: oppositely oriented triangles 
cancel each other. 

So far the label a was fixed. Repeating the same procedure for a, we get a new set of basis 
vectors in V"", V^ a , V a \, V^ a , which must be related to the old ones: 

a a a 5 a a a a 

V = Xa V ' V = H ~ a V ' 

(193) 



a a 



Combining these relations with Eqs. ( 11911) . ( 11921) . we get: 

\Xa\ 2 = d- a ld a = \x a \~ 2 , X a (x?) _1 = 1. 



21 This notation is not standard, but it is convenient for calculations. In the notation used in [551 [751 150]. 
cups and caps are not decorated but each particle has a fictitious degree of freedom: it is considered as either 
a going up or a going down. This prevents noninvariant phases from appearing in formulas. 
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Thus, 

4 = 4, K- a = K* a , \x a \ = l. (194) 

In general, the number x a depends on the arbitrary phases in the definition of cups and caps. 
However, if a = a, then x a is defined uniquely and coincides with the Frobenius-Schur indica- 
tor @ In this case, x a = ±1. 

Now we are in a position to relax the requirement that edges are transversal to the horizontal 
direction. 

Definition-Proposition E.5. A line is an alternating sequence of 2n cups and caps without 
triangle marks. It may be open or closed. Such an object has a canonical normalization given 
by n triangles pointing forward and n triangles pointing backward (relative to a chosen path 
direction). All such decorations are equivalent. 

Let us also renormalize the inner product on splitting spaces: 



where (,J|G^. (195) 



(Triangles may be added in one of the two consistent ways; one may also flip the c-loop to the 
left.) This renormalization does not violate the unitarity of associativity relations. The inner 
product on fusion spaces is defined through the antilinear isomorphism V c afe — > V^ b : ip i— ► ip* . 
When computing either type of inner product, we stack two trees with branches touching each 
other and connect the roots by a loop. Stacking the trees root to root might produce a different 
result, but it turns out to be the same since we can actually rotate vertices by 360°, see below. 
Finally, we rewrite Eq. fl 1 76 [) in a form that is consistent with the new rules: 



a b 



a b 



EE 



c 3 



d, 



c 



c 



44 

a b a b 

ab r- \rab 



where 



'iio = ^ ( i9e ) 



(the elements ip^ G form a complete basis). 
E.2.3 Raising and lowering of indices. Let us define the following linear maps: 



a b 

c 

They are obviously invertible, hence N^ b = N% c = N b5 = = N^ a = N%. Unfortunately, this 
does not save us from the need to to distinguish between lower and upper indices because the 

22 In category theory, x a is not a scalar but rather a morphism from a to a, which are regarded as different 
(albeit isomorphic) objects. Thus the arbitrariness goes away, but the Frobenius-Schur indicator has to be 
defined in a more complicated fashion. 
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Notation used in this paper 


Isotopy-invariant calculus 


Vertices 


ah a 3 ^ a a 
a a a a 


y ^ a -»-<^«- a = a -*T>-«- 5 


Lines 


~ a (x^> a O' a 


"0- a O 


Cancellation rules 


a^aja^ = fl^s|a^ = id a 





Table 5: Notation summary. 



simultaneous raising and lowering of indices on the two ends of a line results in the factor x a . 
For example: 





where r\ G V^ x , £ G V c 



o6 



(198) 



However, as the following theorem shows, we do not have to tie the types of indices to the 
orientation of the corresponding vertex. This way, one obtains an isotopy-invariant calculus for 
planar graphs, see Table [51 We will continue using the previous notation, though. 

Theorem E.6. The maps A ab , B^ b are unitary with respect to the inner product ((-|-)), and the 
following diagram (in which the arrows may be traversed in both directions) commutes: 




(199) 



Note that the commutativity of this diagram is equivalent the "pivotal identity" in Fig [16b. 
The proof of the theorem is preceded by two lemmas. 

Lemma E.7. The composition of any two adjacent arrows in Eq. $199\) is unitary. 

Proof. Due to symmetry, it suffices to consider just one particular case, e.g., the arrows con- 
nected at the upper left corner. Let X = (Bf)~ l Af : V c ab -> V^. For arbitrary elements 
£ eV c ab , 7] G V£ 5 we have: 




X^T) 
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Therefore X 1 = X* . 

Now consider the two paths from top to bottom of the hexagon (j!99p : 



□ 



U h U r : V c ab 




U r : V 



b a 




. (200) 



b a 



From the physical point of view, these operators (which are actually equal) correspond to the 
CPT symmetry. 



Lemma E.8. U, 



i 

rab 



u. 



-1 



Proof. Let £ e V™ and rj e V{- Then 








'\c b( S/^K 










\ b 





Vd a d b d c ((U^r}\0). 



□ 

Proof of Theorem IE. 61 Let us consider the path around the diagram (I199P in the counter- 
clockwise direction, W = U~ l Ui. This operator is a composition of six arrows, therefore it is 
unitary (due to Lemma IE. 7ft . On the other hand, Lemma [E.8I implies that W = Ujlli is Her- 
mitian and positive semidefmite. It follows that W is the identity operator, i.e., the diagram 
commutes. 

Thus, Ui = U r ; let us denote this operator simply by U. It is unitary by Lemma [E. 81 Any 
arrow in the diagram is a composition of U (or U^ 1 ) and some number of arrow pairs. Therefore 
all the arrows are unitary. □ 



E.2.4 Quantum dimension and fusion rules. This is a key identity: 

d a d b = ^N c ah d c . 

c 

To prove it, we use Eq. (I196p . the pivotal property, and Eq. (11951) : 

b 



(201) 



d a db 




a) )b 




} ab 
c,3' 



{ififltY rotating ip? b 



c ' 3 c,j 
counterclockwise 
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Eq. (12011) implies that the matrix N(a) = (N£ b : 5, c G M) has eigenvector v = (d c : c G M), 
and the corresponding eigenvalue is d a . Note that all the entries of v are positive. According 
to the Perron- Frobenius theorem (about matrices with nonnegative entries), all eigenvalues of 
N(a) satisfy the inequality 

\Xj{a)\^d a . (202) 

Thus, the quantum dimension d a is the largest eigenvalue of the matrix N(a). 

The proof of Eq. (12011) also motivates a definition of quantum trace. It is a number assigned 
to any element X G V^ 1 "^™. Such an element acts as an operator in the space y^- a ^ for each 
c; we denote this action by X c . The quantum trace is defined as follows: 




(203) 



E.3 Braiding and topological spin 

Braiding is an additional piece of algebraic structure characterizing anyons. It is defined by a 
set of elements R a b G V^ b a that represent transposition of two particles: 



R 



ab 



b a 

Yv 

a b 



(Ra 



v-1 



OR, 



ab 



a b 

Y 

b a 



(204) 



They must satisfy the Yang-Baxter equation, i.e., a line can be moved over or under a crossing 
between two other lines: 

c b a c b a 

(205) 



a b c a b c 

More generally, a line can slide over or under an arbitrary vertex. It is sufficient to postulate this 
property for three-prong vertices representing splitting and fusion and for lines tilted left and 
right. Furthermore, splitting and fusion are related by the decomposition of the identity (1 1 9 6 j) . 
whereas changing the slope of the intersecting line is equivalent to replacing R with Br 1 . Thus 
the number of independent conditions is reduced to two: 



y z x 



y z x 



1> 



(206) 



The description of braiding in terms of basic data amounts to specifying the action of R a b 
on splitting spaces: 



b a 
< J 



R ab . yab 



ba 



(207) 
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a) b) 
Figure 20: The hexagon equations. 



Note that R^ b is a unitary map, therefore N£ b = N ba . 

To express Eq. (I206p in terms of R^ b , we join the two lines at the bottom of each graph and 
perform equivalence transformations. These include F-moves as well as R-moves — absorbing 
a line crossing by a vertex. Thus we obtain the diagrams in Fig. [201 the bottommost arrow in 
each of them combines an i?-move with the first or the second equation in question. We may 
now forget about the topological meaning of braiding and only keep track of the linear maps 
involved: 



0p v? ® V? 





® T v: r ®vr 




e v? ® v> 



r (id, r yz®R* r ) - swap 

v r 




(20£ 






0« vp ® v; 



r ® v r yz 



),.(id^ z (g)( J R™)- 1 )-swap 




(209) 



These commutative diagrams are known as hexagon equations. They actually look nicer in the 
tensor category formalism, see Eqs. (12741) and (I275p . 
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Note that braiding with label 1 is trivial: 



1 a 



1 a 


1 a 




\y 










i.e. 


a 


a 





B al a 



lu \ 



Pa- 



(210) 



(The proof is analogous to that of Lemma \E. 21 ) Due to this property, we need not worry about 
lines labeled by 1, e.g., ones that are implicitly attached to each cup and cap. 
To each label a we associate a complex number 8 a , called topological spin: 



0a ^ C/- 1 




= d^Tr R aa = d-^dc^R 

c 

It may also be characterized by any of the following relations: 



aa 
c m 



(211) 



a a a a 

Note the following properties of the topological spin: 

\0a\ =1, @a = Qa- 

Indeed, 



(212) 



(213) 



d a 0a0* a 




(I, 



d a a 





d a 9 a . 



This is yet another expression for the topological spin: 



0* 



(214) 



(To prove it, we need to put a left-oriented cap on both sides of the equation and rotate the 
resultant figure-eight on the left-hand side by 90°.) In particular, if a — a, then we have the 
following representation for the invariant scalar R\ a : 

R? = 0> a . (215) 

Remark E.9. Eq. (121 5[) provides a simple physical interpretation for the Frobenius-Schur 
indicator x a . Suppose that the system is rotationally invariant, so that not only the topological 
spin but also the usual spin s a has physical meaning. Of course, s a may take different values 
subject to the constraint e 2mSa = 6 a . Let us consider a pair of identical particles with trivial 
total topological charge. What values does the total angular momentum of this composite 
object take? This question may be answered as follows. Assuming that both particles are in 
the same spin state s a and barring additional "isospin" degrees of freedom, the 180° rotation 
is characterized by the phase factor e 47rs a e «7rs a ^aa _ x ^ Thus, the total angular momentum is 
even if x a = 1 and odd if x a = — 1. 
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It is interesting that the effect of moving one particle around another is fully characterized 
by the topological spin: 

Rh c R c = WT id v^- ( 216 ) 



Proof. Consider the following element of the space V°£ : 




a 



— ^< — RbaRabOaOb- (217) 

PP 

a b 

Acting by it as an operator on the space V^, we obtain the required identity. □ 

Theorem E.10 (Vafa [90]). The topological spins are roots of unity. More specifically, 0™ = 1, 
where the integer n^O depends only on the fusion multiplicities. 

Proof. Let us use this fact from linear algebra: for any integer matrix X there is an integer 

matrix Y (the adjugate) such that YX = (detX)I. Therefore, if Y\ b 0* ab = 1 for all a, then 
£detx = 1 for all b 

To find multiplicative relations between topological spins, we consider determinants of R- 
and F-moves. Such determinants generally depend on the choice of basis in the source and the 
target space. So, let us fix a basis in each space V^ b and replace the linear maps in the hexagon 
equations (12081) and (I209j) by their determinants. Dividing the first equation by the second, we 
get: 

det [®{Rl v Rf) ® id Vir j det ^0 id vr ® (R x q z R z q x )^j = det ^0 id vr ® (R*JR r u 

Note that the determinants here are actually invariant and can be expressed in terms of topo- 
logical spins using Eq. H216[) : 

n [efj n [efj = n \jj r ) 

Y[ 9 f y Ni. + NLN* yv+ K,Ns v = (e x9yBz9v) N ^ , Nxyzv « y, (21s) 

p q 

where we have substituted v for u. 

For our purposes, it is sufficient to consider the case where x = y = a and z = v = a. Thus 
we obtain a system of equations described by a square matrix: 

II ^ = 1. Where X «» = ^a-a-a <U ~ 2J\&JVj, - iV^A^. (219) 

6 
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It remains to show that X is nondegenerate. For this purpose, we may divide the a-th row of 
X by N aa aa so that the matrix becomes 41 — Y; the eigenvalues of Y are bounded by a suitable 
multiplicative matrix norm: 

Y ab = (N^)- 1 {2N b a - a Nl- a + N b aa NL) , (eigenvalue,. (Y) \ ^ max £ \Y ab \ = 3. 

b 

Thus the matrix 41 — Y cannot have a zero eigenvalue. □ 



E.4 Verlinde algebra and topological ^-matrix 

In this section we investigate some properties of fusion, based on the axioms stated above. 
First, let as consider a fusion theory with particle- antiparticle duality, but without braiding. 
The Verlinde algebra is an associative *-algebra spanned by elements e a : a £ M which satisfy 
the following relations: 

e a e b = J2 N ab^ 4 = e s . (220) 

c 

The basis element e a £ Ver may be represented by the matrix N(a) = (N^ b : b,c £ M); this 
representation is faithful. Since N(a) = N(ay, the algebra operation f corresponds to taking 
the adjoint matrix. Therefore Ver is actually a finite-dimensional C*-algebra. Equation (12011) 
says that the linear map Ver — > C : e a i— > d a is a homomorphism. 

If braiding is defined, then N° b = N ba , therefore the Verlinde algebra is commutative. A 
finite-dimensional commutative C*-algebra splits into several copies of C. Therefore we have 
\M\ distinct homomorphisms from Ver to C, which we refer to as fusion characters: 

X, : Ver - C, A,(a)A,(6) = £ iV a c 6 A,(c). (221) 

c 

The vector Vj = (Aj(c) : c £ M) is a common eigenvector of the matrices N(a), the eigen- 
values being equal to Xj(a). Note that the matrices N(a) are normal, therefore eigenvectors 
corresponding to distinct eigenvalues must be orthogonal. Thus the following orthogonality 
condition holds: 

Y^^(a)X k (a)* =0 tij^k. (222) 

a 

We will now construct a map from the set of labels to the set of fusion characters. Specif- 
ically, we are to show that some of the characters are given (up to a constant factor) by the 
columns of the topological S -matrix S = (s ab : a, b £ M), where 




W^A = I^4Tr(«) = ^£A^< 



The normalization factor 1/D with T> = v/5^ tt ^ s chosen so that S is unitary under certain 
conditions. 

Let us study properties of the numbers s ab . First, we observe these symmetries: 

Sab = S ba = S ab - = Sfa = = S* ba = S*g = S^. (224) 
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Indeed, up to normalization we have: 




Factors of the form s ax /si x often arise via this equation: 

X X 




(225) 



X X 



For example, let us repeat the calculation done in proof of Eq. (12011) . but now with an additional 
line passing through the loops: 



0< 






The first and the last expression are equal to Sax ^ hx id- an d N^ b ^-id XJ respectively. Thus, 

S lx „ S la: 



SaxSbx 
Six 



c c 

ab°cx- 



(226) 



Comparing this equation with Eq. (12211) . we conclude that 



— = \j(a) 

Six 



for some j = j(x). 



(227) 



In other words, for each x the map a ^ s ax /si x is a fusion character. 

In the next section we will see that if braiding is sufficiently nontrivial, then the S'-matrix 
is unitary. In this case, the map from labels to fusion characters is one-to-one, and Eq. ( 12261) 
may be cast into a form known as Verlinde formula^ 



SaxSbxSc 
Six 



(228) 



23 Eqs. (|226p . (|228|) were first obtained by Pasquier for statistical models and conformal theories associated 
with Dynkin diagrams [100] . A general derivation in the CFT framework is due to Verlinde [55] . 
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E.5 Braiding nondegeneracy = modularity 

Let us recall the definition of a superselection sector: it is a class of states that can be trans- 
formed one to another by local operators. From the physical perspective, the operators we 
use must respect any unbroken symmetry present in the Hamiltonian, e.g., the U(l) symmetry 
associated with the conservation of electric charge. Another example is the number of fermions 
modulo 2, if that number is not conserved as an integer. (Here we speak about actual fermions 
forming the system rather than effective Majorana modes obtained by a nonlocal transforma- 
tion.) However, in the model studied in this paper the superselection sectors are stable with 
respect to all local operators. More generally, we may consider an arbitrary system built of 
spins or other bosonic degrees of freedom, and ask for properties that are stable to a generic 
perturbation. In this case, the superselection sectors have purely topological nature: some sig- 
nature of a nontrivial excitation 1 ^ 1 will be preserved even if we cut out a piece of material 
containing the quasiparticle. It is quite reasonable to assume that the presence of such an 
excitation can actually be detected by an Aharonov-Bohm measurement, i.e., by moving a test 
particle a around x. Thus we arrive at the condition that the braiding is nondegenerate, which 
must be true for any anyonic system not complicated by external symmetries. 

Definition E.ll. Braiding is said to be nondegenerate if for each label a; ^ 1 there is some 
label a such that the operator R ax R xa is not identity. 

Theorem E.12. In a theory with nondegenerate braiding, the following operator S z acting in 
the space C z = ® b V£ z is unitary: 

b z 

s z y> = ^E d «rr^ ' where ^ ev h z - ( 229 ) 



Note that the standard ^-matrix corresponds to z — 1. The operator S\ differs from S z by the 
orientation of crossings. (This is easy to show by considering a matrix element of S z .) The 
converse of Theorem IE. 121 is also true and can be proved easily: if the standard S"-matrix is 
unitary, then the braiding is nondegenerate. 

Lemma E.13. The operator R ax R xa is trivial for all labels a if and only if x is mapped to the 
trivial fusion character j(x) = i.e., if s ax /s\ x = d a for alia. 

Proof. The triviality of RaxRxa means that for all c, R a c x R x c a = id xa , i.e., «%- = 1 whenever 
N xa 7^ 0. We calculate s xa = s* ax using Eq. (12231) : 

Vs xa = = ^Kadc = d x d a . 

c x a c 

Thus, s ax /si x = d a . Conversely, if the operator R ax R xa is nontrivial, then 

v Res xa = J2 N xa Re \Te) dc < ZX c » dc = dxda > 

c \ x a / c 

therefore s ax /s lx ^ d a . □ 
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The proof of Theorem IE. 121 is based on the following equation, which is useful on its own 
right: 



V 2 



i(«0d'(i) 



I 

X 



(230) 



Indeed, the condition j(x) = j(l) means that the fusion character \j( x ) : a i— > s ax /s\ x is trivial, 
i.e., s ax / s\ x = da- On the other hand, if \jt x ) is nontrivial, then it is orthogonal to the trivial 
character, i.e., ^2 a d a s ax = 0. The above equation combines both cases. 



Proof of Theorem IE. 121 By Lemma [E. 131 the nondegeneracy condition implies that j(x) = 
if and only if x = 1. Thus, the linear combination of loops in Eq. (12301) is a projector onto 
label 1. We apply it as follows: 



V 2 



x y 

i V 

r 3 X y 

dx 

x y 



(231) 



x z x z 




X X 



□ 

Theories with a unitary S'-matrix are often called modular due to their relation to the 
group of diffeomorphisms of the two-dimensional torus considered up to topological equiva- 
lence (namely, isotopy). Let us discuss this subject at a speculative level. The first thing 
to note is that a system of spins or electrons can in principle be put on a torus whether or 
not the braiding is nondegenerate. However, as indicated earlier, degeneracy generally occurs 
due to some external symmetry. If it is a gauge symmetry, then putting the system on the 
torus involves some choice. For example, in the case of U(l) symmetry, one may run an ar- 
bitrary magnetic flux through each basic circle of the torus. Those fluxes can be detected 
by Aharonov-Bohm measurements, but they cannot be changed by intrinsic operations, which 
include splitting, fusion, and moving particles around the torus. Mathematically, such oper- 
ations form a so-called skein algebra, which may be defined in terms of graphs on the torus 
up to equivalence transformations. (Recall that a planar graph without external lines can be 
transformed to a number, therefore the skein algebra of the plane is simply C.) The skein 
algebra is a finite-dimensional C*-algebra; in general it is a direct sum of several blocks, each 
block being isomorphic to the algebra of operators on some finite-dimensional Hilbert space. 
Physically, different blocks correspond to different values of external parameters such as the 
magnetic fluxes. Modular theories are special in that there is only one block. In other words, 
the torus is characterized by a single finite-dimensional space. 
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The Hilbert space of the torus is actually M. = C AI . A basis in this space may be associated 
with effective anyonic charge that is detected by an Aharonov-Bohm measurement along some 
circle. Different circles correspond to different bases. The S-matrix effects a transition between 
two bases. Other important matrices are C = (5ab '■ a, b G M) (which corresponds to a f 80°- 
degree rotation) and T = (6 a d~ a b : a, b e M) (which corresponds to a Dehn twist). We will 
justify this description by showing that S, C, and T obey the modular relations up to a phase 
factor 9: 

(ST) 3 = BC, S 2 = C, C 2 = I, (232) 

where 

Q = V- 1 J2d%. (233) 

a 

As a corollary, note that B is a root of unity. Indeed, 5* 4 = I, therefore B 4 ' M ' = (detT) 12 is a 
product of topological spins, which are roots of unity by Theorem IE. 101 

Let us actually study relations between more general operators acting in the space C z = 
®b (which corresponds to a punctured torus, with z being the anyonic charge of the punc- 
ture). The operator S z has been already defined, whereas 

b z b z b z b z 

b b b b 

Theorem E.14. The operators S z , C z , and T z satisfy the following modular relations: 

S\T Z S Z = T\S\Tl S z = SlC, Cl = 9* z . (235) 

Proof. The calculation of S\T Z S Z parallels the proof of Theorem IE. 121 but the nondegeneracy 
condition is not necessary. First, let us apply the S'-matrix to the row vector with entries d a 6 a : 

^2d a 6 a s aS = -^^d a N c ax ^d c = ^^d x d c ^d c = d x 9* x . 

a a,c x c x 

The result may be written as follows: 



X X 




(236) 



Then we replace the single line by two lines and simplify the result (cf. Eq. (12171) ): 



x y 




(237) 



N8 



Attaching a graph representing if) £ Vg Z to the first and the last expression, we get S\T z S z ifj 
and QT%SlT%if>, respectively. 

Now we show that S z = S\C Z . If ip G V b bz , then 




E.6 Gauge freedom and Ocneanu rigidity 

In the preceding sections, we defined a theory of anyons as a solution to a certain system of 
algebraic equations for F® bc , a a , j3 a , and R" b . Specifically, these are the pentagon equation, the 
triangle and hexagon equations as well as unitarity conditions. (The unitarity conditions are 
algebraic in the real and imaginary part of the corresponding matrices.) The set of solutions 
is some real algebraic variety. This description, however, does not take into account a gauge 
degree of freedom. Indeed, two solutions are physically equivalent if they are related to each 
other by a simultaneous basis change in the spaces V^ b . In the examples studied in this paper, 
quotienting over such transformations makes the set of solutions discrete. In other words, the 
algebraic structure is rigid, with the physical consequence that the properties of anyons are 
stable to small perturbations of the Hamiltonian. Is it true in general? The affirmative answer 
was obtained by Ocneanu, but he did not publish his proof. The only written proof I know 
of is due to Etingof, Nikshych, and Ostrik [30J, but it is more general and hence complex; 
in particular, it does not depend on the unitarity or the pivotal property. Meanwhile, the 
proof of Ocneanu rigidity for unitary theories is not very difficult and may be interesting to a 
mathematical physicist. 

To prove that a solution to an equation system is rigid, one needs to study infinitesimal 
deformations. Let us call a deformation of F, a, (3, R permissible if it satisfies the equations 
in the first order of Taylor expansion. The deformation is called trivial if it can be obtained 
by a change of basis, also in the first order. (By restricting our attention to the first order, we 
potentially make the space of permissible deformations larger and the space of trivial deforma- 
tion smaller.) Our goal is to calculate the quotient, permissible modulo trivial, and show that 
it vanishes. The analysis consists of several steps. 

First, we only consider F and the pentagon equation. The resulting deformation problem 
resembles the definition of H 3 (G,M.) for a group G. In fact, it is exactly that in example [2] on 
page 1721 Associated with an arbitrary fusion theory is a cochain complex defined by Crane and 
Yetter [HI] and independently by Davydov [H2] • It is a sequence of real spaces and linear maps 

C 1 C 2 ^ C 3 ^ - ■ ■ (238) 

whose third cohomology space H 3 (C) == Ker5 3 /Im<5 2 classifies the nontrivial deformations of 
the theory. (The first and the second cohomology spaces have natural interpretation too.) We 
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will use the term "tangent cohomologyT 4 ! proposed by Davydov. 

It is well known that any finite group has trivial real cohomology in all dimensions n ^ 1. 
Likewise, the tangent cohomology vanishes for an arbitrary unitary fusion theory (actually, any 
fusion theory). The proof of this statement is based on the same idea as the corresponding 
proof for groups. In the latter case, one uses averaging over the group, which is well-defined 
because the group is finite. For a unitary fusion theory, one needs to take the quantum trace 
(actually, a partial quantum trace) and average it over the label set with weight d a /T> 2 . 

The next step is to include a, (3 and to consider the triangle equations. This is just a 
technical detail that involves trivial changes to the deformation theory [93J. (The argument is 
quite general and applicable even if the label set is infinite, in which case the cohomology may 
not vanish.) The rigidity of braiding follows from the vanishing of H 2 (C) in conjunction with 
Vafa's theorem (see Theorem IE. 101) . 

E.6.1 Gauge freedom in the description of anyons. An isomorphism between two 
fusion theories, A and A' is given by a bijection between the label sets (we simply assume that 
they are equal) and a collection of unitary map@ : (V')^ b — > V^ b ; these data are enough 
to relate one system of associativity constraints to the other: 




{F'fu C - ^(/ — \y , Ft(j2 T ?® r u) = (E 1 ^® 1 ?)^? - (239) 



nafec 



For example, let A' be obtained from A by changing left and right, i.e., (V') a = V and 
^pi^abc _ ^pd>a^-x_ jf theory A has braiding R, then T ab = (R ab )~ l is an isomorphism (for a 
quick demonstration, apply the Yang-Baxter braid ( 12051) to both trees above). However, we will 
be mostly interested in the case where {V') ab = V^ b so that T ab may be called "basis change"; 
examples can be found in Sections 18.41 and [9j 

By analogy with Eq. (12391) . one can write similar equations for a a vs. a' a and (3 a vs. f3' a . 
However, the physical meaning of a and (3 is rather elusive, so it is not clear whether we obtain 
the most general form of equivalence this way. A reasonable notion of isomorphism between two 
theories, namely tensor functor arises naturally in the categorical formalism (see Sec. IE. 7.3[) . 
To state it in concrete terms, we just need to supplement the basis change with an overall phase 
factor 7 so that 

«. = 1^7 a'a, ft = r^ 7 /t (240) 

Finally, the braiding rules are changed as follows: 

RfTf = r b c a {R!)f. (241) 

Thus we have defined the family of theories A' that are isomorphic to a given theory A. 
However, their parametrization by T ab and 7 is redundant. In particular, F, a, (3, R do not 
change if we replace T ab and 7 with and <p such that 

$ f = A-If , (p = h l7 , (242) 



4 Etingof, Nikshych, and Ostrik call it "Yetter cohomology" in recognition of Yetter's further work in this 



area. 

25 



The direction of these maps is chosen to be consistent with the definition of a tensor functor (see Defini- 
tion HSU)- 
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where h x are arbitrary phase factors. In the categorical language, the tensor functors defined 
by (r, 7) and by ($, cp) are isomorphic. 

Thus we have entities of dimension 0, 1, and 2: theories of anyons, 1-isomorphisms (i.e., iso- 
morphisms between theories), and 2-isomorphisms (i.e., isomorphisms between 1-isomorphisms). 
Let us define three sets indexed by a complementary dimension: 

1. 2-automorphisms h of a given 1-isomorphism (r, 7); 

2. Equivalence classes of 1-automorphisms (r, 7) of a given theory A up to 2-isomorphisms; 

3. Equivalence classes of theories with given fusion multiplicities up to 1-isomorphisms. 

The first set is an Abelian group, the second is a group (see Appendix [Fj , the third does not 
have any special structure. 

E.6.2 Infinitesimal deformations and the vanishing of tangent cohomology. Let us 

define infinitesimal analogues of sets 1, 2, 3; we will eventually prove that they are trivial. The 
constructions involved are typical to cohomology theory, following a pattern that may already 
be recognized by studying set 1. We are eventually interested in set 3, which classifies solutions 
of the pentagon equation up to a basis change. In this subsection we ignore braiding and all 
attributes of the vacuum sector (i.e., 7, a, j3, and the triangle equations); these things will be 
considered later. 

For a fixed fusion theory with associativity constraints F£ bc , let 

h a « 1 - zX a , rf « id^ - iYf, (FX bc « (id v * - xZ* c ) , (243) 

where X a is an infinitely small real number and Y" c ab , Z® bc are infinitely small Hermitian oper- 
ators acting in Vf' and V£ bc , respectively!^! Now, let us substitute the expressions for h and T 
into Eq. ( 1242ft . assuming that $f = id V at. We get: 

yob = ( 6 1 X )>* def ^ _ Xc + Xa)id ^ ab (244) 

An infinitesimal analogue of set 1 is the space of deformations X = (X a : a £ M) such that 
S 1 X = 0. One can deal with Eq. ( 12391) and the pentagon equation in an analogous way. 

To present the results in a more convenient form, let us get rid of the lower index c. (In a 
categorical formulation of the theory, it does not appear at all.) Note that the set of operators 
yob £ Qr a rj c represents the action of a single element Y a ^ £ V® b , which may be constructed as 
follows: 

a b 

a b 



Y ab 



E\ldj b c uJ* ' where «v#K*» = ** ( 245 ) 



a b 



('A 



a6\f 



26 To be rigorous, Z^ c acts in the space @ e V^ b C$> V^ c , the representation of V£ bc corresponding to the {{ab)c) 
label grouping. But that is not so important as we always identify different groupings (= trees) using F. 
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(cf. Eg. I196p . For (t^Jf)" this procedure yields: 

a b a b a b 




(246) 

b 



More generally, we can use the following definition. 

Definition E.15. Let C n be the set of self-adjoint elements in @ an V^"^\ The tangent 
complex of the fusion theory is the sequence of R-linear maps 

n+l 

C° ^ C 1 ^ C 2 ^ C 3 ^ ■ ■ ■ , 5 n = ^(-l) fe /fc , (247) 

k=0 

where the maps fj}: C n — > C n+1 are defined as follows: 



11 0,2 «t!+l 



o\ a„ a n -\-i 













\ J v / a±...a n+ ± 




vX a2...a n+1 


(fn+l X ) ai ... an+1 ~ 


X a \ ...a n 


a 


l 


02 O n + 


L 


ai a„ 



(248) 



(f«x) = y 

\r K I ai...an+l /-^ 



da k d ak+1 





ffe - 


X...a k _ 


[Ca fc+2 ... 


1 





for k = 1, . . . , n. 



Ok a k+l 



(For n = we use = V± = C, hence C° = R. Correspondingly, /q , /i : 1 ^ ^ n '^„°' 
therefore 5° = / ° - /? = 0.) ^ 

Note that C is indeed a cochain complex, i.e., 5 n+1 5 n = 0, which follows from this easily 
verifiable identity: 

fT l fl = fl\\f: fi»0O<ma + l. (249) 

(It is part of structure that makes C into a cosimplicial space.) 

The reader may check that Eq. (12391) becomes Z = S 2 Y . Therefore the infinitesimal analogue 
of set 2 is given by solutions to the equation S 2 Y = modulo elements of the form Y = 5 X. 
Similarly, the pentagon equation may be written as 5 3 Z = 0, and the solutions should be 
considered modulo 5 2 Y. Thus, low-dimensional cohomology of the tangent complex has the 
following meaning: 

• H°{C) = R; 

• if 1 (C) classifies infinitesimal 2-automorphisms of the identity 1-automorphism; 
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• H 2 (C) classifies infinitesimal 1-automorphisms (i.e., basis changes that leave the associa- 
tivity constraints invariant) up to 2-isomorphisms; 

• H 3 (C) classifies infinitesimal deformations of the fusion theory up to 1-isomorphisms (i.e., 
arbitrary basis changes). 

Theorem E.16. H n (C) = for all n>0. 

Proof. We will use a standard method of proving vanishing cohomology results, namely con- 
tracting homotopy. Let 



X n : C n -> C 



n— 1 









X cai ...a n -i 









(250) 



a\ o, n -i 

We will show that 5\ + = 1 or ; more exactly, 

x n+l 6 n + § n-l x n = ^ for n>Q (251) 

If this is true, then any X £ Ker £ n can be represented as £ Im hence H n (C) = 0. 

Equation (12511) is an immediate corollary of this identity: 



,.n+l to 
A Jfc 



idc« if k — 0, 
/^X" iffc,n>0. 



(252) 



Let us rewrite it using graphic notation: 



fc = : 










^ai...a n 









ai...a n 5 



,a a i 



jfe = 1 



a,c,j 



dadai 



c\$j 

X 



0\ 



ca2---a„ 



£> 2 



"''] " " I 

ai ar, 



X 



ca2-..a„ 



oi 



(the case fc > 1 is trivial). The equation for k = is also obvious, whereas the one for k — 1 
follows from this identity: 

01 





a "i 

t = ^ 



c| 




3 /\ 3 

°1 Ol 




where we have used Lemma TE. 71 (Note that we did not use the pivotal property or the positivity 
of the inner product; this sheds some light on why the result holds in a more general setting |30j.) 

□ 
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E.6.3 Technicalities related to the unit and braiding. Let us now take into account 
additional structure that was neglected in the above analysis. We will get more variables 
and more equations, but the old equations will not change. The new degrees of freedom are 
characterized by infinitely small real numbers g, A a , B a and Hermitian operators W^ b which 
are defined as follows: 

7 « 1 - ig, a' a « (1 - iA a )a a , (3' a « (1 - iB a )f3 a , (253) 

{R')f « Rfiidy^-iWf). (254) 

(Recall that 7 is part of the definition of a 1-isomorphism.) It has been previously shown 
that any infinitesimal deformation of the associativity constraints F can be compensated by 
a suitable basis change V . Although a, /3, and R may still remain deformed, the problem is 
reduced to the case where F is fixed. Furthermore, any infinitesimal basis change not affecting 
F has the form (I242p . hence a, (3, or R are not altered either. The only parameter to tune is 
7 « 1 - ig. 

Let us first show that any permissible deformation of a and j3 is trivial, where "permissible" 
means satisfying the triangle equations and "trivial" means satisfying Eq. (12401) . Taken to an 
infinitesimal form, the first triangle equation (see Fig. [T9h) and Eq. ( 12401) read: 

A x = B w for all x, w; g + A a = g + B a = for all a. (255) 

Clearly, the first condition implies the second if we put g = —A\. 

Thus we may assume that the whole fusion theory, i.e., F, a, and (3 are fixed. It remains to 
show that the braiding deformation W vanishes, provided it satisfies an infinitesimal version of 
the hexagon equations. Instead of using the hexagon equations directly, we will rely on the fact 
that braiding defines an isomorphism between the given fusion theory A and the theory A 1 in 
which left and right are changed. Thus (Rf')~ 1 (R')f' is an automorphism of A. It follows that 
5 2 W = 0, hence 

Wf = (6 l X) a * = (X b -X c + X a )id v « b . 

for some X G C 1 . It is important that the right-hand side is a scalar times the identity, which 
enables us to calculate the deformation of R^R^ = j^- easily: 

° n (l-2i(X b -X c + X a )). 



0' a 9' b a a b 

But the topological spin is rigid due to Vafa's theorem, hence X b — X c + X a = 0. 



E.7 Categorical formalism (aside) 

Categories and functors are the language used by mathematicians to describe fusion, braiding, 
and related concepts. I originally tried to write an exposition of the theory of anyons using 
this formalism, but found it too awkward. I still think that functors are necessary for the 
understanding of phase transitions and other advanced properties of anyonic systems, but most 
things can be explained in more elementary terms. This section is a remainder of the abandoned 
plan. Please be warned that these notes are very incomplete, e.g., there is no discussion of 
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duality and related concepts: rigid, pivotal, and spherical categories (not to mention that we 
focus on semisimple categories — this restriction is natural for the intended applications). 

Let us outline the main elements of the theory. An abstraction called tensor category 
generalizes the notion of fusion theory. While anyonic fusion has a compact description in terms 
of basic data, a category is a huge collection of "objects" related by "morphisms". However, 
these relations form a regular structure that does not leave much freedom of choice. A rather 
trivial example is the category Hilb, whose objects are all possible finite-dimensional Hilbert 
spaces and the morphisms are all linear maps. A more interesting category Rep(G) is defined 
as follows: the objects are finite-dimensional unitary representations of a compact group G and 
the morphisms are intertwiners. We may also think of a fusion theory as a category: the objects 
are finite sequences of particle labels, and the morphisms between (61, ... , b m ) and (a\, . . . , ak) 
are fusion/splitting operators, i.e., elements of the space V^ 1 "^. Another way to turn a fusion 
theory into a category will be described later. 

By definition, a tensor category is equipped with an operation □ that is analogous to fusion. 
The role of the vacuum sector is played by a special object 1. In the categories Hilb and Rep(G), 
the operation □ is the usual tensor product and 1 is the set of complex numbers C regarded 
as a one-dimensional space or the trivial representation of G. 

Yet another example: the matrix category Mat mra . In this category!^! objects are m x m 
matrices whose entries are finite-dimensional complex linear spaces, morphisms between matri- 
ces are entrywise linear maps, and 

(A □ B) ft = A jk ® B kl , l jk =1 ^ = ^ (256) 

k I 1 ^ ' 

(If the matrix entries are finite-dimensional Hilbert spaces, we use the notation Mat^ m m because 
in this case each morphism has a Hermitian adjoint.) Note that (A □ B) □ C — A □ (B □ C) 
and AD1 — A— IDA; tensor categories with this property are called strict. In general, 
the equalities are replaced by isomorphisms F, a, and f3 satisfying the pentagon and triangle 
equations. 

A tensor functor is a map from one tensor category to another, a classic example being 
the embedding Rep(G) — > Hilb. Tensor functors have some physical applications, such as 
transformations of particles by global symmetries (see Appendix [F]). They are also related 
to the gauge freedom in the description of anyonic fusion, cf. Eqs. (12391) . (I240p . and (I242p . 
However, our main goal is to understand the status of the isomorphisms F, a, (3. To this end, 
we will find an embedding of an arbitrary tensor category into a strict one by a tensor functor 
that preserves distinction between morphisms (such functors are called faithful). In particular, 
any unitary fusion theory with label set M embeds into the matrix category Ma$ M M . The 
embedding theorem implies 

MacLane's coherence theorem: All morphisms composed of F, a, (3 and having the same 
source and target are equal. 

Indeed, in the strict category the source and the target are the same object, and any composition 
of F, a, and (3 is the identity morphism (i.e., equality). Thus all such compositions in the 

27 The more scientific name for this creature is C'"-mod-C™ 1 , the category of finite-dimensional bimodules over 
the algebra C m (the product of m copies of C). 
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original category are mapped to the same morphism. Since the functor preserves distinction 
between morphisms, the original morphisms are also equal. 

E.7.1 The basics of category theory. A category is a collection of objects and mor- 
phisms. The set of morphisms from A to B in a category C is denoted by C(A, B). 

Morphisms are anything that can be composed: if / : A — > B and g : B — > C, then 
gf : A — > C. It is only required that (hg)f = h(gf) and that every object A has an identity 
morphism id^ such that /id^ = ids/ = / for any f : A —> B. (Example: objects are vertices of 
a given graph, morphisms are paths of arbitrary length, identity morphisms are paths of length 

o.) 

A morphism / : A — > B is called an isomorphism if there is a morphism / 1 : B — > A such 
that = id^ and = id B . Notation for isomorphic objects: A = B. 

For applications in quantum mechanics, the basic example is the category Hilb of finite- 
dimensional Hilbert spaces. Among its objects are the set of complex numbers C and the one- 
dimensional space corresponding to the first level of a harmonic oscillator. These two spaces 
are isomorphic but different. (The isomorphism maps the complex number 1 to a normalized 
wavefunction if), but this map is not canonical because there is no reason to prefer if) over —if).) 
So, it does not generally make sense to ask whether two given objects are equal (unless they are 
the same by definition). Morphisms between two given objects may be compared for equality 
though. For two objects, X and Y, a reasonable question is whether X = Y. Of course, two 
spaces are isomorphic if and only if they have the same dimension. But we also want to keep 
track of isomorphisms, for they may not commute: a sequence of isomorphisms may result in 
a nontrivial automorphism u : X — > X, u ^ idx- (For Hilbert spaces, it is natural to consider 
unitary isomorphisms.) 

The next example has motivation in superselection theory. Let us consider an infinite 
spin system which is almost in the ground state, with excitations being allowed only in some 
finite region. Quantum states of such a system are classified by superselection sectors. For 
each sector a the state belongs to some finite-dimensional Hilbert space X a , which depends on 
specific constraints on the excitations. We are going to consider different sets of constraints 
(called "objects") and transformations from one object to another (called "morphisms"). 

Definition E.17. Let M be some set. The category Vec* M is defined as follows. 

• An object X G Vec M is a collection of finite-dimensional Hilbert spaces (X a : a G M), of 
which only finitely many are nonzero. 

• A morphism f : X — > Y is a collection of linear maps f a : X a — > Y a (for each a G M). 
The identity morphism idx : X — > X consists of the unit operators acting in X a . 

• Additional structure: The set Vec' M (X, Y) of morphisms from X to Y is a complex linear 
space. For any / : X — > Y we define the adjoint morphism y : Y — > X such that 

(P)a = (/«)*■ " 

A unitary isomorphism is an isomorphism / such that f^ 1 = /L 

An simplified version of this construction is the category VecM' we use complex linear spaces 
instead of Hilbert spaces and do not consider adjoint morphisms. 



96 



For each a G M we define the object [a] such that [a] a = C and [a] b = for b ^ a. Objects 
isomorphic to [a] are called simple. For any object Y, 

Y a = Vecl([a],Y). (257) 

From the physical perspective, simple objects correspond to quantum states (e.g., a particle 
pinned to a point), whereas [a] is a reference state in the given superselection sector. For 
a more formal example, consider the category Rep(G) whose objects are finite-dimensional 
unitary representations of a compact group G and whose morphisms are intertwiners. Let M 
be the set of irreps considered up to isomorphism, and let [a] be a particular irrep in the given 
isomorphism class. Then we may identify the category A = Rep(G) with Vexy M this way: 

Y a d ^A([a],Y). 

Remark E.18. Our definition of Vccm and Vec M resembles the construction of a vector space 
using coordinates. An invariant characterization is also possible, though somewhat complicated. 
When the label set M is unspecified, Vecu is called a semisimple C-linear category, and Vec* M 
is called a unitary category. 

The concepts of functor and natural morphism are central in category theory. They are ex- 
tremely general and therefore hard to grasp. We will try to illustrate them by simple examples, 
which are still rather abstract. More meaningful (though less direct) examples can be found in 
the next subsection. 

Definition E.19. A functor JF : A — > B maps each object A of the category A to an object 
J 7 (A) of the category B and each morphism / : A — > A to a morphism J-'(f) : J 7 {A) — > J 7 {A') 
such that 

Hhh) = Hh) Hh), ^(icU) = id^). 

A morphism (or natural transformation) between two functors J 7 , Q : A — > B is a collection ^1 
of morphisms : .F(vl) — > Q(A) such that for any / : A — > A' this diagram commutes 



'a' 



(251 



g(A)-^Lg(A') 



We have already mentioned one example of a functor, namely the embedding JF : Rep(G) — > 
Hilb. Acting on morphisms, it maps the set of intertwiners to the set of all linear maps between 
two representation spaces. A natural morphism h : T ^ J 7 may be constructed as follows: 
is the action of some fixed group element on the representation space of A. Indeed, for this 
particular choice of the functors T = Q the commutative diagram (I258P simply says that T(f) 
is an intertwiner. 

Example E.20. Let H be a finite-dimensional Hilbert space. We define a unitary functor 
7i : Hilb — > Hilb by tensoring with H on the left (notation: Ti = [if®]): 

H(A) = H®A, H(f) = id H ® /. 
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A linear map u : H — > G between two spaces defines a natural morphism U = [it ®] between 
the corresponding functors: 

Ua = w <8> idA- 

The functor 7i = [if®] has some special properties, namely the map / i— > Tt(f) is linear 
and TC(P) = (H(f)Y- Such functors are called unitary. 

Proposition E.21. 

1. Any unitary functor JF : Hilb — > Hilb is naturally isomorphic to the left-tensoring functor 
[H ®] for some space H , namely H = JF(C) . This isomorphism is unitary. 

2. Any natural morphism U : [£?<&]—>•[# <8>] /ias t/je form [u <8>], where u = Uc- 

Proof. The idea is very simple: any space is a direct sum of one-dimensional spaces, any one- 
dimensional space is isomorphic to C, therefore T and U are completely characterized by their 
action on C. Let us go through the detail to see how the formalism works. 

1. Let A be an arbitrary finite-dimensional Hilbert space. Elements of A may be associated 
with morphisms C — > A. An orthonormal basis corresponds to a set of morphisms ipj : C —> A 
such that iplipj = ^fcjidc and Yli^iYj = i°U- Since T is a unitary functor, 

j 

The same is true for the functor H = [J-(C) <£>]. We define a morphism ft a : T{A) — > W(A) as 
follows: 

where h c : JF(C) — > H(C) = JF(C) (g) C is an equality (indeed, X = X £g> C for any space X). 
It is obvious that ft^/iA = id^-M) and ft^A = id-^(A), hence ft. a is a unitary isomorphism. To 
show that ft is natural, consider another object A' with an orthonormal basis {ip'j} and a linear 
map / = J2j k c jki J 'ji J k' where G C. It is easy to check that 

h A >Hf) = ^cjkHi^hc^kY = H(f)h A . 

2. Let Q = [G<2>] and 7i = [H (&]. For an arbitrary object A with an orthonormal basis 
{ipj} we have U A Q(ipj) = 7i(ipj) Uc (due to the naturality of U). Hence 

U A = U A id g{A) = U A = J2 H ^ Ucg ^ = ^^idA- 

□ 

Adding another level of abstraction, we can reformulate Proposition IE.21I as follows: The 
category Fun* 1 (Hilb , Hilb) of unitary functors from Hilb to Hilb is isomorphic to the category 
Hilb. (One may replace Hilb by the category of finite-dimensional complex linear spaces, omit- 
ting the unitarity condition.) By analogy, we obtain the following result. 
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Proposition E.22. The category Fun^(Vec* M , Vec' N ) is isomorphic to the category Mat^ N M 
whose objects are N x M matrices of finite- dimensional Hilbert spaces with a finite number of 
nonzero entries in each column and whose morphisms are entrywise linear maps between such 
matrices. (If M is finite, then Mat NM = Vec NxM .) 

E.7.2 Fusion in categorical terms. Fusion theory may be formulated as additional struc- 
ture on the category C = Vec M . Let us assume that the the label set M contains a special 
element 1 and that V c ab , Ff c , a a , f3 a are defined and satisfy the pentagon and triangle equations 
as well as unitarity conditions. If in addition M is finite and the duality axiom on page ITU 
holds, then the resulting construction is called unitary fusion category. (For a general definition 
of fusion category see Ref. [30J . ) 

The use of categorical formalism has mathematical as well as physical motivation. Let A 
and B be spatially confined excitations such that A is located on the left of B. Each excitation 
may be described by an object in the category C. If we consider both excitations together, we 
will obtain a new object, A □ B, which may be called "physical tensor product". If A and B 
do not have local degrees of freedom, i.e., A = [a], B = [b], then (A □ B) c = V c . A nice 
property of the physical tensor product is that (A □ B) □ C — A □ (B □ C) (provided A, B, 
and C are arranged on the line in that particular order). Mathematically, the operation □ is 
only associative up to a canonical isomorphism, but we will find some abstract representation 
of objects in which the isomorphism becomes an equality. 

The tensor product of two objects, A, B £ C is defined by the equation 

(A □ B) c = V c ab ®A a ® B b . (259) 

a.b 

This operation is neither commutative nor a priory associative. Different ways to multiply 
several objects may be described by parenthesis structures, or trees. For example, 

a b c 

((H □ lb}) □ [c]) u = ©K a6 ®C 

e 
u 

The operation □ is a functor, meaning that for any morphisms / : A —>■ A' and g : B — » B' 
there is a morphism / □ g : A □ B —>■ A' □ B' such that 

(Mi) □ (9291) = if 2 □ □ 9i), icU □ ids = id AaB - (260) 

(The tensor product of morphisms is defined in the obvious way: (fOg) c = ^ a b idyab ® f a ®gb-) 
We may define a unitary isomorphism between (AO B) □ C and A □ (B □ C): 

F ABtC : (A □ B) □ C - A □ (B □ C), (F a ,b,c) u = £ ® id Aa ® id Bb ® id Cc , (261) 

a,b,c 

where F£ bc is the associativity map (11771) . Similarly, 

a A :A^Aa[l], p A : A -»• [1] □ A (262) 




99 



(the definition is obvious). It is important that the isomorphisms Fa,b,Ci a A, Pa are natural, 
i.e., for any morphisms / : A —>■ A', g : B — ► B', and h : C — > C we have 



oca'S = (S □ id[i]) « A , /W = (id[i] □ /) 



(263) 



Proposition E.23. The functor property ^26(A) and the naturality $263\) imply all quadrilateral 
identities of the tree calculus. 

Indeed, there are two kinds of such identities. In the example shown in Fig. [TTb . one of the 
commuting moves occurs inside a subtree descending from the other. The inner move plays the 
role of / in the naturality condition. If both moves occur in disjoint subtrees, one should use 
the functor property. 

E.7.3 Main definitions. Among things we are going to define, tensor category is just a 
general formulation of properties of the operation □. Tensor functor is a new concept though. 
For an elementary but important example, let the label types and fusion spaces be fixed but the 
associativity constraints vary. Given two sets of associativity constraints, we may ask whether 
they are equivalent up to a basis change, cf. Eqs. f!239j) and (1240)) . This type of equivalence 
is a special case of tensor functor. Note that two basis transformations may differ by trivial 
factors, cf. Eq. (j242p ; the corresponding tensor functors are said to be naturally isomorphic. 



Definition E.24. A tensor category C is a category C endowed with a functor □ : C x C — > C, 
a special object 1 e C, and natural isomorphisms 



F a ,b,c ■ (A □ B) □ C -»• A □ (B □ C), 



a A : A^ Ani, (3 A : A-^1UA (264) 



such that the following diagrams (variants of the pentagon equation in Fig. [T7k and the triangle 
equation in Fig. [T9k ) commute: 



((X UY)UZ)UW FxaY ' z ' w , (X □ Y) □ (Z □ W) 



(265) 



xn(Ya(za w)) 



{X □ (Y □ Z)) □ W 



Fx; 



Xn((YDZ) UW) 



xaw 



(xoi)nw Fx,1 - W > XD(1DW) 



(266) 



The tensor category C is called strict if Fa,b,c, a A, Pa are equalities. It is called semisimple 
if the base category C is Vec M . It is called unitary if C = Vec\ 4 , the functor □ is unitary, and 
the isomorphisms F, a, ft are unitary. 
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Semisimple tensor categories admit a compact description in terms of "basic data", which 
is just slightly more general than that of a fusion theory. Indeed, by analogy with Proposi- 
tions [K2U and [R22] one can show that any functor of two arguments, □ : Vecu x VecM — * VecM 
is given by a set of finite-dimensional spaces V® b . However, the unit object is not necessarily 



simple, an example being the category MatM,M or Mat M M with the multiplication rule and the 
unit ( I256p . (Note that these tensor categories are strict.) In physical terms, the nonsimplicity 
of the unit means that the vacuum is not unique; thus simple objects in Mat^ M M are boundaries 
between vacua. 

Definition E.25. Let C = (C, □, 1, F, a, (3) and & = (C, □', 1', F', a', 0') be tensor categories. 
A tensor functor Q : C — > C is a functor Q between the corresponding base categories, plus two 
natural isomorphisms 1^1 



such that the the following diagrams commute: 



7:1'- G(l) 



g((XDY)nz) 

G{Fx,y,z) 

g(xn(yuz)) 



■ g(xoY) a'g(z) 



g{x) ag(Yuz) 



id g(x) n ' r Y,z 



- e(x),g{Y),g(z) 

g{x) □' (g(Y) ag{z)) 



(267) 



(268) 



S(a x ) 

Q{X □ l) 



'g(x) 



^g(x) a v 

id S(X) n '7 

g(x) n>g(i) 



Q(Px) 

g(iax) 



°g(x) 



-l' n'g(x) 

7 n 'idg(x) 

g(i) wg{x) 



(269) 



The tensor functor g = (g,T, 7) is called unitary if all three of its components are unitary. 

This definition seems complicated, but basically it says the following: we may identify 
g(A □ B) with g(A) □' g(B) and with 1' (by means of T and 7) so that 



G{F 



X,Y,Z) 



r g(x),g{Y),g(z): 



Q{oc x ) 



a 



g(x)> 



G{Px) = P'c 



g(xy 



Remark E.26. A lax tensor functor is defined likewise, but Ta,b and 7 are not required to be 
isomorphisms. In physics, this construction describes a situation in which theory C is obtained 
from C as a result of Bose-condensation, cf. note M on page | 



Definition E.27. A morphism between two tensor functors, g = (g,T,j) and Ti. = (H,,Q,(p) 
is a natural morphism h : g — > 7i satisfying these additional commutation relations: 



g(X □ Y) 

hxnY 

H(X □ Y) 



-g(x) wg(Y) 

h x O'h Y 

H(X) u"H{Y) 



H(l) 



1' 



(270) 



2S 



The direction of arrows is a convention, which is slightly inconvenient for our purposes. 
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The morphism h is called a (unitary) isomorphism if all the maps hx are (unitary) isomor- 
phisms. 

Let us see how these definitions work for fusion theories. Let C and C have label sets M and 
M', respectively. We already know that an arbitrary functor Q : C — *■ C is described by some 
matrix (G°, : a G M, a' G M') of finite-dimensional Hilbert spaces. Then the supplementary 
components of the tensor functor, V and 7 are characterized by some unitary linear maps 

r"" : {V')f ® G a a , ® G$, - ® G£, 7 : C - G},. (271) 

Note that since the original 7 (the map form V to £7(1)) is an isomorphism, G\, = if a' 7^ 1'. 
Eqs. (12681) and (I269P impose certain algebraic relations on and 7. A unitary isomorphism 
between two tensor functors is described by a set of unitary maps h a a , : G a a , — > satisfying 
the equations that follow from Eq. (12701) . 

Returning to the example of basis transformations, let C = C, □ = □', but possibly 
F 7^ F', a 1 7^ a, (3' 7^ (3. Thus we fix the label types and fusion spaces while allowing different 
associativity relations. Let us further assume that Q is the identity functor, i.e., G a a = C and 
G a a , = if a' 7^ a. In this case, is just a unitary operator acting in the space Vf 3 = {V')f ) , 
and 7 is a phase factor. Then Eqs. (12681) and (I269P become (12391) and (I240p . respectively. In this 
setting, an isomorphism between two tensor functors is described by unital complex numbers 
h a satisfying Eq. (I242j) . 

E.7.4 The embedding theorem. Now, we will see the "abstract nonsense" in action. 
While we manipulate with definitions, some combinatorial magic happens behind the scenes. 
Basically, we turn pentagons into rectangles; this process may be viewed as a proof that the 
2-skeleton of the Stasheff polytope is simply connected. 

Theorem E.28. For any tensor category C = (C, □, l,F,a,(3) there exists a strict tensor cat- 
egory C' = (C, □', 1') and a faithful tensor functor Q = (G, T, 7) from C to C' (where "faithful" 
means that the images of distinct morphisms are also distinct). 

Proof (borrowed from Ref. [9l]). Let C be the category of functors from C to itself, □' the 
operation of composing two functors (usually denoted by o), and 1' the identity functor. We 

assume that functors act on the left, i.e., (X o Y)(W) = f X(Y(W)). For any object X we set 

G(X) = [Xn], i.e., g(X){W) = X □ W. The isomorphism 

^x,y between two functors may 

also be defined in terms of its action on a test object W: 

Y xx : [X □] o [Y □] -> [(X UY)U\, T XX {W) = F^ YW : X □ (Y □ W) -> (XDY) □ W. 

Finally, we define 

7 : id c -> [1 □], j(W) = @w ■ W->1DW. 

Now we simply rewrite Eqs. (I268j) and (I269p . applying them to XV: 

((xaY)az)aw FxnY - z ' n \ (XDY)n(zn w) x □ (Y □ (Z □ W)) 

{xn(Ynz))nw Fx < Ynz ' w \ xu((yuz)uw) idxaFY < z ' w , x □ (y □ (z □ w)) 
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XDW 

a x n id w 

(Xni)nw ■ 



= XD w 

id x D(3 w 

■XD(1DW) 



XDW 

/3 x oid w 



XDW 

Pxnw 



{iox)ow Fl,x,w ; in (xnw) 

These are the conditions we need to check. But the first and the second of them are identical to 
Eqs. (12641) and ( 12651) . respectively. The third condition is a different type of triangle equation, 
which is shown in Fig. [T9b . It follows from the standard one by Lemma [E.21 

To see that the functor Q is faithful, we set W — 1 and observe that a morphism h : X — ► Y 
is mapped to Q (h)(1) = hDid 1 : X □ 1 Y □ 1. Thus, 



a. 



-} (0(h)(1)) a x = h 



(due to the naturality of a). It follows that distinct morphisms remain distinct. 



□ 



Note that if C = Vec^, then C 



Mat 



t 

M,M 



. The functor Q takes a simple object a to the 



matrix whose [b, c] entry is V^ b . 

E.7.5 Braiding. Let us be brief and just give some definitions. 

Definition E.29. A braiding in a tensor category C = (C, □, l,F,a,/3) is a collection of iso- 
morphisms 

R a ,b ■ An B -> B □ A (272) 

which are natural with respect to A and B, i.e., for any morphisms / : A — ► A', g : B — >• B' 
this square commutes: 

/□s 



An B ■ 

Ra,b 



■A' OB' 



BOA - JLJ + B' □ A' 
It is also required that the following diagrams are commutative: 

Fy,x,z 



(YOX)OZ- 



■YO(XO Z) 



Rx,Y n id z 



(XoY)az 



idyORxz 



Y □ (ZUX) 



Fx; 



XD(YDZ) Rx ' Ynz , (YDZ)nX 



(273) 



(274) 



(ydx)dz Fy,x ' z > Yn(xnz) 



idyOR„\ 



(xnY)nz 



Yn(znx) 



Fx,Y,Z 

XD(YDZ) 



D— 1 

n YUZ,X 



Fy,z,x 

(Ynz)nx 



(275) 
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Definition E.30. A braided tensor functor is a tensor functor that commutes with braiding. 
More exactly, let C = (C, □ , 1, F, a, (3) and C = (C, □', 1', F', a', (3') be tensor categories fur- 
nished with braidings R and R', respectively. A tensor functor Q = (Q, F, 7) from C to C is 
called braided if the following diagram commutes: 



Q{A □ B) 

G{Ra,b) 

Q(B □ A) 



G{A) U'Q{B) 



R 'g(A),g(B) 

G{B) n>g(A) 



(276) 



Example E.31. Let C be an arbitrary braided tensor category and let C be defined as follows: 



C' = C, AW B = B U A, l' = l, 



771/ 

r X,Y,Z 



r Z,Y,Xi 



a 



ft, ft' = a, R'j 



A,B 



R 



B,A- 



Then the identity functor Q = idc supplemented with V a,b = Rab an d 7 — idi is a braided 
tensor functor from C to C. (We omit the proof.) 



Appendix F: Weak symmetry breaking 

A complete, yet to be discovered description of topological quantum order must include fusion 
and braiding rules (characterized by a unitary braided fusion category (UBFC), see Appendix [El) 
as well as symmetries of the underlying microscopic Hamiltonian. The no-symmetry case was 
considered in Sec. IE. 51 On the other hand, in our study of Abelian anyons in the honeycomb 
lattice model (see Sec. 15. 2\\ we observed a rather strange property: the translation by a lat- 
tice vector interchanges two superselection sectors. This phenomenon may be called a "weak 
breaking" of the translational symmetry. 

Let us now consider an arbitrary two-dimensional many-body system with a symmetry 
group G. For simplicity, we assume that the system consists of spins (or other bosonic degrees of 
freedom) rather than fermions. There are several ways in which the ground state can break the 
symmetry that is present in the Hamiltonian. First, the symmetry may be spontaneously broken 
in the usual sense, i.e., there may be a local order parameter. If such an order parameter does 
not exist, we may look for finer signs of symmetry breaking, in particular, for nontrivial action 
of G on superselection sectors. Such an action is described by a homomorphism u>i : G — > Ti, 
where Ti is some finite group defined below. 

But even if G does not permute the superselection sectors, there is still possibility for subtle 
symmetry-breaking properties, some of which were studied by Wen [951 EH] under the rubric 
of "projective symmetry groups" (PSG). While Wen's approach is rather indirect, one may 
instead investigate the action of symmetry transformations on superselection sectors. For a 
simple example, let us consider the Hamiltonian with J m < (note the analogy with Wen's 
model |96j). The ground state of this Hamiltonian may be obtained from the J m > state 
by putting an m-particle on every plaquette (so that any missing m-particle would now be 
regarded as an excitation). Therefore an e-particle picks up the phase factor —1 when it winds 
around a plaquette. It follows that the translations by two basis lattice vectors commute up to 
a minus sign when we consider their action on the nontrivial superselection sector e. 
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Effective dimension 
in which the sym- 
metry is broken 


Mathematical descrintion 


Examnles 





A local order parameter taking 
values in the coset space G/H, 
wneie n is a suogioup oi Lt 


Ferromagnet; Neel phase 


1 


Nontrivial action of the symme- 
try group on superselection sec- 
tors, which is characterized by a 
homomorphism uj\ : G — > Fx 


Abelian phases A x , A y , A z in the 
honeycomb lattice model 




The action of G on each super- 
selection sectoi ti is uesciiDeci oy 
a central extension G a ; the whole 
set of extensions is characterized 

Kv an plpmpnt iii„ FT^( 


Hamiltonian (I3j) with J rn < 0; 
<i nuuiDei oi moueis m [yo, yoj 


3 






4 


Effective topological action for a 
G-gauge field; such actions are 
classified by H A (BG,Z) 


Integer quantum Hall effect 



Table 6: Four levels of symmetry breaking in a two-dimensional quantum system. 



Mathematically, the up-to-sign commutation means that the translational group G = Z x Z 
is replaced by a central extension G e — > G with kernel Z 2 . If general, central extensions 
are classified by the cohomology group H 2 (G, U(l)). For each superselection sector a we get 
an extension G a . The whole collection of extensions is characterized by a cohomology class 
u>2 E H 2 (G,F 2 ), where T 2 is some finite Abelian group (see below). 

However, this is not the end of the story. Even without anyonic excitations (or if anyons 
exist but do not complicate the group action in any way), the system may have nontrivial 
properties such as the integer quantum Hall effect. The latter is described by a Chern-Simons 
term in the effective action for the electromagnetic field. Dijkgraaf and Witten [97] showed 
that a general topological action for a gauge field in 2 + 1 dimensions is characterized by an 
element of the cohomology group H 4 (BG,/L), where BG is the classifying space of G. (The 
group G is assumed to be compact.) An analogue of the Chern-Simons action exists even for 
finite groups. In this case, BG = K(G,1) (where K(G,n) is the Eilenberg-MacLane space) 
and H A (BG, Z) = H 3 (G, U(l)) (where the first H refers to the topological cohomology and the 
second to the group cohomology). 

The four-level classification of symmetry-breaking phenomena is summarized in Table [6j It 
has the flavor of topological obstruction theory (see e.g. [98]). One may hypothesize that there 
is a somehow relevant topological space Y such that iri(Y) = Fx, ^(Y) = F 2 , 7r 3 (Y) = 0, 
vr 4 (y) = Z, and all higher homotopy groups vanish. We now argue that the analogy with 
obstruction theory has precise mathematical meaning in dimensions 1 and 2. 

Suppose that the symmetry remains unbroken in dimension 0, but can be removed (if so 
is desired) by introducing a small perturbation to the Hamiltonian. That is obviously possible 
if the system consists of spins, while the Z, 2 symmetry associated with fermions can never be 
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removed. The properties that survive the perturbation are described by a unitary modular 
category (UMC) A and a chiral central charge c_, but we ignore the latter. 

The action of the symmetry group G on superselection sectors can be described algebraically 
as well as topologically. To begin with algebra, consider the 2-groupoid T = Aut(^4). It has 
a single object, the morphisms are invertible unitary braided tensor functors A — > A, and the 
2-morphisms are natural isomorphisms between such functors@ If the functors are considered 
up to isomorphisms, we obtain the group Ti that was mentioned above. The Abelian group T 2 
consists of all automorphisms of the identity functor. Note that 1^ may act on T 2 in a nontrivial 
way. The 2-groupoid T is characterized by this action and a cohomology class h G H 3 (Ti, T 2 ). 

On the topological side, one can consider the classifying space X = BY, which is a topo- 
logical space with a basepoint. I claim that X is connected, iri(X) = Ti, 7r 2 (X) = T 2 , and 
all higher homotopy groups vanish. Furthermore, the action of G on superselection sectors is 
described by a continuous map / : BG — > X defined up to a homotopy (some explanation is 
given below). One may ask if / is homotopic to a constant map. The first obstruction to such 
a homotopy is given by a homotopy class of maps uo\ : BG — * K(Ti, 1); such classes are in 
one-to-one correspondence with group homomorphisms G — > Y\. If this obstruction vanishes, 
one can define u; 2 G H 2 (BG, T 2 ) = H 2 (G,T 2 ). If uj 2 = 0, then / is homotopic to identity, or 
equivalently, the action of G on A is trivial. 

Note that X can be represented as a fibration with base -^(r^ 1) and fiber K(T 2 ,2). Its 
structure is characterized by the element h G if 3 (Ti,r 2 ). If h ^ 0, then the fibration does 
not have a cross section. In this case, the homomorphism u\ : G — > Y\ must satisfy a certain 
constraint, namely the inverse image of h by u\ must vanish. Thus, some seemingly possible 
types of symmetry breaking in dimension 1 may be forbidden due to an obstruction in dimension 
2. 

In conclusion, I conjecture a "physical" interpretation of the space X and the map /. 
Specifically, X parametrizes quantum states in the universality class described by the unitary 
modular category A. Each point of X is associated with a set of states that differ from each 
other locally; thus all finite-energy excitations of a given gapped Hamiltonian are represented 
by the same point. The original unperturbed system corresponds to some xq G X. To define 
the map BG — > X, let us assume that G is a compact Lie group which acts on each spin 
independently according to a faithful representation G — ► U(m). The Hilbert space of the spin, 
C m may be embedded in a large space C M . Each embedding u has an associated state g(u) G X. 
The space E of embeddings is contractible in the limit M — > oo, therefore its quotient with 
respect to the natural G-action may be regarded as a model of BG. Furthermore, since xq is 
symmetric, the map g : E —>■ X factors through BG. Thus the map / : BG — > X is defined. 

Remark F.l. It is not clear how to define the "extended classifying space" Y that has non- 
trivial homotopy in dimension 4. One may conjecture that Y = X x K(Z, 4), but this definition 
seems contrived. It may well be the case that the very idea to classify the Chern-Simons action 
as symmetry breaking in dimension 4 is wrong. Indeed, the Chern-Simons action is not fully 
described by homotopy theory because the fundamental homology class of the source (physical) 
manifold is involved. 

29 As explained to me by Ezra Getzler, higher groupoids may be nicely defined as Kan complexes [99]; this 
approach also allows one to endow the morphism sets with topology. In the definition of T no topology is 
assumed. It would be interesting to see what changes (if anything) when a natural topology is included. 
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